本文主要是介绍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;
}
这篇关于Disgruntled Judge UVA - 12169的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!