siprng.go 4.6 KB

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