• P1045 [NOIP2003 普及組] 麥森數

解析

看似只是正常的一個高精
然而

暗藏殺機

一開始隨手那么一寫
。。。
(即使用了快速冪)時間復雜度過于感人
后來我們發現：
第一問位數的計算不必真的算出來，只需把2的p次冪轉化為10的k次冪即可（具體請參看高中數學必修1）

(為了大家，我還是敲一下吧。。）
2^p=10的log10(2^p)次方=10^p*log10(2)
k應該等于p*log10（2）

我們就可以使用log10() 函數水掉第一問
至于第二問，既然只問后500位，就可以暴力刨掉后面了~~~
（慷慨的我們又送了十位，不客氣）
PS：別忘了50位一換行（別問我為啥要強調這個。。。）

代碼

#include<cstdio> #include<cmath> #include<cstring> #include<algorithm> #include<string> #include<climits> using namespace std; int p,ans[3100002],res[3100002]; double r=log10(2); int *a=&ans[1],*b=&res[1]; void cheng(int *x,int *y){int l1=*(x-1);int l2=*(y-1);int f[31002]={ };for(int i=1;i<=l1;i++){for(int j=1;j<=l2;j++){f[i+j-1] += *(x+i-1) * *(y+j-1);}}f[0]=l1+l2+5;for(int i=1;i<=f[0];i++){f[i+1]+=f[i]/10;f[i]%=10;}while(f[f[0]]==0) f[0]--;f[0]=min(f[0],510);for(int i=0;i<=f[0];i++){*(x+i-1)=f[i];} } void ksm(int n){while(n){if(n%2==1) cheng(a,b);cheng(b,b);n /= 2;} } int main(){ans[0]=res[0]=ans[1]=1;res[1]=2;scanf("%d",&p);ksm(p);int as=floor(p*r)+1;printf("%d",as);ans[1]--;int place=1;while(ans[place]<0){ans[place]+=10;ans[place+1]--;place++;}for(int i=500;i>=1;i--){if((500-i)%50==0) printf("\n");printf("%d",ans[i]);}return 0; }

AC快樂！

總結

以上是生活随笔為你收集整理的高精度：麦森数*（洛谷P1045）的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。