题目大意:定义无穷数列{Xi}:
1.X0=A;
2.xi = (alpha * xi-1^2 + beta * xi-1 + gamma) mod M, i >= 1.
已知 A, alpha, beta, gamma, M, K.
求 Xk.
输入:A (1 <= A <= 10000), alpha (0 <= alpha <= 100), beta (0 <= beta <= 100), gamma (0 <= gamma <= 100), M (1 <= M <= 1000), k (0 <= k <= 10^9)。所有的数为整型。
输出:Xk
思路:一开始想推公式,推不出。发现M只有1000,明显是一个找循环节的问题。意思就是说在最多1000次运算的时候Xi就能重复。意会一下就出来了~
比较坑的地方是X0是不用模M的,半天才看出来。然后就很水了。。。。
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
int x0, a, b, c, m, k;
int flag[1111], x[1111];
int main() {
int i;
memset(flag, 128 / 3, sizeof(flag));
scanf("%d%d%d%d%d%d", &x0, &a, &b, &c, &m, &k);
if(k == 0){
printf("%d\n", x0);
return 0;
}
x[0] = x0 % m;
for(i = 0; flag[ x[i] ]>1111; i++) {
flag[ x[i] ] = i;
x[i+1] = (a * x[i] * x[i]+ b * x[i] + c) % m;
}
int ans= x[flag[ x[i] ] + (k - flag[ x[i] ]) % (i - flag[ x[i] ]) ];
if (k >= flag[ x[i] ]) printf("%d\n", ans);
else printf("%d\n", x[k]);
return 0;
}