2021牛客寒假算法基础集训营2 I.牛牛的“质因数”

题目链接

思路

比赛的过程中想到了先筛出所有质数，然后通过相互之间的关系进行DP状态的转移，观察题目所给的样例不难得出，一个数的结果就是其任意两因子的结合。例如6可以被分解为2和3，其结果便是2*10+3=23。但此处我选择了选出每个数的最小质因子，然后将这个最小质因子的结果扩大至可以将另一个因子的数值合并。可这样的写法存在一个弊端，即过程中不可取模，因为一旦取模就会影响到最小质因子所对应的因子的位数，结果自然就会出现问题。正确的写法应当与此相反，寻找其最大的质因子，然后将最大质因子所对应的因子扩大至可以将最大质因子的数值合并。这样的写法优越性在于过程中取模并不会影响其正确性，因为被加的最大质因子本身不会因为取模而使位数出现误差，从而解决了超出范围的问题。当然，这样改变写法后似乎不可以再提前打表，否则寻找最大因子的复杂度会远超寻找最小因子的复杂度，导致程序运行超时，所以建议采取埃筛的写法，详见代码。(文中所提及的问题也可以通过优化解决，大家可以尝试对其进行改进)

代码

#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
const int maxn=1e7+5;
typedef long long ll;
const int inf=0x3f3f3f3f;
ll k,n,ans1,dp[maxn];
int cal(int x)
{
    int ans = 10;
    for (int i = 1; i <= x; i++)
        ans *=10;
    return ans;
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);cout.tie(0);
    cin>>n;
    for(int i=2;i<=n;i++)
    {
        if(!dp[i])
        {
            dp[i]=i;
            for(int j=2*i;j<=n;j+=i)
                dp[j]=(dp[j/i]*cal(log10(i))%mod+i)%mod;
        }
    }
    for(int i=2;i<=n;i++)
        ans1=(ans1+dp[i])%mod;
    cout<<ans1<<endl;
}