Online Judge | Problem Set | Authors | Online Contests | User | ||||||
---|---|---|---|---|---|---|---|---|---|---|
Web Board Home Page F.A.Qs Statistical Charts | Current Contest Past Contests Scheduled Contests Award Contest |
0MS 的程序,运用辗转相除法联立方程直接求解,最多一层循环#include "iostream.h" #include "math.h" void Get_Special_Solution(int a,int b,int N,int x[2]) { int x1[2]={0},x2[2]={0},temp,flag1=1,flag2=1; if(N<0) {a=-a;b=-b;N=-N;} if(N==0) {x[0]=-b,x[1]=a;return;} if(abs(a)<abs(b)) {temp=a;a=b;b=temp;} if(a>0 && b<0) {b=-b;flag2=-1;} else if(a<0 && b>0) {a=-a;flag1=-1;} else if(a<0 && b<0) {a=-a;b=-b;} int m=a,n=b,q1,q2,r1,r2,temp1,temp2; q1=a/b,r1=a%b;q2=n/r1;r2=b%r1; m=r1;n=r2; x1[0]=1, x1[1]=-q1; //r1=a-q1*b; x2[0]=-q2, x2[1]=q1*q2+1; //r2=b-q2*r1;两个递推式来实现辗转相除法的求解过程 if(!r2) {x[0]=flag1*x1[0];x[1]=flag2*x1[1];return;} while(r2!=1) { temp1=x1[0];temp2=x1[1]; x1[0]=x2[0];x1[1]=x2[1]; q2=m/n;r2=m%n; x2[0]=-q2*x2[0]+temp1; x2[1]=-q2*x2[1]+temp2; m=n;n=r2; } x[0]=flag1*x2[0];x[1]=flag2*x2[1]; x[0]*=N;x[1]*=N; } void Show_Equation_Solve(int a[]) { int x[2],y[2],z[2],k1,k2; Get_Special_Solution(28,-23,a[0]-a[1],x); while(x[0]>0 && x[1]>0) {x[0]-=23;x[1]-=28;} while(x[0]<=0||x[1]<=0) {x[0]+=23;x[1]+=28;} Get_Special_Solution(33,-28,a[1]-a[2],y); while(y[0]>0 && y[1]>0) {y[0]-=28;y[1]-=33;} while(y[0]<=0||y[1]<=0) {y[0]+=28;y[1]+=33;} if(x[0]!=y[1]) Get_Special_Solution(33,-23,x[0]-y[1],z); while(z[0]>0 && z[1]>0) {z[0]-=23;z[1]-=33;} while(z[0]<=0||z[1]<=0) {z[0]+=23;z[1]+=33;} a[0]=28*z[1]+x[1],a[1]=23*z[1]+x[0],a[2]=28*z[0]+y[0]; } void main() { int i=0,m,days; int a[4]; while(1) { cin>>a[0]>>a[1]>>a[2]>>a[3]; if(a[0]==a[1] && a[1]==a[2] && a[2]==a[3] && a[3]==-1) break; m=a[0]; Show_Equation_Solve(a); cout<<"Case "<<++i<<": the next triple peak occurs in "<<(a[0]*23+21252+m-a[3]-1)%21252+1<<" days."<<endl; } } Followed by:
Post your reply here: |
All Rights Reserved 2003-2013 Ying Fuchen,Xu Pengcheng,Xie Di
Any problem, Please Contact Administrator