CF1794D 题解

一、题目描述：

　　一个正整数 $m$ 可以被唯一分解成 $p_1^{e_1} \times p_2^{e_2} \times ...\times p_k^{e_k}$ 的形式，其中 $p_1,p_2,...,p_k$为互不相同的质数，$e_1,e_2,...,e_k$ 为正整数。

　　定义一个可重集 $f(m)$ 为 $\{p_1,e_1,p_2,e_2,...,p_k,e_k\}$ 。现在给定正整数 $n$ 和大小为 $2\times n$ 集合 $S$ ，求有多少个 $m$ 能使得 $f(m)=S$。

　　由于答案可能很大，请对 $998$ $244$ $353$ 取模。数据范围：$1\leq n\leq 2022$。

---

 二、解题思路：

　　思路不是我的，我是看一个同学的代码看懂的。

　　其实比较容易想到要统计素数个数，合数个数，然后套组合数多重集的板子。

　　但是方案数很多，甚至可能有 $2^n$ 种，所以我考场没做出来。

 

　　说到方案数，其实是一个计数 $dp$，但我对 $dp$ 本来就一窍不通，所以没做出来也很正常。

　　用一个类似背包的数组 $f_{i,j}$ 来表示状态，笼统一点理解：$f_{i,j}$ 表示选了前 $i$ 个数中的 $j$ 个数作底数的方案数。

　　但实际上并不完全是这样，因为有一部分质数也有可能做质数，所以还乘了逆元。看代码就知道了。时间复杂度 $O(n^2)$。

---

三、完整代码：

#include<iostream>
#define N 4050
#define M 1000010
#define lim 1000000
#define ll long long
#define mod 998244353
using namespace std;
ll n,cp,cc,cnt,ans;
ll a[N],s[M],v1[M],v2[M],pri[M];
ll jc[M],inv[M],f[N][N],p[N],c[N];
void Prime()
{
    for(ll i=2;i<=lim;i++)
    {
        if(!v1[i])    pri[++cnt]=i;
        for(ll j=1;j<=cnt;j++)
        {
            if(i*pri[j]>lim)
                break;
            v1[i*pri[j]]=1;
            if(i%pri[j]==0)
                break;
        }
    }
}
ll qsm(ll base,ll q)
{
    ll res=1;
    while(q)
    {
        if(q&1)    res*=base,res%=mod;
        base*=base,base%=mod,q>>=1;
    }
    return res;
}
void pre_work()
{
    Prime();
    v1[1]=jc[0]=1;
    for(ll i=1;i<=lim;i++)
        jc[i]=jc[i-1]*i%mod;
    inv[lim]=qsm(jc[lim],mod-2);
    for(ll i=lim;i>=1;i--)
        inv[i-1]=inv[i]*i%mod;
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);cout.tie(0);
    cin>>n;pre_work();
    for(ll i=1;i<=n*2;i++)
    {
        cin>>a[i],s[a[i]]++;
        if(!v2[a[i]])
            if(!v1[a[i]])    p[++cp]=a[i],v2[a[i]]=1;
            else            c[++cc]=a[i],v2[a[i]]=1;
    }
    for(ll i=0;i<=cp;i++)
        f[i][0]=1;
    for(ll i=1;i<=cp;i++)
        for(ll j=0;j<=n;j++)
        {
            f[i][j]=f[i-1][j-1]*inv[s[p[i]]-1]%mod;
            (f[i][j]+=f[i-1][j]*inv[s[p[i]]])%=mod;
        }
    ans=f[cp][n]*jc[n]%mod;
    for(ll i=1;i<=cc;i++)
        (ans*=inv[s[c[i]]])%=mod;
    cout<<ans<<'\n';
    return 0;
}

---

四、写题心得：

　　感觉有关阶乘，逆元的题都特别有意思，虽然我大多数都不会。不过加油吧！拜拜！