高精度模板

#include <iostream>#include <algorithm>#include <time.h>using namespace std;#define maxsize 1005struct hp{ int len; int s[maxsize+1];};void copy(hp &a,hp &b){a.len=b.len;for(int i=1;i<=b.len;i++)a.s[i]=b.s[i];}void input(hp &a,string str){ int i; while(str[0]=='0' && str.size()!=1) str.erase(0,1); a.len=(int)str.size(); for(i=1;i<=a.len;++i) a.s[i]=str[a.len-i]-48; for (i=a.len+1;i<=maxsize;++i) a.s[i]=0;}void print(const hp &y){ int i; for(i=y.len;i>=1;i--)printf("%d",y.s[i]); cout<<endl;}void plusk(const hp &a,const hp &b,hp &c) //高精度加法{ int i,len; for(i=1;i<=maxsize;i++) c.s[i]=0; if(a.len>b.len) len=a.len; else len=b.len; for(i=1;i<=len;i++) { c.s[i]+=a.s[i]+b.s[i]; if(c.s[i]>=10) { c.s[i]-=10; c.s[i+1]++; } } if(c.s[len+1]>0) len++; c.len=len;}void subtract(const hp &a,const hp &b,hp &c) //高精度减法{ int i,len; for(i=1;i<=maxsize;i++) c.s[i]=0; if(a.len>b.len) len=a.len; else len=b.len; for(i=1;i<=len;i++) { c.s[i]+=a.s[i]-b.s[i]; if(c.s[i]<0) { c.s[i]+=10; c.s[i+1]--; } } while(len>1&&c.s[len]==0) len--; c.len=len;}int compare(const hp &a,const hp &b){ int len; if(a.len>b.len) len=a.len; else len=b.len; while(len>0 && a.s[len]==b.s[len]) len--; if(len==0) return 0; else return a.s[len]-b.s[len];}void multiply(const hp &a,int b,hp &c) //高精度*单精度{ int i,len; for(i=1;i<=maxsize;i++) c.s[i]=0; len=a.len; for(i=1;i<=len;i++) { c.s[i]+=a.s[i]*b; c.s[i+1]+=c.s[i]/10; c.s[i]%=10; } len++; while(c.s[len]>=10) { c.s[len+1]+=c.s[len]/10; c.s[len]%=10; len++; } while(len>1&&c.s[len]==0) len--; c.len=len;}void multiplyh(const hp &a,const hp &b,hp &c) //高精度*高精度{ int i,j,len; for(i=1;i<=maxsize;i++) c.s[i]=0; for(i=1;i<=a.len;i++) for(j=1;j<=b.len;j++) { c.s[i+j-1]+=a.s[i]*b.s[j]; c.s[i+j]+=c.s[i+j-1]/10; c.s[i+j-1]%=10; } len=a.len+b.len+1; while(len>1&&c.s[len]==0) len--; c.len=len;}void divide(const hp &a,int b,hp &c,int &d) //高精度/单精度 　{d为余数}{ int i,len; for(i=1;i<=maxsize;i++) c.s[i]=0; len=a.len; d=0; for(i=len;i>=1;i--) { d=d*10+a.s[i]; c.s[i]=d/b; d%=b; } while(len>1&&c.s[len]==0) len--; c.len=len;}void multiply10(hp &a) //高精度*10{ int i; for(i=a.len;i>=1;i--) a.s[i+1]=a.s[i]; a.s[1]=0; a.len++; while(a.len>1&&a.s[a.len]==0) a.len--;}void divideh(const hp &a,const hp &b,hp &c,hp &d) //高精度/高精度 　{d为余数}{ hp e; int i,len; for(i=1;i<=maxsize;i++) { c.s[i]=0; d.s[i]=0; } len=a.len; d.len=1; for(i=len;i>=1;i--) { multiply10(d); d.s[1]=a.s[i]; while(compare(d,b)>=0) { subtract(d,b,e); d=e; c.s[i]++; } } while(len>1&&c.s[len]==0) len--; c.len=len;}hp sum[1003];int main(){int n;//{n=1002;hp a, b,c,tmp;c.len=0;a.s[1]=1;for(int i=2;i<1002;i++){b.s[i]=0;a.s[i]=0;}for(int i=0;i<1002;i++){sum[i].len=1;memset(sum[i].s,0,sizeof(sum[i].s));}a.len=1;b.s[1]=2 ;sum[1].s[1]=2;sum[2].s[1]=1;b.len=1;int pos=2;if(n>2){for(int i=1;i<n;i++){plusk(a,b,c);//copy(tmp,sum[pos]);copy(sum[pos],c);int pos2;pos2=pos+1;//plusk(sum[pos],a,sum[pos2]);copy(a,b);copy(b,c);pos++;}}//print(sum[pos]);}while(cin>>n){print(sum[n]);}}