Description
Byteotia的领土被占领了,国王Byteasar正在打算组织秘密抵抗运动。国王需要选一些人来进行这场运动,而这些人被分为两部分:一部分成为同谋者活动在被占领区域,另一部分是后勤组织在未被占领的领土上运转。但是这里出现了一个问题: 1. 后勤组织里的任意两人都必须是熟人,以促进合作和提高工作效率。 2. 同谋者的团体中任意两人都不能是熟人。 3. 每一部分都至少要有一个人。国王想知道有多少种分配方案满足以上条件,当然也有可能不存在合理方案。现在国王将这个问题交由你来解决!
Input
第一行一个整数n(2<=n<=5000)表示有n个人参与该抵抗运动,标号为1..n。 之后有n行,第i行的第一个数ki(0<=ki<=n-1)表示i认识ki个人,随后的ki个数表示i的熟人。 p.s.输入满足:如果i是x的熟人,x会在i的序列中出现同时i也会出现在x的熟人序列中。
Output
符合条件的方案总数。
Sample Input
4
2 2 3
2 1 3
3 1 2 4
1 3
Sample Output
HINT
Hint 1和4分到同谋者组织,2和3为后勤组织。 2和4分到同谋者组织,1和3为后勤组织。 4单独分到同谋者组织,1和2、3为后勤组织。
题解:
考虑构造一组可行解,把每个点拆成两个点x0,x1,x0表示后勤组织,x1表示同谋者。
若x与y认识,则x1向y0连边。
若x与y不认识,则x0向y1连边。
如此求出一组2-SAT的可行解,如果无解则答案为0。
若有解,那么最多只能把一个人从后勤组织改为同谋者,也最多只能把一个人从同谋者改为后勤组织,也可以将一个在同谋者另一个在后勤组织的两个人交换。
显然不能直接暴力搞
先预处理出每个点和它相冲突的点的个数和其中任意一个冲突点的编号
冲突点定义为假如把这个点放到对面去,对面的点中会和这个点发生冲突的点
举个例子,假如两个人u和v,u在后勤,v在同谋者,而u,v互相认识,则v是u的矛盾点
或者是两个人a和b,a在同谋者,b在后勤,而a,b不认识,这b是a的矛盾点
先讨论要交换的情况
显然,当一个点的矛盾点数量超过2,他既不能直接去对面,也不能和对面的某个人直接交换
当一个点u的矛盾点为1时,设它的矛盾点为v
假如v也只有一个矛盾点,且v的矛盾点是u,显然这两人可以交换(注意不要重复计算)
假如v没有矛盾点,显然可以直接交换
设在后勤组织中没有矛盾点的个数为t0,同谋者中没有矛盾点的个数为t1,显然这两类点可以任取一对来交换,方案数为t0*t1
再讨论直接去对面的情况
假如这个点没有矛盾点且他所在的组织超过1人(要保证每一边至少有一人),则他可以去对面
然后初始解假如每一边都至少有一人,那么这也是一个合法方案
code:
1 #include<cstdio> 2 #include<iostream> 3 #include<cmath> 4 #include<cstring> 5 #include<algorithm> 6 #define maxn 5005 7 using namespace std; 8 char ch; 9 bool ok; 10 void read(int &x){ 11 for (ok=0,ch=getchar();!isdigit(ch);ch=getchar()) if (ch=='-') ok=1; 12 for (x=0;isdigit(ch);x=x*10+ch-'0',ch=getchar()); 13 if (ok) x=-x; 14 } 15 int n,k,x,con[maxn],cnt0,cnt1,list0[maxn],list1[maxn],num[maxn],boom[maxn]; 16 bool g[maxn][maxn],bo[maxn]; 17 struct Graph{ 18 int tot,now[maxn<<1],son[maxn*maxn],pre[maxn*maxn]; 19 int idx,dfn[maxn<<1],low[maxn<<1],top,stack[maxn],cnt,bel[maxn<<1]; 20 bool in[maxn<<1]; 21 void put(int a,int b){pre[++tot]=now[a],now[a]=tot,son[tot]=b;} 22 void dfs(int u){ 23 dfn[u]=low[u]=++idx,stack[++top]=u,in[u]=1; 24 for (int p=now[u],v=son[p];p;p=pre[p],v=son[p]) 25 if (!dfn[v]) dfs(v),low[u]=min(low[u],low[v]); 26 else if (in[v]) low[u]=min(low[u],dfn[v]); 27 if (dfn[u]==low[u]){ 28 int v; ++cnt; 29 do{v=stack[top--],in[v]=0,bel[v]=cnt;}while (v!=u); 30 } 31 } 32 }G; 33 int main(){ 34 read(n); 35 for (int i=1;i<=n;i++){ 36 read(k),con[i]=k; 37 while (k--) read(x),g[i][x]=1; 38 } 39 for (int i=1;i<=n;i++) for (int j=1;j<=n;j++) if (i!=j){ 40 if (g[i][j]) G.put((i<<1)+1,j<<1); else G.put(i<<1,(j<<1)+1); 41 } 42 for (int i=2;i<=(n<<1)+1;i++) if (!G.dfn[i]) G.dfs(i); 43 for (int i=1;i<=n;i++){ 44 if (G.bel[i<<1]==G.bel[(i<<1)+1]){puts("0");return 0;} 45 if (G.bel[i<<1]<G.bel[(i<<1)+1]) list0[++cnt0]=i; 46 else bo[i]=1,list1[++cnt1]=i; 47 } 48 int ans=(cnt0&&cnt1); 49 for (int i=1;i<=cnt0;i++) for (int j=1;j<=cnt1;j++) 50 if (g[list0[i]][list1[j]]) num[list0[i]]++,boom[list0[i]]=list1[j]; 51 for (int i=1;i<=cnt1;i++) for (int j=1;j<=cnt0;j++) 52 if (!g[list1[i]][list0[j]]) num[list1[i]]++,boom[list1[i]]=list0[j]; 53 for (int i=1;i<=n;i++) if (num[i]==1){ 54 if (num[boom[i]]==1&&boom[i]>i&&boom[boom[i]]==i) ans++; 55 else if (!num[boom[i]]) ans++; 56 } 57 int t0=0,t1=0; 58 for (int i=1;i<=n;i++) if (!num[i]){ 59 if ((bo[i]&&cnt1>1)||(!bo[i]&&cnt0>1)) ans++; 60 if (bo[i]) t1++; else t0++; 61 } 62 ans+=t1*t0; 63 printf("%d ",ans); 64 return 0; 65 } 66