题目链接

\(Description\)

给定\(n,k\),求 满足对于所有\(i\),\(|a_i-i|\neq k\)的排列的个数。

\(2\leq n\leq 2000,\quad 1\leq k\leq n-1\)。

\(Solution\)

容斥。则\(Ans=\sum_{i=0}^n(-1)^ig(i)(n-i)!\),其中\(g(i)\)为至少有\(i\)个位置满足\(|a_i-i|=k\)的排列数。

考虑如何计算\(g(x)\)。每个\(i\)向\(i+k\)和\(i-k\)连边,可以得到一张二分图,\(g(x)\)就是在这张二分图上选\(x\)个匹配的方案数。

我们还可以发现,图中的匹配形成了\(2k\)条互不相交的链,且每条链上的数模\(k\)相同(也就是模\(k\)不同的数是互不影响的,所以枚举模数就可以得到所有链了)。

如果只有一条长\(l\)的链,那么就是对\(l-1\)个点DP,\(f[i][j][0/1]\)表示当前第\(i\)个点,已经选了\(j\)个匹配,这个点选不选(这个点选要求上一个点没选)。

对于\(2k\)条链也是一样的,只要在链之间加一个点,并强制它不能选(\(f[i][j][1]=0\)),就可以把这些链合在一起DP啦。

这样复杂度\(O(n^2)\),然而好多dalao用NTT\(n\log n\)过掉了orz。


因为想跑快点(个人习惯)写的有点乱,代码也可以看他的:http://www.cnblogs.com/wxjor/p/9476998.html。

//20ms	256KB
#include <cstdio>
#include <algorithm>
#define mod 924844033
typedef long long LL;
const int N=4005;//2n int ban[N],f[2][N][2],fac[N]; int main()
{
int n,K,cnt=0; scanf("%d%d",&n,&K);
for(int i=1; i<=K; ++i)
{
ban[++cnt]=1;
for(int j=i+K; j<=n; j+=K) ban[++cnt]=0;
ban[++cnt]=1;//就是两条链啊
for(int j=i+K; j<=n; j+=K) ban[++cnt]=0;
}
int p=0; f[p][0][0]=1;
for(int i=0; i<cnt; ++i,p^=1)
for(int j=0; j<=i; ++j)
{
f[p^1][j][0]=(f[p][j][0]+f[p][j][1])%mod;
if(!ban[i+1]) f[p^1][j+1][1]=f[p][j][0];
f[p][j][0]=f[p][j][1]=0;
}
fac[0]=1;
for(int i=1; i<=n; ++i) fac[i]=1ll*fac[i-1]*i%mod;
LL ans=0;
for(int i=0; i<=n; ++i)
ans+=i&1?mod-1ll*fac[n-i]*(f[p][i][0]+f[p][i][1])%mod:1ll*fac[n-i]*(f[p][i][0]+f[p][i][1])%mod;
printf("%lld\n",ans%mod); return 0;
}

