如何优化

P1045 [NOIP2003 普及组] 麦森数

CJXJR @ 2019-07-18 18:26:09

#include<bits/stdc++.h>
using namespace std;
int c[10000005],l=1;
void work(int x)
{
    if(x==0) return ;
    c[1]=2;
    for(int i=1;i<x;i++)
    {
        for(int j=1;j<=l;j++)
        c[j]*=2;
        for(int k=1;k<=l;k++)
        {
            if(c[k]>=10)
            {
                c[k+1]+=c[k]/10;
                c[k]%=10;
            }
        }
        while(c[l+1]) l++;
        for(int k=1;k<=l;k++)
        {
            if(c[k]>=10)
            {
                c[k+1]+=c[k]/10;
                c[k]%=10;
            }
        }
    }
    while(c[l]==0 && l>1) l--;
}
int main()
{
    int p;
    cin>>p;
    work(p);
    cout<<l<<endl;
    int t=0;
    for(int i=500;i>=2;i--)
    {
        printf("%d",c[i]);
        t++;
        if(t==50) printf("\n"),t=0;
    }
    printf("%d",c[1]-1);
}

by ud2_ @ 2019-07-18 18:38:38

l = p * 0.3010299956639812 + 1;

by 枫初音斗颂皮 @ 2019-07-18 18:42:24

@龚宸进天才 快速幂不会吗?不会的话 右转快速幂模板题


|