siprng.go 4.2 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267
  1. // Copyright 2014 The Gogs Authors. All rights reserved.
  2. // Use of this source code is governed by a MIT-style
  3. // license that can be found in the LICENSE file.
  4. package captcha
  5. import (
  6. "crypto/rand"
  7. "encoding/binary"
  8. "io"
  9. "sync"
  10. )
  11. // siprng is PRNG based on SipHash-2-4.
  12. type siprng struct {
  13. mu sync.Mutex
  14. k0, k1, ctr uint64
  15. }
  16. // siphash implements SipHash-2-4, accepting a uint64 as a message.
  17. func siphash(k0, k1, m uint64) uint64 {
  18. // Initialization.
  19. v0 := k0 ^ 0x736f6d6570736575
  20. v1 := k1 ^ 0x646f72616e646f6d
  21. v2 := k0 ^ 0x6c7967656e657261
  22. v3 := k1 ^ 0x7465646279746573
  23. t := uint64(8) << 56
  24. // Compression.
  25. v3 ^= m
  26. // Round 1.
  27. v0 += v1
  28. v1 = v1<<13 | v1>>(64-13)
  29. v1 ^= v0
  30. v0 = v0<<32 | v0>>(64-32)
  31. v2 += v3
  32. v3 = v3<<16 | v3>>(64-16)
  33. v3 ^= v2
  34. v0 += v3
  35. v3 = v3<<21 | v3>>(64-21)
  36. v3 ^= v0
  37. v2 += v1
  38. v1 = v1<<17 | v1>>(64-17)
  39. v1 ^= v2
  40. v2 = v2<<32 | v2>>(64-32)
  41. // Round 2.
  42. v0 += v1
  43. v1 = v1<<13 | v1>>(64-13)
  44. v1 ^= v0
  45. v0 = v0<<32 | v0>>(64-32)
  46. v2 += v3
  47. v3 = v3<<16 | v3>>(64-16)
  48. v3 ^= v2
  49. v0 += v3
  50. v3 = v3<<21 | v3>>(64-21)
  51. v3 ^= v0
  52. v2 += v1
  53. v1 = v1<<17 | v1>>(64-17)
  54. v1 ^= v2
  55. v2 = v2<<32 | v2>>(64-32)
  56. v0 ^= m
  57. // Compress last block.
  58. v3 ^= t
  59. // Round 1.
  60. v0 += v1
  61. v1 = v1<<13 | v1>>(64-13)
  62. v1 ^= v0
  63. v0 = v0<<32 | v0>>(64-32)
  64. v2 += v3
  65. v3 = v3<<16 | v3>>(64-16)
  66. v3 ^= v2
  67. v0 += v3
  68. v3 = v3<<21 | v3>>(64-21)
  69. v3 ^= v0
  70. v2 += v1
  71. v1 = v1<<17 | v1>>(64-17)
  72. v1 ^= v2
  73. v2 = v2<<32 | v2>>(64-32)
  74. // Round 2.
  75. v0 += v1
  76. v1 = v1<<13 | v1>>(64-13)
  77. v1 ^= v0
  78. v0 = v0<<32 | v0>>(64-32)
  79. v2 += v3
  80. v3 = v3<<16 | v3>>(64-16)
  81. v3 ^= v2
  82. v0 += v3
  83. v3 = v3<<21 | v3>>(64-21)
  84. v3 ^= v0
  85. v2 += v1
  86. v1 = v1<<17 | v1>>(64-17)
  87. v1 ^= v2
  88. v2 = v2<<32 | v2>>(64-32)
  89. v0 ^= t
  90. // Finalization.
  91. v2 ^= 0xff
  92. // Round 1.
  93. v0 += v1
  94. v1 = v1<<13 | v1>>(64-13)
  95. v1 ^= v0
  96. v0 = v0<<32 | v0>>(64-32)
  97. v2 += v3
  98. v3 = v3<<16 | v3>>(64-16)
  99. v3 ^= v2
  100. v0 += v3
  101. v3 = v3<<21 | v3>>(64-21)
  102. v3 ^= v0
  103. v2 += v1
  104. v1 = v1<<17 | v1>>(64-17)
  105. v1 ^= v2
  106. v2 = v2<<32 | v2>>(64-32)
  107. // Round 2.
  108. v0 += v1
  109. v1 = v1<<13 | v1>>(64-13)
  110. v1 ^= v0
  111. v0 = v0<<32 | v0>>(64-32)
  112. v2 += v3
  113. v3 = v3<<16 | v3>>(64-16)
  114. v3 ^= v2
  115. v0 += v3
  116. v3 = v3<<21 | v3>>(64-21)
  117. v3 ^= v0
  118. v2 += v1
  119. v1 = v1<<17 | v1>>(64-17)
  120. v1 ^= v2
  121. v2 = v2<<32 | v2>>(64-32)
  122. // Round 3.
  123. v0 += v1
  124. v1 = v1<<13 | v1>>(64-13)
  125. v1 ^= v0
  126. v0 = v0<<32 | v0>>(64-32)
  127. v2 += v3
  128. v3 = v3<<16 | v3>>(64-16)
  129. v3 ^= v2
  130. v0 += v3
  131. v3 = v3<<21 | v3>>(64-21)
  132. v3 ^= v0
  133. v2 += v1
  134. v1 = v1<<17 | v1>>(64-17)
  135. v1 ^= v2
  136. v2 = v2<<32 | v2>>(64-32)
  137. // Round 4.
  138. v0 += v1
  139. v1 = v1<<13 | v1>>(64-13)
  140. v1 ^= v0
  141. v0 = v0<<32 | v0>>(64-32)
  142. v2 += v3
  143. v3 = v3<<16 | v3>>(64-16)
  144. v3 ^= v2
  145. v0 += v3
  146. v3 = v3<<21 | v3>>(64-21)
  147. v3 ^= v0
  148. v2 += v1
  149. v1 = v1<<17 | v1>>(64-17)
  150. v1 ^= v2
  151. v2 = v2<<32 | v2>>(64-32)
  152. return v0 ^ v1 ^ v2 ^ v3
  153. }
  154. // rekey sets a new PRNG key, which is read from crypto/rand.
  155. func (p *siprng) rekey() {
  156. var k [16]byte
  157. if _, err := io.ReadFull(rand.Reader, k[:]); err != nil {
  158. panic(err.Error())
  159. }
  160. p.k0 = binary.LittleEndian.Uint64(k[0:8])
  161. p.k1 = binary.LittleEndian.Uint64(k[8:16])
  162. p.ctr = 1
  163. }
  164. // Uint64 returns a new pseudorandom uint64.
  165. // It rekeys PRNG on the first call and every 64 MB of generated data.
  166. func (p *siprng) Uint64() uint64 {
  167. p.mu.Lock()
  168. if p.ctr == 0 || p.ctr > 8*1024*1024 {
  169. p.rekey()
  170. }
  171. v := siphash(p.k0, p.k1, p.ctr)
  172. p.ctr++
  173. p.mu.Unlock()
  174. return v
  175. }
  176. func (p *siprng) Int63() int64 {
  177. return int64(p.Uint64() & 0x7fffffffffffffff)
  178. }
  179. func (p *siprng) Uint32() uint32 {
  180. return uint32(p.Uint64())
  181. }
  182. func (p *siprng) Int31() int32 {
  183. return int32(p.Uint32() & 0x7fffffff)
  184. }
  185. func (p *siprng) Intn(n int) int {
  186. if n <= 0 {
  187. panic("invalid argument to Intn")
  188. }
  189. if n <= 1<<31-1 {
  190. return int(p.Int31n(int32(n)))
  191. }
  192. return int(p.Int63n(int64(n)))
  193. }
  194. func (p *siprng) Int63n(n int64) int64 {
  195. if n <= 0 {
  196. panic("invalid argument to Int63n")
  197. }
  198. max := int64((1 << 63) - 1 - (1<<63)%uint64(n))
  199. v := p.Int63()
  200. for v > max {
  201. v = p.Int63()
  202. }
  203. return v % n
  204. }
  205. func (p *siprng) Int31n(n int32) int32 {
  206. if n <= 0 {
  207. panic("invalid argument to Int31n")
  208. }
  209. max := int32((1 << 31) - 1 - (1<<31)%uint32(n))
  210. v := p.Int31()
  211. for v > max {
  212. v = p.Int31()
  213. }
  214. return v % n
  215. }
  216. func (p *siprng) Float64() float64 { return float64(p.Int63()) / (1 << 63) }