高精度Operator用法及压位
看到一道BZOJ1089，有高精度各种操作，可以当高精度模板了
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 | #include<cstdio>#define rad 10000//压位 using namespace std;int n,d,i,now;struct data{int v[2005],l;}f[30];//v为值，l为长度 void print(data a){printf("%d",a.v[a.l]);for(i=a.l-1;i;i--)printf("%04d",a.v[i]);}data operator*(data a,data b){//高精度乘法 data c; for(int i=1;i<=a.l+b.l;i++)c.v[i]=0; for(int i=1;i<=a.l;i++) for(int j=1;j<=b.l;j++) c.v[i+j-1]+=a.v[i]*b.v[j]; c.l=a.l+b.l-1; for(int i=1;i<=c.l;i++){ c.v[i+1]+=c.v[i]/rad; c.v[i]%=rad; } if(c.v[c.l+1])c.l++; return c;}data operator^(data a,int p){//高精度求快速幂 data c; c.l=1;c.v[1]=1; for(;p;p>>=1,a=a*a)if(p&1)c=c*a; return c;}data operator-(data a,data b){//高精度减高精度 for(int i=1;i<=a.l;i++){ a.v[i]-=b.v[i]; if(a.v[i]<0){ a.v[i]+=rad; a.v[i+1]--; } } while(a.l>1&&!a.v[a.l])a.l--; return a;}data operator+(data a,int p){//高精度加单精度 a.v[1]+=p; int now=1; while(a.v[now]>=rad){ a.v[now+1]+=a.v[now]/rad; a.v[now]%=rad; now++; } if(now>a.l)a.l=now; return a;}int main(){ scanf("%d%d",&n,&d); if(!d){puts("1");return 0;} f[0].l=f[0].v[1]=1; for(i=1;i<=d;i++)f[i]=(f[i-1]^n)+1; print(f[d]-f[d-1]);} |
2015年5月16日 14:09
然而你这个高精度连负数都不兹磁
2015年5月16日 14:11
而且还不兹磁高精度加高精度和高精度除法，而且高精度乘法还是暴力实现的。。。
2015年5月16日 18:24
胡主力您是很强嘛