Disgruntled Judge UVA - 12169 (暴力或者拓展欧几里得)

这里数据太水了吧，暴力竟然都可以过。
然后第二种方法是拓展欧几里得。
题意：给你T组数据 x[1],x[3],x[5]…x[2*T-1]

x[i]=(a*x[i-1]+b)%10005; 让你输出x[2],x[4],…x[2*T];

枚举a，然后通过x[1],x[3]求b，然后验证符不符合。

-10001k+(a+1)b=x3-a*a*x1; 就是拓展欧几里得，记得用gcd求出来的一组解还要乘上c/gcd；

这里b只用一组解就行了，因为b=y0+10001*k，之后会被约掉

wa，这里有一个坑点，就是会报PE，因为UVA 要用%lld,那为什么不报WA呢。。。。。。。。

#include <iostream>
#include <cstring>
#include <cmath>
#include <cstdio>
#include <algorithm>

using namespace std;
const int maxn = 20005;
#define inf 10001
#define ll long long
void gcd(ll a,ll b,ll &d,ll &x1,ll &y1)
{
    if(!b){
        x1=1;y1=0;d=a;
        return ;
    }
    else {
        gcd(b,a%b,d,y1,x1);
        y1-=x1*(a/b);
        return ;
    }
}
ll x[maxn];
int main()
{
    //freopen("E:\\input.txt","r",stdin);
        int t;
        scanf("%d",&t);
        for(int i=1;i<=t;i++)
            scanf("%I64d",&x[2*i-1]);
        int flag=0;
        for(ll a=0;a<=10000;a++)
        {
            ll temp=(x[3]-a*a*x[1]);
            ll x1,y1,d;
            gcd(10001,a+1,d,y1,x1);//这里要对应好
            if(temp%d) continue;
            else {
                ll k=x1*temp/d;//这里temp/b不能加abs
                for(int i=2;i<=2*t;i++)
                {
                    ll tt=(a*x[i-1]+k)%inf;
                    if(i&1)
                    {
                        if(tt!=x[i])
                            break;
                    }
                    x[i]=tt;
                    if(i==(2*t))
                        flag=1;
                }
                if(flag==1) break;
            }
            if(flag==1) break;
        }
        for(int i = 2 ; i <= 2 * t ; i+=2)
              printf("%lld\n",x[i]);
    return 0;
}