一樣搬之前初三寫的。。。
其實還有個val.2,搬到文章裏去了函數
@spa
在作一道斯特林數的時候被卡高精。。。因而滾來寫一些簡單的高精code
struct bigint{ int length,num[maxn]; bigint(){ memset(num,0,sizeof(num)); length=0; } };
注意:length存0開始的位數,實際有效位數爲0~length-1
構造函數能寫就寫it
inline void change(char a[],bigint &b){ b.length=strlen(a); for(int i=0;i<b.length;i++){ b.num[i]=a[b.length-i-1]-'0';//注意要減去‘0’ } }
inline void add(bigint &a,bigint &b){ int ll=max(a.length,b.length); int exd=0; for(int i=0;i<ll;i++){ int tmp=a.num[i]+b.num[i]+exd; a.num[i]=tmp%10; exd=tmp/10; } a.length=ll; if(exd){ a.num[ll]=exd; a.length++; } }
變量維護進位,模擬便可class
inline void mul(bigint &a,int b){ for(int i=0;i<a.length;i++){ a.num[i]*=b; } for(int i=0;i<a.length;i++){ a.num[i+1]+=a.num[i]/10; a.num[i]%=10; }//此時ai可能不止一位 int m=a.num[a.length]; while(m){ a.num[(++a.length)-1]=m%10; m/=10; }//最高位長度處理 }
inline void div(bigint &a,int b){ int d=0; for(int i=a.length-1;i>=0;i--){ int tmp=d*10+a.num[i]; a.num[i]=tmp/b; d=tmp%b; } while(a.num[a.length-1]==0){ a.length--; } //去除高位0 }
高精和高精乘除與壓位高精就以後再寫吧 (變量
(咕構造函數