Gym 100641A Continued Fractions 模拟|签到

题目链接：点我

题目大意：给你一个连分式的公式，然后给你两个实数的组成数字，求出其和 差 积 商 同样以连分式的形式表现出来。

题目思路：模拟即可

代码参考：来自队友@BinacsLee@Vectorhr

AC代码：

/*
2017年8月15日09:28:08
A题模拟 
WA:头文件改了一下就过了
RE：换了另外一个GCD模板 
AC
*/
#include <cstdio>
#include <algorithm>//很难受，加了这个头文件就过了 
#include <cmath>

using namespace std;
typedef long long ll;
const int maxn=20;
ll a[maxn],b[maxn];
int n1,n2;

ll gcd(ll a,ll b){  
    if(b) return gcd(b,a%b);
	return a;
}
void simple(ll &fz,ll &fm){
	ll g=gcd(fz,fm);
	fz/=g;
	fm/=g;
}

void add(ll fz1,ll fm1,ll fz2,ll fm2,ll &fz,ll &fm){
	ll g=gcd(fm1,fm2);
	fz=fm1/g*fz2+fm2/g*fz1;
	fm=fm1/g*fm2;
	simple(fz,fm);
}

void sub(ll fz1,ll fm1,ll fz2,ll fm2,ll &fz,ll &fm){
	ll g=gcd(fm1,fm2);

	fz=fm2/g*fz1-fm1/g*fz2;
	fm=fm1/g*fm2;
	simple(fz,fm);
}

void mul(ll fz1,ll fm1,ll fz2,ll fm2,ll &fz,ll &fm){
	
	ll g1=gcd(fm1,fz2);
	ll g2=gcd(fm2,fz1);
	if(g1!=1) {
		fm1/=g1;
		fz2/=g1;
	}
	if(g2!=1) {
		fm2/=g2;
		fz1/=g2;
	}
	
	fm=fm1*fm2;
	fz=fz1*fz2;
	simple(fz,fm);
	
}

void div(ll fz1,ll fm1,ll fz2,ll fm2,ll &fz,ll &fm){
	
	ll t=fz2;
	fz2=fm2;
	fm2=t;
	mul(fz1,fm1,fz2,fm2,fz,fm);	
} 

void print(ll fz,ll fm){
	ll t;
	if((fz<0&&fm>0)||(fz>0&&fm<0)){
		printf("%I64d",fz/fm-1);
		fm=abs(fm);
		fz=fm-abs(fz%fm);
	}
	else{
		printf("%I64d",fz/fm);
		fm=abs(fm);
		fz=abs(fz);
		fz=fz%fm;
	}
	while(fz!=1){
		printf(" %I64d",fm/fz);
		//ll t=fm-fz;
		t=fz;
		fz=fm%fz;
		fm=t;
	}
	if(fm) printf(" %I64d",fm);
	printf("\n");
}

int main(){
	int kase=0;
	while(scanf("%d%d",&n1,&n2)==2&&n1&&n2){
		
		for(int i=1;i<=n1;i++){
			scanf("%I64d",&a[i]);
		}
		
		for(int i=1;i<=n2;i++){
			scanf("%I64d",&b[i]);
		}
		
		ll fz1,fm1,fz2,fm2,fenzi,fenmu,tmp;
		if(n1==1) fenzi=1,fenmu=a[n1]; 
		else{
			fenzi=1;fenmu=a[n1];
			for(int i=n1-1;i>=1;i--){
				tmp=fenzi+a[i]*fenmu;
				fenzi=fenmu;
				fenmu=tmp;
			}
		}
		fz1=fenmu;
		fm1=fenzi;
		simple(fz1,fm1);
		
		//printf("r1:%d/%d \n",fz1,fm1);
		if(n2==1) fenzi=1,fenmu=b[n2]; 
		else{
			fenzi=1;fenmu=b[n2];
			for(int i=n2-1;i>=1;i--){
				ll tmp=fenzi+b[i]*fenmu;
				fenzi=fenmu;
				fenmu=tmp;
			}
		}
		fz2=fenmu;
		fm2=fenzi;
		simple(fz2,fm2);
		
		printf("Case %d:\n",++kase);
		
		ll rfz,rfm;
		add(fz1,fm1,fz2,fm2,rfz,rfm);
		if(rfz==0) printf("0\n");
		else{
			if(rfm==1) printf("%I64d\n",rfz);
			else if(rfm==-1) printf("%I64d\n",-rfz);
			else{
				print(rfz,rfm);
			}
		}
		
		sub(fz1,fm1,fz2,fm2,rfz,rfm);
		if(rfz==0) printf("0\n");
		else{
			if(rfm==1) printf("%I64d\n",rfz);
			else if(rfm==-1) printf("%I64d\n",-rfz);
			else{
				print(rfz,rfm);
			}
		}
		
		mul(fz1,fm1,fz2,fm2,rfz,rfm);
		if(rfz==0) printf("0\n");
		else{
			if(rfm==1) printf("%I64d\n",rfz);
			else if(rfm==-1) printf("%I64d\n",-rfz);
			else{
				print(rfz,rfm);
			}
		}
		
		div(fz1,fm1,fz2,fm2,rfz,rfm);
		if(rfz==0) printf("0\n");
		else{
			if(rfm==1) printf("%I64d\n",rfz);
			else if(rfm==-1) printf("%I64d\n",-rfz);
			else{
				print(rfz,rfm);
			}
		}
	//	printf("r2:%d/%d \n",fz2,fm2);
	
		
	}
	return 0;
}