codeforces gym 100548F. Color(容斥原理)

传送门
题意简述:对n个排成一排的物品涂色,有m种颜色可选。
要求相邻的物品颜色不相同,且总共恰好有K种颜色,问所有可行的方案数。( n , m ≤ 1 e 9 , k ≤ 1 e 6 n,m\le1e9,k\le1e6 n,m1e9,k1e6


思路:
容斥原理套路:
先不考虑是否选全 k k k种颜色,方案数为 C m k ∗ k ∗ ( k − 1 ) n − 1 C_m^k*k*(k-1)^{n-1} Cmkk(k1)n1
然后枚举剩下的至少有几种颜色没选来容斥掉非法情况:
于是 A n s = C m k ∑ i = k 1 ( − 1 ) k − i C k i i ( i − 1 ) n − 1 Ans=C_m^k\sum_{i=k}^1(-1)^{k-i}C_k^ii(i-1)^{n-1} Ans=Cmki=k1(1)kiCkii(i1)n1
代码:

#include<bits/stdc++.h>
#define ri register int
using namespace std;
const int N=1e6+5,mod=1e9+7;
typedef long long ll;
inline int read(){
	int ans=0;
	char ch=getchar();
	while(!isdigit(ch))ch=getchar();
	while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
	return ans;
}
inline int add(const int&a,const int&b){return a+b>=mod?a+b-mod:a+b;}
inline int dec(const int&a,const int&b){return a>=b?a-b:a-b+mod;}
inline int mul(const int&a,const int&b){return (ll)a*b%mod;}
inline int ksm(int a,int p){int ret=1;for(;p;p>>=1,a=mul(a,a))if(p&1)ret=mul(ret,a);return ret;}
int n,m,k,C[N],inv[N],mult;
inline void init(){
	inv[1]=1;
	for(ri i=2;i<=1000000;++i)inv[i]=mul(inv[mod-mod/i*i],mod-mod/i);
}
inline void Init(){
	C[0]=1,mult=1;
	for(ri i=1;i<=k;++i)C[i]=mul(mul(C[i-1],k-i+1),inv[i]);
	for(ri i=1;i<=k;++i)mult=mul(mult,mul(m-i+1,inv[i]));
}
int main(){
	init();
	for(ri ans,tt=1,up=read();tt<=up;++tt){
		n=read(),m=read(),k=read(),ans=0;
		Init();
		for(ri i=k,tmp;i;--i){
			tmp=mul(C[i],mul(i,ksm(i-1,n-1)));
			ans=(k-i)&1?dec(ans,tmp):add(ans,tmp);
		}
		cout<<"Case #"<<tt<<": "<<mul(mult,ans)<<'\n';
	}
	return 0;
}
  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值