天天看點

AtCoder Grand Contest 005 d

題意:

給出n和k,問有多少個長度為n的排列p使得對于任意的i,有 ∣ p i − i ∣ ≠ k |p_i−i|≠k ∣pi​−i∤​=k。

n &lt; = 2000 n&lt;=2000 n<=2000

題解:

考慮容斥,即計算至少有i個位置滿足 ∣ p i − i ∣ = k |p_i-i|=k ∣pi​−i∣=k。

關鍵就是要求出強制 k k k個滿足測方案數。

将位置和數值建成二分圖,相當于求選擇k對合法比對的方案數。

容易發現,将二分圖建出來後連出的邊是一條條鍊,那麼将這些鍊連在一起,dp即可。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
const LL mod=924844033;
LL n,k;
LL f[4010][2010][2];
LL vis[2010],a[4010],N=0;
void add(LL x,LL tim)
{
	LL i=x;a[N+1]=1;
	//printf("N:%d\n",N);
	while(1)
	{
		tim^=1;if(tim) vis[i]=1;N++;
		if(i+k>n) break;i+=k;
	}
}
LL fac[2010];
int main()
{
	scanf("%lld %lld",&n,&k);
	memset(vis,0,sizeof(vis));
	for(LL i=1;i<=n;i++) if(vis[i]==0) add(i,0),add(i,1);
	for(LL i=0;i<=N;i++) f[i][0][0]=1;
	for(LL i=1;i<=N;i++)
		for(LL j=1;j<=n;j++)
		{
			f[i][j][0]=(f[i-1][j][0]+f[i-1][j][1])%mod;
			if(a[i]!=1) f[i][j][1]=f[i-1][j-1][0];
		}
	//for(int i=1;i<=N;i++) printf("%lld ",a[i]);printf("\n");
	//for(LL i=0;i<=n;i++) printf("%lld:%lld\n",i,f[N][i][0]+f[N][i][1]);
	LL ans=0;
	fac[0]=1;for(LL i=1;i<=n;i++) fac[i]=fac[i-1]*i%mod;
	for(LL i=0;i<=n;i++)
		(ans+=(i&1?-1:1)*fac[n-i]%mod*(f[N][i][0]+f[N][i][1])%mod)%=mod;
	printf("%lld",(ans+mod)%mod);
}