Disgruntled Judge UVA - 12169

Disgruntled Judge UVA - 12169

题意：给定2*T大小的数组，给定奇数项，符合关系x[i] = (a · x[i−1] + b) mod 10001.求各个偶数项。

题解：范围10000，可以通过暴力求解过，但是时间复杂度太高，根据公式可以采用拓展欧几里得求解。

   可以得到  x2 = (a * x1 + b) % 10001;

x3 = (a * x2 + b) % 10001;

   联立得   x3 - a * a * x1 = (a + 1) * b + 10001 * (-k);

   然后枚举a的值,通过奇数项求得b，看看是否有冲突，无冲突即可求得答案。（注意使用long long）

AC代码：

#include<iostream>
#include<cstdio>
using namespace std;

const int maxn = 105;
const int MOD = 10001;

long long t,n;
long long a[maxn*2];
void gcd(long long a,long long b,long long &d,long long &x,long long &y)
{
    if(b==0)
    {
        x=1;
        d=a;
        y=0;
    }
    else
    {
        gcd(b,a%b,d,y,x);
        y-=x*(a/b);
        return ;
    }
}
int main()
{
    scanf("%d",&t);
    for(int i=1;i<=t;i++)
    {
        scanf("%d",&a[i*2-1]);
    }
    long long d,b,k;
    for(int i=0;;i++)
    {
        gcd(MOD,i+1,d,k,b);
        long long tt = a[3]-i*i*a[1];
        if(tt%d)continue;
        b=b*tt/d;
        int ok = 1;
        for(int j=2;j<=2*t;j++)
        {
            if(j&1)
            {
                if(a[j]!=(a[j-1]*i+b)%MOD)
                {
                    ok=0;
                    break;
                }
            }
            else
            {
                a[j]=(i*a[j-1]+b)%MOD;
            }
        }
        if(ok==1)break;
    }
    for(int i=2;i<=2*t;i+=2)
    {
        cout<<a[i]<<endl;
    }
    return 0;
}