pku 2115 C Looooops（解模线性方程）

原创于 2009-08-13 16:39:00 发布 · 789 阅读

1 ·

CC 4.0 BY-SA版权

文章标签：

#c #编译器

杂七杂八专栏收录该内容

17 篇文章

订阅专栏

本文介绍了一个关于64位整数左移操作的常见错误及其正确做法，并提供了一个使用64位整数解决线性方程问题的示例代码。文章详细解释了为何在某些情况下直接进行左移操作会导致问题，并给出了正确的解决方案。

__int64 t=1<<k;

本来t有64位，k可以取（0-62），但这里因为编译器会把常量1当成int处理，所以当k>30就会出现问题。

正确的写法是 __int64 t=(__int64)1<<k;

#include <iostream> using namespace std; __int64 exGcd(__int64 a,__int64 b,__int64& x,__int64& y) { if(b==0) { x=1; y=0; return a; } __int64 r=exGcd(b,a%b,x,y); __int64 t=x; x=y; y=t-a/b*y; return r; } void liner_equation_solver(__int64 a,__int64 b,__int64 n) { __int64 x,y,d; a=a%n; if(a<0) a+=n; b=b%n; if(b<0) b+=n; d=exGcd(a,n,x,y); if(b%d) printf("FOREVER/n"); else { x=x*b/d; x=x%n; if(x<0) x+=n; while(x>=n/d) { x-=n/d; } printf("%I64d/n",x); } } int main() { __int64 A,B,C,k,t; while(scanf("%I64d%I64d%I64d%I64d",&A,&B,&C,&k)) { if(A==0&&B==0&&C==0&&k==0) return 0; t=(__int64)1<<k; liner_equation_solver(C,B-A,t); } return 0; }