Uva 12169 Disgruntled Judge(枚举)

链接：https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&category=326&problem=2290&mosmsg=Submission+received+with+ID+2196098

大意：丧心病狂的出题人设计了一个长度为2t的数列x，满足x[i+1]=a*x[i]+b mod 10001，将x[1],x[3]....作为输入，x[2],x[4]....作为输出,其中a  b和x[i]均在[0,10001]中,给出输入,求输出.

分析：首先10001=73*137，可以枚举a，相应地算出b，b满足x[3]-a*a*x[1]=(a+1)*b mod 10001，用中国剩余定理算一下b，然后就可以推出所有的x[i]，检查下是否与输入矛盾，不矛盾的话就是一组解。

#include<iostream>
#include<cstring>
using namespace std;
typedef long long ll;
int x[2*100+5];
int inv(int a,int p){
    for(int i=0;i<p;i++){
        if((a*i)%p==1)return i;
    }
}
int main(){
    int t;
    cin>>t;
    memset(x,-1,sizeof(x));
    for(int i=0;i<t;i++){
        cin>>x[2*i+1];
    }
    for(int a=0;a<10001;a++){
        bool ok=true;
        int b=(1096*((x[3]-(ll)a*a*x[1])%10001)*inv(a+1,73))%10001+(8906*(((x[3]-(ll)a*a*x[1])%10001)*inv(a+1,137)))%10001;
        b%=10001;
        for(int i=2;i<=2*t;i++){
            int q=(a*x[i-1]+b)%10001;
            if(i%2&&x[i]!=q){
                ok=false;break;
            }
            x[i]=q;
        }
        if(ok){
            for(int i=2;i<=2*t;i+=2)
                cout<<x[i]<<endl;
        }
    }
    return 0;
}