AGC 005D.~K Perm Counting(容斥 DP 二分图)的更多相关文章

  1. AGC 005 D - ~K Perm Counting

    D - ~K Perm Counting 链接 题意: 求有多少排列对于每个位置i都满足$|ai−i|!=k$.n<=2000 分析: 容斥+dp. $answer = \sum\limits_ ...

  2. [Agc005D]K Perm Counting

    [Agc005D] K Perm Counting Description 糟糕爷特别喜爱排列.他正在构造一个长度为N的排列.但是他特别讨厌正整数K.因此他认为一个排列很糟糕,当且仅当存在至少一个i( ...

  3. HDU 5794 A Simple Chess (容斥+DP+Lucas)

    A Simple Chess 题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5794 Description There is a n×m board ...

  4. 【BZOJ3622】已经没有什么好害怕的了 容斥+DP

    [BZOJ3622]已经没有什么好害怕的了 Description Input Output Sample Input 4 2 5 35 15 45 40 20 10 30 Sample Output ...

  5. BZOJ3622(容斥+dp)

    思路 "恰k个"考虑求至少k.k+1--个容斥 题面说所有数字都不同,可以将所求转化为糖比药多的组数恰为\((n+k)/2\)的方案数 \(f[i][j]\)数组我觉得更好的理解方 ...

  6. [CF1086E]Beautiful Matrix(容斥+DP+树状数组)

    给一个n*n的矩阵,保证:(1)每行都是一个排列 (2)每行每个位置和上一行对应位置不同.求这个矩阵在所有合法矩阵中字典序排第几.考虑类似数位DP的做法,枚举第几行开始不卡限制,那么显然之前的行都和题 ...

  7. $bzoj2560$ 串珠子 容斥+$dp$

    正解:容斥+$dp$ 解题报告: 传送门$QwQ$ $umm$虽然题目蛮简练的了但还是有点难理解,,,我再抽象一点儿,就说有$n$个点,点$i$和点$j$之间有$a_{i,j}$条无向边可以连,问有多 ...

  8. [AGC005D] ~K Perm Counting [dp]

    题面 传送门 思路 首先可以明确的一点是,本题中出现不满足条件的所有的数,都是分组的 只有模$K$意义下相同的数之间才会出现不满足条件的情况,而且仅出现在相邻的情况 那么我们考虑把这个性质利用起来 我 ...

  9. 【XSY3156】简单计数II 容斥 DP

    题目大意 定义一个序列的权值为:把所有相邻的相同的数合并为一个集合后,所有集合的大小的乘积. 特别的,第一个数和最后一个数是相邻的. 现在你有 \(n\) 种数,第 \(i\) 种有 \(c_i\) ...

随机推荐

  1. pyhon----函数和方法的关系

    1.如果使用类名调用,为函数,需要手动传self 2.如果使用对象调用,为方法,不用手动传self class Foo(object): def __init__(self): self.name=& ...

  2. cf1140E 回文串+染色方案dp

    有点硬核的dp..要用到一个结论.. /* 把原串拆成奇偶串,再拆成极大连续的-1串:该串两端都是非-1数,中间都是-1,并且下标要么都是偶数,要么都是技术 然后对所有这些串进行dp,dp[i][0] ...

  3. hdu4009最小树形图板子题

    /*调了一下午的最小树形图,昨天刚刚看懂模板..最小树形图,就是有向图的最小生成树,很神奇==*/ #include<iostream> #include<cstring> # ...

  4. GAN-生成手写数字-Keras

    from keras.models import Sequential from keras.layers import Dense from keras.layers import Reshape ...

  5. Python函数系列之eval()

    1.作用:将字符串str当成有效的表达式来求值并返回计算结果. 2.语法:eval(source[, globals[, locals]])  3.说明:参数:source:一个Python表达式或函 ...

  6. Leetcode刷题第20天

    一.找树左下角的值 题目:513. Find Bottom Left Tree Value C++ Soution 1: /** * Definition for a binary tree node ...

  7. 记录mysql正在执行的SQL语句

    show variables like "general_log%"; SET GLOBAL general_log = 'ON';SET GLOBAL general_log = ...

  8. Linux(centos)系统各个目录的作用详解 推荐

    文件系统的类型 LINUX有四种基本文件系统类型:普通文件.目录文件.连接文件和特殊文件,可用file命令来识别. 普通文件:如文本文件.C语言元代码.SHELL脚本.二进制的可执行文件等,可用cat ...

  9. Task任务的屏障机制

    Barrier 是 .Net 提供的一直并发的机制,它允许多个任务同步他们不同阶段的并发工作. 这里的关键点是[多个任务]和[不同阶段]. 假设有4个相同的任务(Task),每个任务都有4个阶段(Ph ...

  10. Node 杂技

    1.关于require 当文件夹a中含有index.js时,在b.js中如果有require("文件夹a的路径"),则将会自动执行index.js的语句