30分求调,玄关

P1045 [NOIP2003 普及组] 麦森数

Foraver @ 2024-02-25 19:24:45

#include<bits/stdc++.h>
#define int long long
using namespace std;
int f[100000],tot=500,a;
struct big
{
    int num[10005]={0};
    int size;
};
big input(string a)
{
    big b;
    b.size=a.size();
    for(int i=a.size()-1;i>=0;i--)
    {
        b.num[b.size-i-1]=a[i]-'0';
    }
    return b;
}
void clear(big &a)
{
    for(int i=a.size-1;i>=1;i--)
    {
        if(a.num[i]==0) a.size--;
        else return ;
    }
}
big mul(big a,big b)
{
    big c;
    c.size=a.size+b.size+2;
    for(int i=0;i<=501;i++)
    {
        for(int j=0;j<=501;j++)
        {
            c.num[i+j]+=a.num[i]*b.num[j];
        }
    }
    for(int i=0;i<=501;i++)
    {
        c.num[i+1]+=c.num[i]/10;
        c.num[i]%=10;
    }
    return c;
}
big sub(big a,big b)
{
    big c;
    c.size=a.size;
    int flag=0;
    for(int i=0;i<=c.size;i++)
    {
        c.num[i]=a.num[i]-b.num[i]-flag;
        flag=(c.num[i]<0);
        c.num[i]=(c.num[i]+10)%10;
    }
    return c;
}
big YI=input("1");
void print(big b)
{
    int cnt=0;
    b=sub(b,YI);
    clear(b);
    cout<<floor(a*log10(2))+1<<"\n";
    for(int i=0;i<500;i++)f[tot--]=b.num[i];
    for(int i=1;i<=500;i++)
    {
        printf("%lld",f[i]);
        cnt++;
        if(cnt==50)
        {
            printf("\n");
            cnt=0;
        }
    }
}
bool cmp(big a,big b)
{
    clear(a);
    clear(b);
    if(a.size>b.size) return 1;
    if(a.size<b.size) return 0;
    for(int i=a.size-1;i>=0;i--)
    {
        if(a.num[i]>b.num[i]) return 1;
        if(a.num[i]<b.num[i]) return 0;
    }
    return 1;
}
signed main()
{
    cin>>a;
    big x=input("2"),sum=input("1");
    for(int i=1;i<=a;i++)sum=mul(x,sum);
    print(sum);
    return 0;
}

by Foraver @ 2024-02-25 19:25:07

TLE+RE


by Chenyufeng040525 @ 2024-03-06 22:06:33

@Foraver

TLE需要你使用快速幂或者增大每次乘的数(至少是65536),两种方法题解都有可以看看怎么写。

RE是因为你sub函数中c.side因为mul函数加得太大了,导致下标访问超过数组长度,或许你可以只取500位就行了。


by Foraver @ 2024-03-08 19:45:14

@Chenyufeng040525 谢谢


|