BZOJ-3122-随机数生成器-SDOI2013-BSGS

描述

---

分析

• 关键就是对式子的变形, 将原递推公式转化为通项公式, 中间会用到等比数列求和公式. 然后切记此时方程两边同乘(1-a), 然后就化简开了.
• 之后得到的式子是 a^(n-1) = ... (mod p) 然后就可以用BSGS求解了. 注意最后加1.
• 公式明天补上
• 几种特判情况 :
  
  • 如果a = 0, 这时 X = b, 如果t = b输出1, 否则输出-1
  • 如果a = 1, 为什么这也要特判呢, 因为等比数列求和公式里的q一定不为1.
    此时Xi = Xi-1 + b是等差数列 Xn = X1 + (n-1) * b (mod p), i = (Xn - X1 + b) * b^(-1)
    这时还需要看如果b = 0, 输出-1; 如果Xn-X1+b = 0, 求出 i = 0, 但是想想实际情况应该是第 p 天才看到那一页
• 现在写代码有点不走脑子, 调了半天发现我图省事把变量全开long long结果就是有个地方开了int还没转换.
• 感觉最近睡的有点少, 该充电了.

---

代码

https://code.csdn.net/snippets/624452