(JZ5230)2019.01.27【NOIP提高组】模拟B组 1.队伍统计

【NOIP2017模拟A组模拟8.5】队伍统计

Description

现在有n个人要排成一列，编号为1->n 。但由于一些不明原因的关系，人与人之间可能存在一些矛盾关系，具体有m条矛盾关系(u,v),表示编号为u的人想要排在编号为v的人前面。要使得队伍和谐，最多不能违背k条矛盾关系（即不能有超过k条矛盾关系(u,v),满足最后v排在了u前面）。问有多少合法的排列。答案对10^9+7取模。

Input

输入文件名为count.in。
第一行包括三个整数n,m,k。
接下来m行，每行两个整数u,v,描述一个矛盾关系(u,v)。
保证不存在两对矛盾关系(u,v),(x,y)，使得u=x且v=y 。

Output

输出文件名为count.out。
输出包括一行表示合法的排列数。

Sample Input
输入1：
4 2 1
1 3
4 2

输入2：
10 12 3
2 6
6 10
1 7
4 1
6 1
2 4
7 6
1 4
10 4
10 9
5 9
8 10

Sample Output
输出1：
18

输出2：
123120

Data Constraint

对于30%的数据，n<=10
对于60%的数据，n<=15
对应100%的数据，n,k<=20,m<=n*(n-1)，保证矛盾关系不重复。

纪中题解：

 状压dp，从前往后逐个确定排列。设? ? [?]表示当前选了?集合这些人，?表示违反了多少矛盾关系。
 暴力的转移复杂度是?(2^??²?)。
 可以发现因为矛盾关系没有重复的，设 ?[?] 表示 要排在? 后面的人的集合，那么每次只要求出? & ?[?]的1的个数就能快速更新?了。
 复杂度?(2^???)

科普：（摘自Tony_Double_Sky的博客）

   状态压缩动态规划，就是我们俗称的状压DP，是利用计算机二进制的性质来描述状态的一种DP方式

位运算：

1. ’&’符号，x&y，会将两个十进制数在二进制下进行与运算，然后返回其十进制下的值。例如3(11)&2(10)=2(10)。
2. ’|’符号，x|y，会将两个十进制数在二进制下进行或运算，然后返回其十进制下的值。例如3(11)|2(10)=3(11)。
3. ’^’符号， x^y，会将两个十进制数在二进制下进行异或运算，然后返回其十进制下的值。例如3(11)^2(10)=1(01)。
4. ’<<’符号，左移操作，x<<2，将x在二进制下的每一位向左移动两位，最右边用0填充，x<<2相当于让x乘以4。相应的，’>>’是右移操作，x>>1相当于给x/2，去掉x二进制下的最有一位。

story↗ time↘(请脑补卡面来打zio豉油（粤语）的音效)：

   一开始写这题的时候，没想到dp。就想着打个全排列，然后判断，水个30分。用C++打着打着，发现不对啊！我不会C++的字符串函数啊！（到后面数字位数会很大，怕数组装不下）就只能中途转回p。
   辛辛苦苦打完之后，发现过不了第二个样例，因为全排列10个数会超时，所以想着只能水10分了。结果结果出来是个0┓( ´∀` )┏问了问别人，据说用C++，再加上剪枝就能水30分，看来在C++的学习中，我还有很长一段路要走呢

#include<cstdio>
using namespace std;
int n,m,k,u,v,ans;
int a[21],g[1<<20],p[21],f[1<<20][21];
inline void dfs(int x,int y,int z){
    if(z>n) return;
    g[x]=y;
    dfs(x,y,z+1);
    dfs(x+p[z],y+1,z+1);
}

int main(){
    freopen("count.in","r",stdin);
    freopen("count.out","w",stdout);
    
    scanf("%d%d%d",&n,&m,&k);
    for(int i=p[0]=1; i<=n; i++) p[i]=p[i-1]<<1;
    dfs(0,0,0);
    for(int i=1; i<=m; i++){
        scanf("%d%d",&u,&v);
        a[v]|=p[u-1];
    }
    f[0][0]=1;
    for(int s=0; s<p[n]; s++)
        for(int j=1; j<=n; j++)
            if(p[j-1]&s){
                int sum=g[s&a[j]];
                for(int i=sum;i<=k;i++) f[s][i]=(f[s][i]+f[s-p[j-1]][i-sum])%1000000007;
            }
    for(int i=0; i<=k ; i++) ans=(ans+f[p[n]-1][i])%1000000007;
    printf("%d",ans);
    return 0;
    
    fclose(stdin); fclose(stdout);
}