Disgruntled Judge UVA - 12169 拓展欧几里得 +暴力

参考了 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;
}