[USACO20OPEN] Exercise G

题面:Luogu
题解:
大概是和这个一样的
只是那个是求方案数,这个是求和
所以我们在转移时

\[f_{i,j}=\sum f_{i-1,j-p_i^{k}}*p_i^k~(k \ge 0,p_i^k\le j)
\]

#include<bits/stdc++.h>
using namespace std;
template<typename T>
inline void read(T& x)
{
    x=0;char c=getchar();
    while(!isdigit(c)) c=getchar();
    while(isdigit(c)) x=x*10+c-'0',c=getchar();
}
#define ll long long
#define maxn 10005
int prime[maxn],vis[maxn],pnum;
void eular(int maxnum=10000)
{
    for(int i=2;i<=maxnum;++i)
    {
        if(!vis[i]) prime[++pnum]=i;
        for(int j=1;j<=pnum&&i*prime[j]<=maxnum;++j)
        {
            vis[i*prime[j]]=1;
            if(!(i%prime[j])) break;
        }
    }
}
ll dp[maxn],m;
int main()
{
    int n;read(n);read(m);
    eular(n);
    dp[0]=1;
    for(int i=1;i<=pnum;++i)
        for(int j=n;j>=prime[i];--j)
            for(int k=prime[i];k<=j;k*=prime[i])
                dp[j]=(dp[j]+dp[j-k]*k)%m;
    ll sum=1;
    for(int i=1;i<=n;++i) sum+=dp[i],sum%=m;
    printf("%lld\n",sum);
    return 0;
}

原文链接: https://www.cnblogs.com/123789456ye/p/12739188.html

欢迎关注

微信关注下方公众号,第一时间获取干货硬货;公众号内回复【pdf】免费获取数百本计算机经典书籍;

也有高质量的技术群,里面有嵌入式、搜广推等BAT大佬

    [USACO20OPEN] Exercise G

原创文章受到原创版权保护。转载请注明出处:https://www.ccppcoding.com/archives/343440

非原创文章文中已经注明原地址,如有侵权,联系删除

关注公众号【高性能架构探索】,第一时间获取最新文章

转载文章受原作者版权保护。转载请注明原作者出处!

(0)
上一篇 2023年3月2日 上午2:07
下一篇 2023年3月2日 上午2:07

相关推荐