UVA_485 Pascal's Triangle of Death_高精度

6701158	485	Pascal's Triangle of Death	Accepted	C++	0.130	2008-10-04 17:05:05
6701129	485	Pascal's Triangle of Death	Wrong answer	C++	0.130	2008-10-04 16:49:39
6701125	485	Pascal's Triangle of Death	Wrong answer	C++	0.130	2008-10-04 16:48:23
6701122	485	Pascal's Triangle of Death	Wrong answer	C++	0.120	2008-10-04 16:46:54

没看清楚题目的意思结果WA了几次，题目是说如果出现10^60以上的数的时候这行也要输出完以后才exit -_-!

1. #include <iostream>
2. using namespace std;
3. #define DIGIT   4
4. #define DEPTH   10000
5. #define MAX     100
6. typedef int bignum_t[MAX+1];
7. int read(bignum_t a,istream& is=cin){
8.     char buf[MAX*DIGIT+1],ch;
9.     int i,j;
10.     memset((void*)a,0,sizeof(bignum_t));
11.     if (!(is>>buf)) return 0;
12.     for (a[0]=strlen(buf),i=a[0]/2-1;i>=0;i--)
13.         ch=buf[i],buf[i]=buf[a[0]-1-i],buf[a[0]-1-i]=ch;
14.     for (a[0]=(a[0]+DIGIT-1)/DIGIT,j=strlen(buf);j<a[0]*DIGIT;buf[j++]='0');
15.     for (i=1;i<=a[0];i++)
16.         for (a[i]=0,j=0;j<DIGIT;j++)
17.             a[i]=a[i]*10+buf[i*DIGIT-1-j]-'0';
18.     for (;!a[a[0]]&&a[0]>1;a[0]--);
19.     return 1;
20. }
22. void write(const bignum_t a,ostream& os=cout){
23.     int i,j;
24.     for (os<<a[i=a[0]],i--;i;i--)
25.         for (j=DEPTH/10;j;j/=10)
26.             os<<a[i]/j%10;
27. }
29. int comp(const bignum_t a,const bignum_t b){
30.     int i;
31.     if (a[0]!=b[0])
32.         return a[0]-b[0];
33.     for (i=a[0];i;i--)
34.         if (a[i]!=b[i])
35.             return a[i]-b[i];
36.     return 0;
37. }
39. int comp(const bignum_t a,const int b){
40.     int c[12]={1};
41.     for (c[1]=b;c[c[0]]>=DEPTH;c[c[0]+1]=c[c[0]]/DEPTH,c[c[0]]%=DEPTH,c[0]++);
42.     return comp(a,c);
43. }
45. int comp(const bignum_t a,const int c,const int d,const bignum_t b){
46.     int i,t=0,O=-DEPTH*2;
47.     if (b[0]-a[0]<d&&c)
48.         return 1;
49.     for (i=b[0];i>d;i--){
50.         t=t*DEPTH+a[i-d]*c-b[i];
51.         if (t>0) return 1;
52.         if (t<O) return 0;
53.     }
54.     for (i=d;i;i--){
55.         t=t*DEPTH-b[i];
56.         if (t>0) return 1;
57.         if (t<O) return 0;
58.     }
59.     return t>0;
60. }
62. void add(bignum_t a,const bignum_t b){
63.     int i;
64.     for (i=1;i<=b[0];i++)
65.         if ((a[i]+=b[i])>=DEPTH)
66.             a[i]-=DEPTH,a[i+1]++;
67.     if (b[0]>=a[0])
68.         a[0]=b[0];
69.     else
70.         for (;a[i]>=DEPTH&&i<a[0];a[i]-=DEPTH,i++,a[i]++);
71.     a[0]+=(a[a[0]+1]>0);
72. }
74. void add(bignum_t a,const int b){
75.     int i=1;
76.     for (a[1]+=b;a[i]>=DEPTH&&i<a[0];a[i+1]+=a[i]/DEPTH,a[i]%=DEPTH,i++);
77.     for (;a[a[0]]>=DEPTH;a[a[0]+1]=a[a[0]]/DEPTH,a[a[0]]%=DEPTH,a[0]++);
78. }
80. void sub(bignum_t a,const bignum_t b){
81.     int i;
82.     for (i=1;i<=b[0];i++)
83.         if ((a[i]-=b[i])<0)
84.             a[i+1]--,a[i]+=DEPTH;
85.     for (;a[i]<0;a[i]+=DEPTH,i++,a[i]--);
86.     for (;!a[a[0]]&&a[0]>1;a[0]--);
87. }
89. void sub(bignum_t a,const int b){
90.     int i=1;
91.     for (a[1]-=b;a[i]<0;a[i+1]+=(a[i]-DEPTH+1)/DEPTH,a[i]-=(a[i]-DEPTH+1)/DEPTH*DEPTH,i++);
92.     for (;!a[a[0]]&&a[0]>1;a[0]--);
93. }
95. void sub(bignum_t a,const bignum_t b,const int c,const int d){
96.     int i,O=b[0]+d;
97.     for (i=1+d;i<=O;i++)
98.         if ((a[i]-=b[i-d]*c)<0)
99.             a[i+1]+=(a[i]-DEPTH+1)/DEPTH,a[i]-=(a[i]-DEPTH+1)/DEPTH*DEPTH;
100.     for (;a[i]<0;a[i+1]+=(a[i]-DEPTH+1)/DEPTH,a[i]-=(a[i]-DEPTH+1)/DEPTH*DEPTH,i++);
101.     for (;!a[a[0]]&&a[0]>1;a[0]--);
102. }
104. void mul(bignum_t c,const bignum_t a,const bignum_t b){
105.     int i,j;
106.     memset((void*)c,0,sizeof(bignum_t));
107.     for (c[0]=a[0]+b[0]-1,i=1;i<=a[0];i++)
108.         for (j=1;j<=b[0];j++)
109.             if ((c[i+j-1]+=a[i]*b[j])>=DEPTH)
110.                 c[i+j]+=c[i+j-1]/DEPTH,c[i+j-1]%=DEPTH;
111.     for (c[0]+=(c[c[0]+1]>0);!c[c[0]]&&c[0]>1;c[0]--);
112. }
114. void mul(bignum_t a,const int b){
115.     int i;
116.     for (a[1]*=b,i=2;i<=a[0];i++){
117.         a[i]*=b;
118.         if (a[i-1]>=DEPTH)
119.             a[i]+=a[i-1]/DEPTH,a[i-1]%=DEPTH;
120.     }
121.     for (;a[a[0]]>=DEPTH;a[a[0]+1]=a[a[0]]/DEPTH,a[a[0]]%=DEPTH,a[0]++);
122.     for (;!a[a[0]]&&a[0]>1;a[0]--);
123. }
125. void mul(bignum_t b,const bignum_t a,const int c,const int d){
126.     int i;
127.     memset((void*)b,0,sizeof(bignum_t));
128.     for (b[0]=a[0]+d,i=d+1;i<=b[0];i++)
129.         if ((b[i]+=a[i-d]*c)>=DEPTH)
130.             b[i+1]+=b[i]/DEPTH,b[i]%=DEPTH;
131.     for (;b[b[0]+1];b[0]++,b[b[0]+1]=b[b[0]]/DEPTH,b[b[0]]%=DEPTH);
132.     for (;!b[b[0]]&&b[0]>1;b[0]--);
133. }
135. void div(bignum_t c,bignum_t a,const bignum_t b){
136.     int h,l,m,i;
137.     memset((void*)c,0,sizeof(bignum_t));
138.     c[0]=(b[0]<a[0]+1)?(a[0]-b[0]+2):1;
139.     for (i=c[0];i;sub(a,b,c[i]=m,i-1),i--)
140.         for (h=DEPTH-1,l=0,m=(h+l+1)>>1;h>l;m=(h+l+1)>>1)
141.             if (comp(b,m,i-1,a)) h=m-1;
142.             else l=m;
143.     for (;!c[c[0]]&&c[0]>1;c[0]--);
144.     c[0]=c[0]>1?c[0]:1;
145. }
147. void div(bignum_t a,const int b,int& c){
148.     int i;
149.     for (c=0,i=a[0];i;c=c*DEPTH+a[i],a[i]=c/b,c%=b,i--);
150.     for (;!a[a[0]]&&a[0]>1;a[0]--);
151. }
153. void sqrt(bignum_t b,bignum_t a){
154.     int h,l,m,i;
155.     memset((void*)b,0,sizeof(bignum_t));
156.     for (i=b[0]=(a[0]+1)>>1;i;sub(a,b,m,i-1),b[i]+=m,i--)
157.         for (h=DEPTH-1,l=0,b[i]=m=(h+l+1)>>1;h>l;b[i]=m=(h+l+1)>>1)
158.             if (comp(b,m,i-1,a)) h=m-1;
159.             else l=m;
160.     for (;!b[b[0]]&&b[0]>1;b[0]--);
161.     for (i=1;i<=b[0];b[i++]>>=1);
162. }
164. int length(const bignum_t a){
165.     int t,ret;
166.     for (ret=(a[0]-1)*DIGIT,t=a[a[0]];t;t/=10,ret++);
167.     return ret>0?ret:1;
168. }
170. int digit(const bignum_t a,const int b){
171.     int i,ret;
172.     for (ret=a[(b-1)/DIGIT+1],i=(b-1)%DIGIT;i;ret/=10,i--);
173.     return ret%10;
174. }
176. int zeronum(const bignum_t a){
177.     int ret,t;
178.     for (ret=0;!a[ret+1];ret++);
179.     for (t=a[ret+1],ret*=DIGIT;!(t%10);t/=10,ret++);
180.     return ret;
181. }
183. void comp(int* a,const int l,const int h,const int d){
184.     int i,j,t;
185.     for (i=l;i<=h;i++)
186.         for (t=i,j=2;t>1;j++)
187.             while (!(t%j))
188.                 a[j]+=d,t/=j;
189. }
191. void convert(int* a,const int h,bignum_t b){
192.     int i,j,t=1;
193.     memset(b,0,sizeof(bignum_t));
194.     for (b[0]=b[1]=1,i=2;i<=h;i++)
195.         if (a[i])
196.             for (j=a[i];j;t*=i,j--)
197.                 if (t*i>DEPTH)
198.                     mul(b,t),t=1;
199.     mul(b,t);
200. }
202. void combination(bignum_t a,int m,int n){
203.     int* t=new int[m+1];
204.     memset((void*)t,0,sizeof(int)*(m+1));
205.     comp(t,n+1,m,1);
206.     comp(t,2,m-n,-1);
207.     convert(t,m,a);
208.     delete []t;
209. }
211. void permutation(bignum_t a,int m,int n){
212.     int i,t=1;
213.     memset(a,0,sizeof(bignum_t));
214.     a[0]=a[1]=1;
215.     for (i=m-n+1;i<=m;t*=i++)
216.         if (t*i>DEPTH)
217.             mul(a,t),t=1;
218.     mul(a,t);
219. }
221. #define SGN(x) ((x)>0?1:((x)<0?-1:0))
222. #define ABS(x) ((x)>0?(x):-(x))
224. int read(bignum_t a,int &sgn,istream& is=cin){
225.     char str[MAX*DIGIT+2],ch,*buf;
226.     int i,j;
227.     memset((void*)a,0,sizeof(bignum_t));
228.     if (!(is>>str)) return 0;
229.     buf=str,sgn=1;
230.     if (*buf=='-') sgn=-1,buf++;
231.     for (a[0]=strlen(buf),i=a[0]/2-1;i>=0;i--)
232.         ch=buf[i],buf[i]=buf[a[0]-1-i],buf[a[0]-1-i]=ch;
233.     for (a[0]=(a[0]+DIGIT-1)/DIGIT,j=strlen(buf);j<a[0]*DIGIT;buf[j++]='0');
234.     for (i=1;i<=a[0];i++)
235.         for (a[i]=0,j=0;j<DIGIT;j++)
236.             a[i]=a[i]*10+buf[i*DIGIT-1-j]-'0';
237.     for (;!a[a[0]]&&a[0]>1;a[0]--);
238.     if (a[0]==1&&!a[1]) sgn=0;
239.     return 1;
240. }
242. struct bignum{
243.     bignum_t num;
244.     int sgn;
245. public:
246. inline bignum(){memset(num,0,sizeof(bignum_t));num[0]=1;sgn=0;}
247. inline int operator!(){return num[0]==1&&!num[1];}
248. inline bignum& operator=(const bignum& a){memcpy(num,a.num,sizeof(bignum_t));sgn=a.sgn;return *this;}
249. inline bignum& operator=(const int a){memset(num,0,sizeof(bignum_t));num[0]=1;sgn=SGN(a);add(num,sgn*a);return *this;};
250. inline bignum& operator+=(const bignum& a){if(sgn==a.sgn)add(num,a.num);else if(sgn&&a.sgn){int ret=comp(num,a.num);if(ret>0)sub(num,a.num);else if(ret<0){bignum_t t;
251.     memcpy(t,num,sizeof(bignum_t));memcpy(num,a.num,sizeof(bignum_t));sub(num,t);sgn=a.sgn;}else memset(num,0,sizeof(bignum_t)),num[0]=1,sgn=0;}else if(!sgn)memcpy(num,a.num,sizeof(bignum_t)),sgn=a.sgn;return *this;}
252. inline bignum& operator+=(const int a){if(sgn*a>0)add(num,ABS(a));else if(sgn&&a){int ret=comp(num,ABS(a));if(ret>0)sub(num,ABS(a));else if(ret<0){bignum_t t;
253.     memcpy(t,num,sizeof(bignum_t));memset(num,0,sizeof(bignum_t));num[0]=1;add(num,ABS(a));sgn=-sgn;sub(num,t);}else memset(num,0,sizeof(bignum_t)),num[0]=1,sgn=0;}else if(!sgn)sgn=SGN(a),add(num,ABS(a));return *this;}
254. inline bignum operator+(const bignum& a){bignum ret;memcpy(ret.num,num,sizeof(bignum_t));ret.sgn=sgn;ret+=a;return ret;}
255. inline bignum operator+(const int a){bignum ret;memcpy(ret.num,num,sizeof(bignum_t));ret.sgn=sgn;ret+=a;return ret;}
256. inline bignum& operator-=(const bignum& a){if(sgn*a.sgn<0)add(num,a.num);else if(sgn&&a.sgn){int ret=comp(num,a.num);if(ret>0)sub(num,a.num);else if(ret<0){bignum_t t;
257.     memcpy(t,num,sizeof(bignum_t));memcpy(num,a.num,sizeof(bignum_t));sub(num,t);sgn=-sgn;}else memset(num,0,sizeof(bignum_t)),num[0]=1,sgn=0;}else if(!sgn)add(num,a.num),sgn=-a.sgn;return *this;}
258. inline bignum& operator-=(const int a){if(sgn*a<0)add(num,ABS(a));else if(sgn&&a){int ret=comp(num,ABS(a));if(ret>0)sub(num,ABS(a));else if(ret<0){bignum_t t;
259.     memcpy(t,num,sizeof(bignum_t));memset(num,0,sizeof(bignum_t));num[0]=1;add(num,ABS(a));sub(num,t);sgn=-sgn;}else memset(num,0,sizeof(bignum_t)),num[0]=1,sgn=0;}else if(!sgn)sgn=-SGN(a),add(num,ABS(a));return *this;}
260. inline bignum operator-(const bignum& a){bignum ret;memcpy(ret.num,num,sizeof(bignum_t));ret.sgn=sgn;ret-=a;return ret;}
261. inline bignum operator-(const int a){bignum ret;memcpy(ret.num,num,sizeof(bignum_t));ret.sgn=sgn;ret-=a;return ret;}
262. inline bignum& operator*=(const bignum& a){bignum_t t;mul(t,num,a.num);memcpy(num,t,sizeof(bignum_t));sgn*=a.sgn;return *this;}
263. inline bignum& operator*=(const int a){mul(num,ABS(a));sgn*=SGN(a);return *this;}
264. inline bignum operator*(const bignum& a){bignum ret;mul(ret.num,num,a.num);ret.sgn=sgn*a.sgn;return ret;}
265. inline bignum operator*(const int a){bignum ret;memcpy(ret.num,num,sizeof(bignum_t));mul(ret.num,ABS(a));ret.sgn=sgn*SGN(a);return ret;}
266. inline bignum& operator/=(const bignum& a){bignum_t t;div(t,num,a.num);memcpy(num,t,sizeof(bignum_t));sgn=(num[0]==1&&!num[1])?0:sgn*a.sgn;return *this;}
267. inline bignum& operator/=(const int a){int t;div(num,ABS(a),t);sgn=(num[0]==1&&!num[1])?0:sgn*SGN(a);return *this;}
268. inline bignum operator/(const bignum& a){bignum ret;bignum_t t;memcpy(t,num,sizeof(bignum_t));div(ret.num,t,a.num);ret.sgn=(ret.num[0]==1&&!ret.num[1])?0:sgn*a.sgn;return ret;}
269. inline bignum operator/(const int a){bignum ret;int t;memcpy(ret.num,num,sizeof(bignum_t));div(ret.num,ABS(a),t);ret.sgn=(ret.num[0]==1&&!ret.num[1])?0:sgn*SGN(a);return ret;}
270. inline bignum& operator%=(const bignum& a){bignum_t t;div(t,num,a.num);if (num[0]==1&&!num[1])sgn=0;return *this;}
271. inline int operator%=(const int a){int t;div(num,ABS(a),t);memset(num,0,sizeof(bignum_t));num[0]=1;add(num,t);return t;}
272. inline bignum operator%(const bignum& a){bignum ret;bignum_t t;memcpy(ret.num,num,sizeof(bignum_t));div(t,ret.num,a.num);ret.sgn=(ret.num[0]==1&&!ret.num[1])?0:sgn;return ret;}
273. inline int operator%(const int a){bignum ret;int t;memcpy(ret.num,num,sizeof(bignum_t));div(ret.num,ABS(a),t);memset(ret.num,0,sizeof(bignum_t));ret.num[0]=1;add(ret.num,t);return t;}
274. inline bignum& operator++(){*this+=1;return *this;}
275. inline bignum& operator--(){*this-=1;return *this;};
276. inline int operator>(const bignum& a){return sgn>0?(a.sgn>0?comp(num,a.num)>0:1):(sgn<0?(a.sgn<0?comp(num,a.num)<0:0):a.sgn<0);}
277. inline int operator>(const int a){return sgn>0?(a>0?comp(num,a)>0:1):(sgn<0?(a<0?comp(num,-a)<0:0):a<0);}
278. inline int operator>=(const bignum& a){return sgn>0?(a.sgn>0?comp(num,a.num)>=0:1):(sgn<0?(a.sgn<0?comp(num,a.num)<=0:0):a.sgn<=0);}
279. inline int operator>=(const int a){return sgn>0?(a>0?comp(num,a)>=0:1):(sgn<0?(a<0?comp(num,-a)<=0:0):a<=0);}
280. inline int operator<(const bignum& a){return sgn<0?(a.sgn<0?comp(num,a.num)>0:1):(sgn>0?(a.sgn>0?comp(num,a.num)<0:0):a.sgn>0);}
281. inline int operator<(const int a){return sgn<0?(a<0?comp(num,-a)>0:1):(sgn>0?(a>0?comp(num,a)<0:0):a>0);}
282. inline int operator<=(const bignum& a){return sgn<0?(a.sgn<0?comp(num,a.num)>=0:1):(sgn>0?(a.sgn>0?comp(num,a.num)<=0:0):a.sgn>=0);}
283. inline int operator<=(const int a){return sgn<0?(a<0?comp(num,-a)>=0:1):(sgn>0?(a>0?comp(num,a)<=0:0):a>=0);}
284. inline int operator==(const bignum& a){return (sgn==a.sgn)?!comp(num,a.num):0;}
285. inline int operator==(const int a){return (sgn*a>=0)?!comp(num,ABS(a)):0;}
286. inline int operator!=(const bignum& a){return (sgn==a.sgn)?comp(num,a.num):1;}
287. inline int operator!=(const int a){return (sgn*a>=0)?comp(num,ABS(a)):1;}
288. inline int operator[](const int a){return digit(num,a);}
289. friend inline istream& operator>>(istream& is,bignum& a){read(a.num,a.sgn,is);return is;}
290. friend inline ostream& operator<<(ostream& os,const bignum& a){if(a.sgn<0)os<<'-';write(a.num,os);return os;}
291. friend inline bignum sqrt(const bignum& a){bignum ret;bignum_t t;memcpy(t,a.num,sizeof(bignum_t));sqrt(ret.num,t);ret.sgn=ret.num[0]!=1||ret.num[1];return ret;}
292. friend inline bignum sqrt(const bignum& a,bignum& b){bignum ret;memcpy(b.num,a.num,sizeof(bignum_t));sqrt(ret.num,b.num);ret.sgn=ret.num[0]!=1||ret.num[1];b.sgn=b.num[0]!=1||ret.num[1];return ret;}
293. inline int length(){return ::length(num);}
294. inline int zeronum(){return ::zeronum(num);}
295. inline bignum C(const int m,const int n){combination(num,m,n);sgn=1;return *this;}
296. inline bignum P(const int m,const int n){permutation(num,m,n);sgn=1;return *this;}
297. };
298. bignum tmp[10000],k;
299. int main()
300. {
301.     int i,len=3;
302.     bool flag=false;
303.     cout<<1<<endl;
304.     cout<<1<<' '<<1<<endl;
305.     tmp[0]=1;
306.     tmp[1]=2;
307.     tmp[2]=1;
308.     while(1)
309.         {
310.             for(i=0;i<len-1;i++){if(tmp[i].length()>=61)flag=true;cout<<tmp[i]<<' ';}cout<<1<<endl;
311.             tmp[0]=1;
312.             len++;
313.             for(i=len-2;i>=1;i--)tmp[i]+=tmp[i-1];
314.             tmp[len-1]=1;
315.             if(flag)break;
316.         }
317.     return 0;
318. }