【学习笔记】[ABC236Ex] Distinct Multiples

典中典题目。

考虑怎么转化成普通的容斥形式🤔

发现可以将$S$拆分成若干个集合相加的形式，然后分开求$\text{lcm}$算贡献。

感觉和斯特林数有关系。

破防了，题解告诉我们这个容斥系数是$(-1{)}^{|S|-1}(|S|-1)!$😅

可以通过对边集容斥得到这个优美的结果。

可以参考洛谷题解。

发现可以将$a_i=a_j$的两个点之间连边，这样我们转化成了图计数的问题。

设$C(E)$表示连上边集$E$后的连通块结果，$F(P)$表示连通块为$P$时的答案，也就是连通块方案数的乘积。

那么根据经典的容斥原理，答案是${\sum }_E(-1{)}^{|E|}F(C(E))$。

改变求和顺序得到：${\sum }_{P}F(P){\sum }_{C(E)=P}(-1{)}^{|E|}$

首先要做过这道题 [集训队作业2013] 城市规划 。

记$G(x)$表示简单无向图的$EGF$，$F(x)$表示简单连通图的$EGF$，发现$G(x)={\sum }_{n=0}^{\infty }{2}^{\left(\genfrac{}{}{0ex}{}{n}{2}\right)}\frac{x^n}{n!}({\sum }_{i=0}^{\left(\genfrac{}{}{0ex}{}{n}{2}\right)}(-1{)}^i\left(\genfrac{}{}{0ex}{}{\left(\genfrac{}{}{0ex}{}{n}{2}\right)}{i}\right))=1+x$，$G(x)=1+\frac{F(x)}{1!}+\frac{F(x{)}^2}{2!}+...+\frac{F(x{)}^n}{n!}=\text{exp(F(x))}$，手算便知$F(x)$，这样我们得到了容斥系数。最后$DP$即可。

暴力实现，复杂度$O(3^n)$。

为什么是圆排列？我感觉有不用多项式的方法，但是我想不到。

$\text{remark}$ 感觉对多项式的掌握还差的很远。。。

#include<bits/stdc++.h>
#define ll long long
#define pb push_back
#define fi first
#define se second
#define inf 0x3f3f3f3f
using namespace std;
const int mod=998244353;
int n;
ll m,fac[16],a[16],f[1<<16],g[1<<16];
ll gcd(ll x,ll y){
    return y==0?x:gcd(y,x%y);
}
ll lcm(ll x,ll y){
    ll z=gcd(x,y);
    if((__int128)x/z*y>m)return m+1;
    return x/z*y;
}
int main(){
    ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);
    cin>>n>>m;for(int i=0;i<n;i++)cin>>a[i];
    fac[0]=1;for(int i=1;i<n;i++)fac[i]=fac[i-1]*i%mod;
    for(int i=1;i<1<<n;i++){
        ll tmp=0;int sz=0;
        for(int j=0;j<n;j++){
            if(i>>j&1){
                sz++;
                if(tmp==0)tmp=a[j];
                else tmp=lcm(tmp,a[j]);
            }
        }
        g[i]=m/tmp%mod*fac[sz-1]%mod*((sz%2==0)?-1:1);
    }
    f[0]=1;
    for(int i=1;i<1<<n;i++){
        int k=0;while(!(i>>k&1))k++;
        for(int j=i;j;j=(j-1)&i){
            if(j>>k&1)f[i]=(f[i]+f[i-j]*g[j])%mod;
        }
    }
    cout<<(f[(1<<n)-1]+mod)%mod;
}