【题解】[AGC005D] ~K Perm Counting

题意

如果一个排列 P P P 满足对于所有的 i i i 都有 ∣ P i − i ∣ ≠ k |P_i-i|\neq k Pii=k ,则称排列 P P P 是合法的。求有多少种合法的排列。

Solution:

本题的限制条件很像错排,但是区别在于并不是一一对应的。

考虑容斥。问题转化成了求满足其中 K K K 个条件的排列数。

考虑这样一个序列:1 1+k 1+2k ... 1+mk 。在二分图中是这样的:

请添加图片描述
考虑其中一条链。为了使结构统一,我们把原始排列拆分成 2 n 2n 2n 个点,这样做的好处是将题意限制转化成了 不能选择两条相邻的边 。(这样点和边形成了映射关系)。

d p [ i ] [ j ] [ 0 / 1 ] dp[i][j][0/1] dp[i][j][0/1] 表示处理到前 i i i 个点,其中满足 j j j 个条件,i->i-1 是否连边。状态转移方程为:

  1. d p [ i ] [ j ] [ 0 ] = d p [ i − 1 ] [ j ] [ 0 ] + d p [ i − 1 ] [ j ] [ 1 ] dp[i][j][0]=dp[i-1][j][0]+dp[i-1][j][1] dp[i][j][0]=dp[i1][j][0]+dp[i1][j][1]
  2. d p [ i ] [ j ] [ 1 ] = d p [ i − 1 ] [ j − 1 ] [ 0 ] ( j > = 1 ) dp[i][j][1]=dp[i-1][j-1][0](j>=1) dp[i][j][1]=dp[i1][j1][0](j>=1)

当然可以用 背包2k 条链合并。我们可以把 2k 条链合并在一起进行 d p dp dp ,对于链首的情况要特殊讨论。

时间复杂度 O ( n 2 ) O(n^2) O(n2)

#include<bits/stdc++.h> #define ll long long using namespace std; const int mod=924844033; const int mx=4005; int n,m,k,dp[mx][mx][2],a[mx]; ll res,fac[mx]; int main() { scanf("%d%d",&n,&k); fac[0]=1; for(int i=1;i<=n;i++) fac[i]=fac[i-1]*i%mod; for(int i=1;i<=k;i++) { for(int j=i;j<=n;j+=k) { a[++m]=(j==i); } for(int j=i;j<=n;j+=k) { a[++m]=(j==i); } } dp[0][0][0]=1; for(int i=1;i<=m;i++) { for(int j=0;j<=i;j++) { dp[i][j][0]=(dp[i-1][j][0]+dp[i-1][j][1])%mod; dp[i][j][1]=(a[i]==1||j==0)?0:dp[i-1][j-1][0]; } } for(int i=0;i<=n;i++) { if(i&1) { res=(res-1ll*(dp[m][i][0]+dp[m][i][1])*fac[n-i]%mod)%mod; } else { res=(res+1ll*(dp[m][i][0]+dp[m][i][1])*fac[n-i]%mod)%mod; } } if(res<0) res+=mod; printf("%lld",res); }

__EOF__

本文作者仰望星空的蚂蚁
本文链接https://www.cnblogs.com/cqbzly/p/17530296.html
关于博主:评论和私信会在第一时间回复。或者直接私信我。
版权声明:本博客所有文章除特别声明外,均采用 BY-NC-SA 许可协议。转载请注明出处!
声援博主:如果您觉得文章对您有帮助,可以点击文章右下角推荐一下。您的鼓励是博主的最大动力!
posted @   仰望星空的蚂蚁  阅读(7)  评论(0编辑  收藏  举报  
相关博文:
阅读排行:
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
· DeepSeek 开源周回顾「GitHub 热点速览」
点击右上角即可分享
微信分享提示