POJ2115 扩展欧几里得
2013-01-07 16:31
330 查看
扩展欧几里得模板题 根据题意理出二元一次方程 A+X*C-B=Y*2^k 移项可得X*C+Y*2^k=B-A
然后扩展欧几里得 就得出答案了
#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
void exgcd(long long a,long long b,long long &d,long long &x,long long &y)
{
if(b==0)
{
x=1;
y=0;
d=a;
return;
}
else
{
exgcd(b,a%b,d,x,y);
long long temp=x;
x=y;
y=temp-(a/b)*y;
}
}
int main()
{
long long a,b,c,k,m[35],gc,x,y;
m[1]=2;
for(int i=2; i<=32; i++)
m[i]=2*m[i-1];
while(~scanf("%lld%lld%lld%lld",&a,&b,&c,&k)&&(a+b+c+k))
{
exgcd(c,m[k],gc,x,y);
if((b-a)%gc==0)
{
x*=((b-a)/gc);
x=(x%(m[k]/gc)+m[k]/gc)%(m[k]/gc);
printf("%lld\n",x);
}
else
printf("FOREVER\n");
}
return 0;
}
然后扩展欧几里得 就得出答案了
#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
void exgcd(long long a,long long b,long long &d,long long &x,long long &y)
{
if(b==0)
{
x=1;
y=0;
d=a;
return;
}
else
{
exgcd(b,a%b,d,x,y);
long long temp=x;
x=y;
y=temp-(a/b)*y;
}
}
int main()
{
long long a,b,c,k,m[35],gc,x,y;
m[1]=2;
for(int i=2; i<=32; i++)
m[i]=2*m[i-1];
while(~scanf("%lld%lld%lld%lld",&a,&b,&c,&k)&&(a+b+c+k))
{
exgcd(c,m[k],gc,x,y);
if((b-a)%gc==0)
{
x*=((b-a)/gc);
x=(x%(m[k]/gc)+m[k]/gc)%(m[k]/gc);
printf("%lld\n",x);
}
else
printf("FOREVER\n");
}
return 0;
}
相关文章推荐
- POJ2115 C Looooops 模线性方程(扩展欧几里得)
- POJ2115 C Looooops(扩展欧几里得)
- 【扩展欧几里得】poj2115 C Looooops
- POJ2115:C Looooops(扩展欧几里得)
- POJ2115(扩展欧几里得)
- poj2115 C Looooops(扩展欧几里得)
- POJ2115 - C Looooops(扩展欧几里得)
- POJ2115 C Looooops[扩展欧几里得]
- poj2115(扩展欧几里得运用)
- poj2115(扩展欧几里得)
- noip2012 同余方程 (扩展欧几里得)
- 费小或者扩展欧几里得求逆元
- POJ2142_The Balance_扩展欧几里得
- 扩展欧几里得
- 扩展欧几里得 【记录】
- POJ 1061 青蛙的约会(扩展欧几里得)
- poj1061青蛙的约会--扩展欧几里得
- codeforces 7c line(扩展欧几里得)
- 扩展欧几里得--让你一次刷个够
- 扩展欧几里得