【洛谷P4916】魔力环（Burnside）（组合数学）

传送门

---

题解：

看到旋转同构，考虑Burnside。

旋转同构中常见的phi的转化就不提了。

直接考虑循环个数为$d$的时候的答案。

显然只有当$d\mid n$且$d\mid m$的时候答案非0。

特判掉$n=m$和$m=0$的情况。

现在环上必然有黑色和白色。

设$N=n/d,M=m/d$，接下来分类讨论。

如果$M\le K$，随便放就行了，反正有白色，黑色连不起来。

否则，我们断环为链，把所有白色拿出来，枚举断掉的地方首尾链接起来有多少黑色的，然后把剩下的黑色分成若干组，每组可以为空但是不允许超过$K$。直接枚举有多少组超过$K$暴力容斥即可。

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

using std::cerr;
using std::cout;

cs int mod=998244353;
inline int add(int a,int b){a+=b-mod;return a+(a>>31&mod);}
inline int dec(int a,int b){a-=b;return a+(a>>31&mod);}
inline int mul(int a,int b){ll r=(ll)a*b;return r>=mod?r%mod:r;}
inline int power(int a,int b,int res=1){
	for(;b;b>>=1,a=mul(a,a))(b&1)&&(res=mul(res,a));
	return res;
}
inline void Inc(int &a,int b){a+=b-mod;a+=a>>31&mod;}
inline void Dec(int &a,int b){a-=b;a+=a>>31&mod;}
inline void Mul(int &a,int b){a=mul(a,b);}

int n,m,K;

cs int N=1e5+7;

int p[N],pc;
bool mark[N];
int phi[N];

inline void sieves(int n){
	phi[1]=1;
	for(int re i=2;i<=n;++i){
		if(!mark[i])p[++pc]=i,phi[i]=i-1;
		for(int re j=1;i*p[j]<=n;++j){
			mark[i*p[j]]=true;
			if(i%p[j]){phi[i*p[j]]=phi[i]*(p[j]-1);}
			else {
				phi[i*p[j]]=phi[i]*p[j];
				break;
			}
		}
	}
}

int fac[N<<1],ifac[N<<1];
inline void init_inv(int n){
	fac[0]=fac[1]=1;
	for(int re i=2;i<=n;++i)fac[i]=mul(fac[i-1],i);
	ifac[n]=power(fac[n],mod-2);
	for(int re i=n-1;~i;--i)ifac[i]=mul(ifac[i+1],i+1);
}

inline int C(int n,int m){return n>=0&&m>=0&&n>=m?mul(fac[n],mul(ifac[m],ifac[n-m])):0;}
inline int calc(int m,int n){return C(n+m-1,n-1);}

inline int Put(int m,int n,int k){
	if(m<0)return 0;
	int res=0;
	for(int re i=0;i<=n;++i)
	if(m>=i*(k+1)){
		int val=mul(calc(m-i*(k+1),n),C(n,i));
		(i&1)?Dec(res,val):Inc(res,val);
	}else break;
	return res;
}

inline int calc(int d){
	int n=::n/d,m=::m/d;
	if(m<=K)return C(n,m);
	int res=0;
	for(int re i=0;i<=K;++i)
	Inc(res,mul(i+1,Put(m-i,n-m-1,K)));
	return res;
}

signed main(){
#ifdef zxyoi
//	freopen("loop.in","r",stdin);
#endif
	scanf("%d%d%d",&n,&m,&K);
	if(n==m){
		cout<<(K==n)<<"\n";
		return 0;
	}
	if(!m)return cout<<"1\n",0;
	sieves(n);init_inv(n<<1);
	int g=std::__gcd(n,m),ans=0;
	for(int re i=1;i*i<=g;++i)if(g%i==0){
		Inc(ans,mul(phi[i],calc(i)));
		if(i*i!=g)Inc(ans,mul(calc(g/i),phi[g/i]));
	}
	cout<<power(n,mod-2,ans)<<"\n"; 
	return 0;
}