【容斥原理】【DP】AGC005D ~K Perm Counting

最新推荐文章于 2021-08-09 09:45:41 发布

616156

最新推荐文章于 2021-08-09 09:45:41 发布

阅读量347

点赞数

CC 4.0 BY-SA版权

分类专栏：容斥原理 DP

本文链接：https://blog.youkuaiyun.com/qq_34454069/article/details/82790017

DP 同时被 2 个专栏收录

100 篇文章

订阅专栏

容斥原理

9 篇文章

订阅专栏

本文探讨了在解决复杂组合计数问题时，如何运用容斥原理与动态规划技巧。通过一个具体的算法实现案例，详细解析了如何将互相冲突的位置进行分类，并利用动态规划求解特定条件下方案的数量。最后，通过迭代计算得出最终的答案。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

分析：

比较简单（板）的容斥题。

设枚举出i个非法位置的方案数为 $f_i$
答案就是 $∑i=0i≤n(−1)ifi∗(n−i)!\sum_{i=0}^{i\leq n}(-1)^if_i*(n-i)!$

问题就在于如何求 $f_i$ 。显然，可以把互相可能矛盾的位置分为一类，然后每一类搞一个 $n^2$ 的DP求出来有k个矛盾的方案数。然后所有的一起求个背包。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define SF scanf
#define PF printf
#define MAXN 2010
#define MOD 924844033
using namespace std;
typedef long long ll;
ll dp[MAXN][MAXN][2];
void prepare(){
	dp[0][0][0]=1;
	dp[0][0][1]=1;
	for(int i=1;i<=2000;i++)
		for(int j=0;j<=i;j++){
			dp[i][j][0]=dp[i-1][j][1];
			if(i>=1&&j>=1)
				dp[i][j][0]+=dp[i-1][j-1][0];
			dp[i][j][0]%=MOD;
			
			dp[i][j][1]=dp[i][j][0];
			if(i>=1&&j>=1)
				dp[i][j][1]+=dp[i-1][j-1][1];
			dp[i][j][1]%=MOD;
		}
}
int n,k;
ll sum[MAXN],fac[MAXN];
int main(){
	prepare();
	SF("%d%d",&n,&k);
	sum[0]=1;
	for(int i=1;i<=k*2&&i<=n;i++){
		int lenx=(n-i)/(k*2)+1;
		int fir=i-k;
		int las=(lenx-1)*k*2+i+k;
		int kd;
		if(fir>=1&&las<=n)
			kd=1;
		else if(fir>=1||las<=n)
			kd=0;
		else{
			kd=1;
			lenx--;
		}
		/*PF("%d:{%d %d (%d %d)}\n",i,lenx,kd,fir,las);
		for(int j=0;j<=lenx;j++){
			PF("[%d %lld]\n",j,dp[lenx][j][kd]);
		}*/
		for(int i=n;i>=0;i--)
			for(int j=1;j<=lenx&&i-j>=0;j++){
				sum[i]+=sum[i-j]*dp[lenx][j][kd];
				sum[i]%=MOD;
			}
		/*PF("now,sum:\n");
		for(int i=1;i<=n;i++)
			PF("{%d %lld}\n",i,sum[i]);
		PF("-----------\n");*/
	}
	fac[0]=1;
	for(int i=1;i<=n;i++)
		fac[i]=fac[i-1]*i%MOD;
	ll ans=0;
	for(int i=0;i<=n;i++){
		if(i%2==0)
			ans=(ans+sum[i]*fac[n-i])%MOD;
		else
			ans=(ans-sum[i]*fac[n-i])%MOD;
	}
	PF("%lld",(ans+MOD)%MOD);
}