http://poj.org/problem?id=1001
1 //10000000 100000
2
3 #include<iostream>
4 #include<cstring>
5 using namespace std;
6
7 char ts[7],tts[7]; //用于处理输入，删去点，签到0，点后0，统计小数点后数的个数保存至origu
8 short orig[7],ans[2500],t_ans[2500]; // t_ans存乘法运算的中间过程
9 int b,u,origu,BITS,origBITS;
10
11 void my_pow(int x){
12 if(x==1)
13 return ;
14
15
16 my_pow(x/2);
17 //接下来ans*ans
18 memset(t_ans,0,sizeof(t_ans));
19 for(int i=BITS-1;i>=0;i--){ //模拟乘法
20 int tpos=BITS-1-i; //tpos即 这次乘的结果存在t_ans哪一位 ,注意ans是高位在前， t_ans是高位在后的
21 for(int j=BITS-1;j>=0;j--){
22 t_ans[tpos]+=ans[j]*ans[i];
23 if(t_ans[tpos]>9){ //进位， 这个进位应该进一位就可以了，没有必要再向前推了，因为某一位算出来的最大也到不了100（？）
24 t_ans[tpos+1]+=t_ans[tpos]/10;
25 t_ans[tpos]%=10;
26 }
27 tpos++;
28 }
29 }
30
31 int _BITS=0; // 统计 2500-t_ans的位数
32 for(int i=2500-1;i>=0;i--)
33 if(t_ans[i]==0)
34 _BITS++;
35 else
36 break;
37 BITS=2500-_BITS; //更新BITS
38
39 for(int i=BITS-1,j=0;i>=0;i--) //高位在前存到ans里去
40 ans[j++]=t_ans[i];
41 u*=2; //因为是自己×自己 所以小数点后数字个数增加一倍
42
43
44 if(x%2){ //如果^b的b是奇数，还得再×一下初始的数即orig
45 //ans*orig
46
47 memset(t_ans,0,sizeof(t_ans));
48
49 for(int i=origBITS-1;i>=0;i--){
50 int tpos=origBITS-1-i;
51 for(int j=BITS-1;j>=0;j--){
52 t_ans[tpos]+=ans[j]*orig[i];
53 if(t_ans[tpos]>9){
54 t_ans[tpos+1]+=t_ans[tpos]/10;
55 t_ans[tpos]%=10;
56 }
57 tpos++;
58 }
59 }
60
61 _BITS=0;
62 for(int i=2500-1;i>=0;i--)
63 if(t_ans[i]==0)
64 _BITS++;
65 else
66 break;
67 BITS=2500-_BITS;
68
69 for(int i=BITS-1,j=0;i>=0;i--)
70 ans[j++]=t_ans[i];
71 u=u+origu;
72 }
73
74 }
75
76 void _print(){
77 if(u<BITS){
78 int _del=0;
79 for(int i=BITS-1;i>=BITS-u;i--){
80 if(ans[i]==0)
81 _del++;
82 else
83 break;
84 }
85 BITS-=_del;
86 u-=_del;
87 for(int i=0;i<BITS-u;i++)
88 printf("%d",ans[i]);
89 if(u){
90 printf(".");
91 for(int i=BITS-u;i<BITS;i++)
92 printf("%d",ans[i]);
93 }
94 printf("\n");
95 }
96 else{
97 int _del=0;
98 for(int i=BITS-1;i>=0;i--){
99 if(ans[i]==0)
100 _del++;
101 else
102 break;
103 }
104 BITS-=_del;
105 u-=_del;
106 printf(".");
107 for(int i=0;i<u-BITS;i++)
108 printf("0");
109 for(int i=0;i<BITS;i++)
110 printf("%d",ans[i]);
111 printf("\n");
112 }
113 }
114
115 int main(){
116 while(scanf("%s%d",ts,&b)==2){
117
118 u=0;
119
120 bool fg=1; //删点后0
121 for(int i=5;i>=0;i--)
122 if(fg&&ts[i]=='0')
123 ts[i]='\0';
124 else
125 fg=0;
126
127
128 for(int i=strlen(ts)-1;i>=0;i--) //删点，统计点后数个数u origu
129 if(ts[i]!='.')
130 u++;
131 else
132 break;
133 origu=u;
134
135
136 for(int i=strlen(ts)-u-1;ts[i]!='\0';i++) // 删点，点后前移
137 ts[i]=ts[i+1];
138
139
140 memset(tts,0,sizeof(tts)); //删前导0
141 int pos=0;
142 fg=1;
143 for(int i=0;ts[i]!='\0';i++)
144 if(ts[i]=='0'&&fg)
145 pos++;
146 else
147 fg=0;
148 for(int i=pos,j=0;ts[i]!='\0';i++)
149 tts[j++]=ts[i];
150
151 BITS=origBITS=0;
152 memset(ans,-1,sizeof(ans));
153 memset(orig,-1,sizeof(orig));
154 for(int i=0;tts[i]!='\0';i++){ //转到int数组
155 ans[i]=orig[i]=(int)(tts[i]-'0');
156 origBITS++;BITS++;
157 }
158
159 my_pow(b); //快速幂
160 _print(); //输出
161 }
162 return 0;
163 }
但是杭电这道同样的过不了，据说杭电会有00000.这样的测试数据？但是0能记成这种记法吗？还有人家明明说了R>0.0 , n>0为什么还非得找R=0N=0的特例啊，不管了，反正俺poj过了