poj1417 带权并查集+0/1背包

题意:有一个岛上住着一些神和魔,并且已知神和魔的数量,现在已知神总是说真话,魔总是说假话,有 n 个询问,问某个神或魔(身份未知),问题是问某个是神还是魔,根据他们的回答,问是否能够确定哪些是神哪些是魔。

对于这些问题,我们只需要发现,如果回答对方是魔,那么即可以判断出这两个不是同一种族,而如果回答对方是神,那么说明这两个是同一种族,那么就可以用带权并查集合并这些神和魔,然后记录两种分别多少个,这样当所有询问都处理完时我们就可以得到一系列的集合,每个集合分别有它的两个种族的人数,但是此时对于每个集合,这两个人数我们并不知道分别哪个是神哪个是魔,这时候就需要用到0/1背包的方法, dp[i][j] 代表处理到第 i 个集合时共 j 个人数(神或魔)的情况数,它从 dp[i-1][j-num1[i]] 和 dp[i-1][j-num2[i]] 转移过来,这里 num1、num2 就是一个并查集的两种人数。转移时顺便记录它是从哪个值转移过来的,便于最后输出。这样只要最后处理完所有集合时正好神的人数或魔的人数的情况正好是一种,那么就说明可行。按照记录的转移遍历回去输出所有解就行。加了个小优化就是当某个集合的两种数量相等的时候,直接可以判断否,因为这时这两个数量等效。

  1 #include<stdio.h>
  2 #include<string.h>
  3 
  4 int fa[605],num[605],n,num1[605],num2[605],p[605],dp[605][605],fat[605][605],num3[605];
  5 int p1,p2;
  6 char s[5];
  7 
  8 int mmax(int a,int b){
  9     return a>b?a:b;
 10 }
 11 
 12 void init(){
 13     for(int i=1;i<=p1+p2;i++){
 14         fa[i]=i;
 15         num1[i]=1;
 16     }
 17     memset(num,0,sizeof(num));
 18     memset(num2,0,sizeof(num2));
 19 }
 20 
 21 int find(int x){
 22     int r=x,t1,t2,c=0;
 23     while(r!=fa[r]){
 24         c+=num[r];
 25         r=fa[r];
 26     }
 27     while(r!=x){
 28         t1=fa[x];
 29         t2=c-num[x];
 30         num[x]=c%2;
 31         fa[x]=r;
 32         x=t1;
 33         c=t2;
 34     }
 35     return r;
 36 }
 37 
 38 int main(){
 39     while(scanf("%d%d%d",&n,&p1,&p2)!=EOF&&n!=0||p1!=0||p2!=0){
 40         int i;
 41         init();
 42         for(i=1;i<=n;i++){
 43             int a,b,v;
 44             scanf("%d%d%s",&a,&b,s);
 45             if(s[0]=='y')v=0;
 46             else v=1;
 47             int x=find(a),y=find(b);
 48             if(x!=y){
 49                 num[x]=((num[b]+v-num[a])%2+2)%2;
 50                 if(num[x]==0){
 51                     num1[y]+=num1[x];
 52                     num2[y]+=num2[x];
 53                 }
 54                 else{
 55                     num1[y]+=num2[x];
 56                     num2[y]+=num1[x];
 57                 }
 58                 fa[x]=y;
 59             }
 60         }
 61         bool f=1;
 62         if(p1==p2)printf("no
");
 63         else{
 64             int cnt=0,j;
 65             for(i=1;i<=p1+p2;i++){
 66                 if(fa[i]==i){
 67                     p[++cnt]=i;
 68                     if(num1[i]==num2[i])f=0;
 69                 }
 70             }
 71             if(f){
 72                 memset(dp,0,sizeof(dp));
 73                 dp[0][0]=1;
 74                 for(i=1;i<=cnt;i++){
 75                     for(j=0;j<=p1;j++){
 76                         if(dp[i-1][j]){
 77                             dp[i][j+num1[p[i]]]+=dp[i-1][j];
 78                             fat[i][j+num1[p[i]]]=j;
 79                             dp[i][j+num2[p[i]]]+=dp[i-1][j];
 80                             fat[i][j+num2[p[i]]]=j;
 81                         }
 82                     }
 83                 }
 84                 if(dp[cnt][p1]!=1)f=0;
 85             }
 86             if(!f)printf("no
");
 87             else{
 88                 int father=p1;
 89                 for(i=cnt;i>=1;i--){
 90                     if(father-fat[i][father]==num1[p[i]]){
 91                         num3[p[i]]=0;
 92                     }
 93                     else num3[p[i]]=1;
 94                     father=fat[i][father];
 95                 }
 96                 for(i=1;i<=p1+p2;i++){
 97                     int x=find(i);
 98                     if(num[i]==num3[x])printf("%d
",i);
 99                 }
100                 printf("end
");
101             }
102         }
103     }
104     return 0;
105 }
View Code
原文地址:https://www.cnblogs.com/cenariusxz/p/4790016.html