解线性同余方程的应用(一)

 例题一 

题目如图所示



典型的解同余方程组的题目，代码如下

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
long long exgcd(long long a,long long b,long long &x,long long &y)
{
if(b==0)
{
x=1;
y=0;
return a;
}
long long ans=exgcd(b,a%b,x,y);
long long temp=x;
x=y;
y=temp-a/b*y;
return ans;
}
int main()
{
long long i,n,a1,r1,a2,r2,a,b,c,x0,y0;
while(scanf("%lld",&n)!=EOF)
{
bool ifhave=1;
scanf("%lld%lld",&a1,&r1);
for(i=1;i<n;i++)
{
scanf("%lld%lld",&a2,&r2);
a=a1,b=a2,c=r2-r1;
long long g=exgcd(a,b,x0,y0);
if(c%g!=0)
{
ifhave=0;
}
int t=b/g;
x0=(x0*(c/g)%t+t)%t;
r1=a1*x0+r1;
a1=a1*a2/g;
}
if(ifhave==0)
{
printf("-1\n");
continue;
}
printf("%lld\n",r1);
}
return 0;
}

例题2

c looooops(pku 2115)

对于循环语句

for(variable=A; variable!=B;variable+=C )

statement;

已知所有的数均为k进制无符号整数类型，给出A,B,C和k的值，计算并输出statement执行的次数， 如果执行的次数为无限次，那么直接输出"FOREVER"。

输入：输入数据有多组，每组数据占一行，有四个整数A,B,C,k（0<=A,B,C<2^k）,输入为0 0 0 0时程序结束。

输出：对于每组数据，输出执行的次数，如果为无限次，那么直接输出"FOREVER";

分析：

题意不难理解，只是利用了 k位存储系统 的数据特性进行循环。

例如int型是16位的，那么int能保存2^16个数据，即最大数为65535（本题默认为无符号），

当循环使得i超过65535时，则i会返回0重新开始计数

如i=65534，当i+=3时，i=1

其实就是 i=(65534+3)%(2^16)=1

一元同余方程求解，假设算法恰好执行了X步，则A+CX ≡B(mod M)，其中M=1<<k.将其变形后，用扩展欧几里得算法求解。

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
long long exgcd(long long a,long long b,long long &x,long long &y)
{
if(b==0)
{
x=1;
y=0;
return a;
}
long long ans=exgcd(b,a%b,x,y);
long long temp=x;
x=y;
y=temp-a/b*y;
return ans;
}
int main()
{
long long a,b,c,k;
long long x,y;
while(scanf("%lld%lld%lld%lld",&a,&b,&c,&k)!=EOF)

ae5c
{
if(a==0&&b==0&&c==0&&k==0)
break;
long long M=1<<k;
long long g=exgcd(c,M,x,y);///求未知数前的系数的最大公因数
if((b-a)%g!=0)
cout<<"FOREVER"<<endl;
else
{
long long ans=x*(b-a)/g;
long long temp=M/g;
ans=(ans%temp+temp)%temp;
cout<<ans<<endl;
}
}
return 0;
}