NOI2015程式自動分析(並查集)【做題報告】
這是一道國賽水題
題意:(來自洛谷)
題目描述
在實現程式自動分析的過程中,常常需要判定一些約束條件是否能被同時滿足。
考慮一個約束滿足問題的簡化版本:假設x1,x2,x3...代表程式中出現的變數,給定n個形如xi=xj或xi≠xj的變數相等/不等的約束條件,請判定是否可以分別為每一個變數賦予恰當的值,使得上述所有約束條件同時被滿足。例如,一個問題中的約束條件為:x1=x2,x2=x3,x3=x4,x4≠x1,這些約束條件顯然是不可能同時被滿足的,因此這個問題應判定為不可被滿足。
現在給出一些約束滿足問題,請分別對它們進行判定。
輸入輸出格式
輸入格式:
從檔案prog.in中讀入資料。
輸入檔案的第1行包含1個正整數t,表示需要判定的問題個數。注意這些問題之間是相互獨立的。
對於每個問題,包含若干行:
第1行包含1個正整數n,表示該問題中需要被滿足的約束條件個數。接下來n行,每行包括3個整數i,j,e,描述1個相等/不等的約束條件,相鄰整數之間用單個空格隔開。若e=1,則該約束條件為xi=xj;若�e=0,則該約束條件為xi≠xj;
輸出格式:
輸出到檔案 prog.out 中。
輸出檔案包括t行。
輸出檔案的第 k行輸出一個字串“ YES” 或者“ NO”(不包含引號,字母全部大寫),“ YES” 表示輸入中的第k個問題判定為可以被滿足,“ NO” 表示不可被滿足。
輸入輸出樣例
輸入樣例#1: 複製2 2 1 2 1 1 2 0 2 1 2 1 2 1 1輸出樣例#1: 複製
NO YES輸入樣例#2: 複製
2 3 1 2 1 2 3 1 3 1 1 4 1 2 1 2 3 1 3 4 1 1 4 0輸出樣例#2: 複製
YES NO
說明
【樣例解釋1】
在第一個問題中,約束條件為:x1=x2,x1≠x2。這兩個約束條件互相矛盾,因此不可被同時滿足。
在第二個問題中,約束條件為:x1=x2,x1=x2。這兩個約束條件是等價的,可以被同時滿足。
【樣例說明2】
在第一個問題中,約束條件有三個:x1=x2,x2=x3,x3=x1。只需賦值使得x1=x1=x1,即可同時滿足所有的約束條件。
在第二個問題中,約束條件有四個:x1=x2,x2=x3,x3=x4,x4≠x1。由前三個約束條件可以推出x1=x2=x3=x4,然而最後一個約束條件卻要求x1≠x4,因此不可被滿足。
【資料範圍】
【時限2s,記憶體512M】
這道題思路還是很水的,還是帶權並查集,val值1和0進行合併;然後路徑壓縮更新;,十分鐘出程式碼
然後。。。
然後清一下陣列,陣列開大些,然後。。。
這是為什麼呢
看這組數
8 11 0
8 13 0
11 13 0
按我的方式就會WA
所以我後桌大佬的思路就來了
先把所有的相等的進行合併
然後再判定不相等
最後。。
考前摸jx,暴力能AC;
考前摸jx,高一拿省一!
1 #include<cstdio> 2 #include<algorithm> 3 #include<cstring> 4 using namespace std; 5 int T,n,x,y; 6 bool ans=0; 7 int f[101111*2],opt[101111*2],val[101111*2]; 8 struct node{ 9 int val,hs,no; 10 }nd[101111*2]; 11 bool cmp1(node a,node b) 12 {return a.val<b.val;} 13 bool cmp2(node a,node b) 14 {return a.no<b.no;} 15 int find(int x) 16 { 17 if(f[x]==x)return x; 18 int ff=find(f[x]); 19 val[x]=(val[x]+val[f[x]])%2; 20 return f[x]=ff; 21 } 22 int main() 23 { 24 //freopen("testdata.in","r",stdin); 25 scanf("%d",&T); 26 while(T--) 27 { 28 memset(val,0,sizeof(val)); 29 ans=0; 30 scanf("%d",&n); 31 for(int i=1;i<=n;i++) 32 { 33 scanf("%d%d%d",&x,&y,&opt[i]); 34 nd[i].val=x; 35 nd[i+n].val=y; 36 nd[i].no=nd[i+n].no=i; 37 } 38 sort(nd+1,nd+2*n+1,cmp1); 39 int tt=0; 40 for(int i=1;i<=2*n;i++) 41 { 42 if(nd[i].val==nd[i-1].val)nd[i].hs=tt; 43 else nd[i].hs=++tt; 44 f[tt]=tt; 45 } 46 sort(nd+1,nd+2*n+1,cmp2); 47 int cnt=0; 48 for(int i=1;i<=2*n;i+=2) 49 { 50 cnt=i/2+1; 51 if(opt[cnt]==0)continue; 52 x=nd[i].hs,y=nd[i+1].hs; 53 int fx=find(x),fy=find(y); 54 if(fx==fy) 55 { 56 continue; 57 }else 58 { 59 val[fx]=(val[x]+val[y])%2; 60 f[fx]=fy; 61 } 62 } 63 for(int i=1;i<=2*n;i+=2) 64 { 65 cnt=i/2+1; 66 if(opt[cnt]==1)continue; 67 x=nd[i].hs,y=nd[i+1].hs; 68 int fx=find(x),fy=find(y); 69 if(fx==fy) 70 { 71 if(val[x]==val[y]) 72 { 73 ans=1; 74 break; 75 } 76 } 77 } 78 if(ans==1)printf("NO\n"); 79 else printf("YES\n"); 80 } 81 return 0; 82 }