参考了 https://www.cnblogs.com/Przz/p/5409721.html 大佬的题解 他写的推导非常清晰
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod =10001;
long long f[10005];
void gcd(ll a,ll b,ll&d,ll&x,ll&y){//拓展欧几里得
if(!b){
d=a;
x=1,y=0;
return ;
}
else {
gcd(b,a%b,d,y,x);
y-=x*(a/b);
return ;
}
}
int main(){
int t;
while(cin>>t){
memset(f,0,sizeof(f));
for(int i = 1; i <= 2*t; i+=2)
scanf("%I64d",&f[i]);
for(int a=0;a<mod;a++){
ll k,b,d;
gcd(mod,a+1,d,k,b);
ll c=f[3] - a * a * f[1];//手推公式发现ax+by=c 这里的c是这个
if(c%d)continue;//不是倍数没有整数解
b=b*c/d;//算出解b进行下面的验证
bool flag=1;
for(int i=2;i<=2*t;i++){
int temp=(a*f[i-1]+b)%10001;
if(i%2){
if(temp!=f[i]){
flag=0;
break;
}
}
else f[i]=temp;
}
if(flag)break;
}
for(int i = 2; i <= t*2; i += 2) cout << f[i] << "\n";
}
return 0;
}