google code jam 2010资格赛B题 Fair Warning

B题本身不是很难，最大公约数大家一般都能想到，关键是这里全是大整数。我昨天做题的时候临时写了一个大整数模板，为了方便我用不断做减法实现的求模运算，结果小数据还凑合，但大数据跑到第5组数据就卡着不动了。

今天早上起来看官方分析，很明显的表示以后会出大整数方面的题目，让大家做好准备。

根据提示找了到gmp类库，结果和VS兼容不好，搜到几个别人编译好的dll，但是好几年前的了，只能在vc6下。

不断找不断找，发现有个mpir类库，和gmp基本一致，但支持ms下的vs。

下面就是利用这个类库写的B题的代码，速度很快的说～～

#include "stdafx.h" #include <iostream> using namespace std; #include "mpir.h" #pragma comment(lib,"mpir.lib") mpz_t bignums[1000],zero; char largeNum[1000]; void GCD(mpz_t a,mpz_t b,mpz_t ans)//返回类型不能用 mpz_t 类型，所以我把常规的GCD运算修改了下 { if(mpz_cmp(b,zero)==0) { mpz_set(ans,a); return; } mpz_mod(a,a,b); GCD(b,a,ans); } int main() { freopen("e://B-large-practice.in","r",stdin); freopen("e://1.out","w",stdout); mpz_init(zero); mpz_set_ui(zero,0); int C,N,minIndex; scanf("%d",&C); for(int i=1;i<=C;i++) { scanf("%d",&N); for(int j=0;j<N;j++) { scanf("%s",largeNum); mpz_init(bignums[j]); mpz_set_str(bignums[j],largeNum,10); } minIndex=0; for(int j=1;j<N;j++) if(mpz_cmp(bignums[j],bignums[minIndex])<0) minIndex=j; for(int j=0;j<N;j++) { if(j!=minIndex) mpz_sub(bignums[j],bignums[j],bignums[minIndex]); } mpz_t offence; mpz_init(offence); mpz_set(offence,bignums[minIndex]); mpz_sub(bignums[minIndex],bignums[minIndex],bignums[minIndex]); mpz_t gcd; mpz_init(gcd); GCD(bignums[0],bignums[1],gcd); for(int j=2;j<N;j++) GCD(gcd,bignums[j],gcd); mpz_mod(offence,offence,gcd); mpz_t ans; mpz_init(ans); if(mpz_cmp(zero,offence)!=0) mpz_sub(ans,gcd,offence); gmp_printf("Case #%d: %Zd/n",i,ans); } return 0; }