【2-sat】专题~
2-sat是一个逻辑性很强的算法,但是其套路比较固定,所以不是很热,题目很少,但也不乏AC后让人大呼爽快的好题,下面放出两篇极品论文还有几道题目的题解以供交流~~
2-sat学习:
强连通判是否存在解:
HDU 3062 Party [入门]
2-sat入门题,建图后求强连通判断是否有解即可
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<cstdio>
2 #include<cstring>
3 usingnamespace std;
4 #define MIN(x,y) (x<y?x:y)
5 #define N 2005
6
7 struct Edge{
8 int vtx,next;
9 }E[999999];
10 int head[N],size;
11 int dfn[N],low[N],blg[N],step,scc;
12 int stk[N],top;
13 bool ins[N];
14
15 void Init(){
16 memset(head,-1,sizeof(head));
17 memset(dfn,-1,sizeof(dfn));
18 memset(ins,false,sizeof(ins));
19 size=step=scc=0; top=-1;
20 }
21
22 void Insert(int u,int v){
23 E[size].vtx=v;
24 E[size].next=head[u];
25 head[u]=size++;
26 }
27
28 void Tarjan(int u){
29 stk[++top]=u; ins[u]=true;
30 dfn[u]=low[u]=step++;
31 for(int i=head[u];i!=-1;i=E[i].next){
32 int v=E[i].vtx;
33 if(dfn[v]==-1){
34 Tarjan(v);
35 low[u]=MIN(low[u],low[v]);
36 }elseif(ins[v]){
37 low[u]=MIN(low[u],dfn[v]);
38 }
39 }
40 if(low[u]==dfn[u]){
41 for(int v=-1;v!=u;top--){
42 v=stk[top];
43 ins[v]=false;
44 blg[v]=scc;
45 }
46 scc++;
47 }
48 }
49
50 int main(){
51 int n,m;
52 while(~scanf("%d%d",&n,&m)){
53 Init();
54 for(int i=0;i<m;i++){
55 int a1,a2,c1,c2;
56 scanf("%d%d%d%d",&a1,&a2,&c1,&c2);
57 Insert((a1<<1)+c1,(a2<<1|1)-c2);
58 Insert((a2<<1)+c2,(a1<<1|1)-c1);
59 }
60 for(int i=0;i<2*n;i++){
61 if(dfn[i]==-1) Tarjan(i);
62 }
63 bool flag=true;
64 for(int i=0;i<n;i++){
65 if(blg[i<<1]==blg[i<<1|1]){
66 flag=false;
67 break;
68 }
69 }
70 puts(flag?"YES":"NO");
71 }
72 }
HDU 1824 Let's go home [基础]
把每队的两名队员当作一个整体,队长和队员为相对状态,建图求强连通即可
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<cstdio>
2 #include<cstring>
3 #include<cstdlib>
4 usingnamespace std;
5 #define MIN(x,y) (x>y?y:x)
6 #define N 2005
7
8 struct Edge{
9 int vtx,next;
10 }E[999999];
11 int head[N],size;
12 int dfn[N],low[N],blg[N],step,scc;
13 int stk[N],top;
14 bool ins[N];
15
16 void Init(){
17 memset(head,-1,sizeof(head));
18 memset(dfn,-1,sizeof(dfn));
19 memset(ins,false,sizeof(ins));
20 size=step=scc=0; top=-1;
21 }
22
23 void Insert(int u,int v){
24 E[size].vtx=v;
25 E[size].next=head[u];
26 head[u]=size++;
27 }
28
29 void Tarjan(int u){
30 stk[++top]=u; ins[u]=true;
31 dfn[u]=low[u]=step++;
32 for(int i=head[u];i!=-1;i=E[i].next){
33 int v=E[i].vtx;
34 if(dfn[v]==-1){
35 Tarjan(v);
36 low[u]=MIN(low[u],low[v]);
37 }elseif(ins[v]){
38 low[u]=MIN(low[u],dfn[v]);
39 }
40 };
41 if(low[u]==dfn[u]){
42 for(int v=-1;v!=u;top--){
43 v=stk[top];
44 ins[v]=false;
45 blg[v]=scc;
46 }
47 scc++;
48 }
49 }
50
51 int B[3005];
52 int main(){
53 int n,m;
54 while(~scanf("%d%d",&n,&m)){
55 for(int i=0;i<n;i++){
56 int l,s1,s2;
57 scanf("%d%d%d",&l,&s1,&s2);
58 B[l]=i<<1;
59 B[s1]=i<<1|1;
60 B[s2]=i<<1|1;
61 }
62 Init();
63 for(int i=0;i<m;i++){
64 int a,b;
65 scanf("%d%d",&a,&b);
66 Insert(B[a],B[b]^1);
67 Insert(B[b],B[a]^1);
68 }
69 for(int i=0;i<n*2;i++){
70 if(dfn[i]==-1) Tarjan(i);
71 }
72 bool flag=true;
73 for(int i=0;i<n;i++){
74 if(blg[i<<1]==blg[i<<1|1]){
75 puts("no");
76 main();
77 exit(0);
78 }
79 }
80 puts("yes");
81 }
82 }
POJ 3207 Ikki’s Story IV - Panda’s Trick [基础]
判断两段是否相交,如果相交的话不能在同侧,在两侧为相对状态,建图求强连通
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<cstdio>
2 #include<cstring>
3 usingnamespace std;
4 #define MIN(x,y) (x>y?y:x)
5 #define N 1005
6
7 struct Edge{
8 int vtx,next;
9 Edge(){}
10 Edge(int v,int n):
11 vtx(v),next(n){}
12 }E[999999];
13 int head[N],size;
14 int dfn[N],low[N],blg[N],step,scc;
15 int stk[N],top;
16 bool ins[N];
17
18 void Init(){
19 memset(head,-1,sizeof(head));
20 memset(dfn,-1,sizeof(dfn));
21 memset(ins,false,sizeof(ins));
22 size=step=scc=0; top=-1;
23 }
24
25 void Insert(int u,int v){
26 E[size]=Edge(v,head[u]);
27 head[u]=size++;
28 }
29
30 void Tarjan(int u){
31 stk[++top]=u; ins[u]=true;
32 dfn[u]=low[u]=step++;
33 for(int i=head[u];i!=-1;i=E[i].next){
34 int v=E[i].vtx;
35 if(dfn[v]==-1){
36 Tarjan(v);
37 low[u]=MIN(low[u],low[v]);
38 }elseif(ins[v]){
39 low[u]=MIN(low[u],dfn[v]);
40 }
41 }
42 if(low[u]==dfn[u]){
43 for(int v=-1;v!=u;top--){
44 v=stk[top];
45 ins[v]=false;
46 blg[v]=scc;
47 }
48 scc++;
49 }
50 }
51
52 int l[505],r[505];
53 int main(){
54 int n,m;
55 while(~scanf("%d%d",&n,&m)){
56 for(int i=0;i<m;i++){
57 scanf("%d%d",&l[i],&r[i]);
58 if(l[i]>r[i]){
59 int tmp=l[i]; l[i]=r[i]; r[i]=tmp;
60 }
61 }
62 Init();
63 for(int i=0;i<m;i++){
64 for(int j=i+1;j<m;j++){
65 int cnt=0;
66 if(l[i]>l[j]&&l[i]<r[j]) cnt++;
67 if(r[i]>l[j]&&r[i]<r[j]) cnt++;
68 if(cnt==1){
69 Insert(i<<1,j<<1|1);
70 Insert(j<<1|1,i<<1);
71 Insert(i<<1|1,j<<1);
72 Insert(j<<1,i<<1|1);
73 }
74 }
75 }
76 for(int i=0;i<2*m;i++){
77 if(dfn[i]==-1) Tarjan(i);
78 }
79 bool flag=true;
80 for(int i=0;i<m;i++){
81 if(blg[i<<1]==blg[i<<1|1]){
82 flag=false;
83 break;
84 }
85 }
86 puts(flag?"panda is telling the truth...":"the evil panda is lying again");
87 }
88 }
POJ 3678 Katu Puzzle [中等]
2-sat建图题,把每个值是1(a)和0(~a)为两种状态,分清楚各种操作的本质就很简单了
AND 结果为1:建边 ~x->x,~y->y (两个数必须全为1)
AND 结果为0:建边 y->~x,x->~y (两个数至少有一个为0)
OR 结果为1:建边 ~x->y,~y->x (两个数至少有一个为1)
OR 结果为0:建边 x->~x,y->~y (两个数必须全为0)
XOR 结果为1:建边 x->~y,y->~x,~y->x,~x->y (两个数必须不同)
XOR 结果为0:建边 x->y,y->x,~x->~y,~y->~x (两个数必须相同)
然后求强联通即可
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<cstdio>
2 #include<cstring>
3 usingnamespace std;
4 #define MIN(x,y) (x>y?y:x)
5 #define N 2005
6
7 struct Edge{
8 int vtx,next;
9 Edge(){}
10 Edge(int v,int n):
11 vtx(v),next(n){}
12 }E[999999];
13 int head[N],size;
14 int dfn[N],low[N],blg[N],step,scc;
15 int stk[N],top;
16 bool ins[N];
17
18 void Init(){
19 memset(head,-1,sizeof(head));
20 memset(dfn,-1,sizeof(dfn));
21 memset(ins,false,sizeof(ins));
22 size=step=scc=0; top=-1;
23 }
24
25 void Insert(int u,int v){
26 E[size]=Edge(v,head[u]);
27 head[u]=size++;
28 }
29
30 void Tarjan(int u){
31 stk[++top]=u; ins[u]=true;
32 dfn[u]=low[u]=step++;
33 for(int i=head[u];i!=-1;i=E[i].next){
34 int v=E[i].vtx;
35 if(dfn[v]==-1){
36 Tarjan(v);
37 low[u]=MIN(low[u],low[v]);
38 }elseif(ins[v]){
39 low[u]=MIN(low[u],dfn[v]);
40 }
41 }
42 if(low[u]==dfn[u]){
43 for(int v=-1;v!=u;top--){
44 v=stk[top];
45 ins[v]=false;
46 blg[v]=scc;
47 }
48 scc++;
49 }
50 }
51
52 int main(){
53 int n,m;
54 while(~scanf("%d%d",&n,&m)){
55 Init();
56 for(int i=0;i<m;i++){
57 int a,b,c;
58 char cmd[5];
59 scanf("%d%d%d %s",&a,&b,&c,cmd);
60 if(cmd[0]=='A'){
61 if(c==1){
62 Insert(a<<1|1,a<<1);
63 Insert(b<<1|1,b<<1);
64 }else{
65 Insert(a<<1,b<<1|1);
66 Insert(b<<1,a<<1|1);
67 }
68 }elseif(cmd[0]=='O'){
69 if(c==1){
70 Insert(a<<1|1,b<<1);
71 Insert(b<<1|1,a<<1);
72 }else{
73 Insert(a<<1,a<<1|1);
74 Insert(b<<1,b<<1|1);
75 }
76 }elseif(cmd[0]=='X'){
77 if(c==1){
78 Insert(a<<1,b<<1|1);
79 Insert(a<<1|1,b<<1);
80 Insert(b<<1,a<<1|1);
81 Insert(b<<1|1,a<<1);
82 }else{
83 Insert(a<<1,b<<1);
84 Insert(a<<1|1,b<<1|1);
85 Insert(b<<1,a<<1);
86 Insert(b<<1|1,a<<1|1);
87 }
88 }
89 }
90 for(int i=0;i<n*2;i++){
91 if(dfn[i]==-1) Tarjan(i);
92 }
93 bool flag=true;
94 for(int i=0;i<n;i++){
95 if(blg[i<<1]==blg[i<<1|1]){flag=false;break;}
96 }
97 puts(flag?"YES":"NO");
98 }
99 }
二分+强联通判解:
HDU 3622 Bomb Game [中等]
比较裸的 二分+2-sat
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<cstdio>
2 #include<cstring>
3 usingnamespace std;
4 #define MIN(x,y) (x>y?y:x)
5 #define N 205
6
7 struct Edge{
8 int vtx,next;
9 }E[999999];
10 int head[N],size;
11 int dfn[N],low[N],blg[N],step,scc;
12 int stk[N],top;
13 bool ins[N];
14
15 void Init(){
16 memset(head,-1,sizeof(head));
17 memset(dfn,-1,sizeof(dfn));
18 memset(ins,false,sizeof(ins));
19 size=step=scc=0; top=-1;
20 }
21
22 void Insert(int u,int v){
23 E[size].vtx=v;
24 E[size].next=head[u];
25 head[u]=size++;
26 }
27
28 void Tarjan(int u){
29 stk[++top]=u; ins[u]=true;
30 dfn[u]=low[u]=step++;
31 for(int i=head[u];i!=-1;i=E[i].next){
32 int v=E[i].vtx;
33 if(dfn[v]==-1){
34 Tarjan(v);
35 low[u]=MIN(low[u],low[v]);
36 }elseif(ins[v]){
37 low[u]=MIN(low[u],dfn[v]);
38 }
39 }
40 if(low[u]==dfn[u]){
41 for(int v=-1;v!=u;top--){
42 v=stk[top];
43 ins[v]=false;
44 blg[v]=scc;
45 }
46 scc++;
47 }
48 }
49
50 inline double len(int ax,int ay,int bx,int by){
51 return (double)((ax-bx)*(ax-bx)+(ay-by)*(ay-by));
52 }
53
54 int lx[105],ly[105],rx[105],ry[105];
55 int main(){
56 int n;
57 while(~scanf("%d",&n)){
58 for(int i=0;i<n;i++){
59 scanf("%d%d%d%d",&lx[i],&ly[i],&rx[i],&ry[i]);
60 }
61 double low=0.0,high=10000.0;
62 double mid=(low+high)/2.0,ans=mid;
63 while(high-low>1e-4){
64 Init();
65 double ll=mid*mid;
66 for(int i=0;i<n;i++){
67 for(int j=0;j<n;j++){
68 if(i==j) continue;
69 if(len(lx[i],ly[i],lx[j],ly[j])<ll){
70 Insert(i<<1,j<<1|1);
71 Insert(j<<1,i<<1|1);
72 }
73 if(len(rx[i],ry[i],rx[j],ry[j])<ll){
74 Insert(i<<1|1,j<<1);
75 Insert(j<<1|1,i<<1);
76 }
77 if(len(lx[i],ly[i],rx[j],ry[j])<ll){
78 Insert(i<<1,j<<1);
79 Insert(j<<1|1,i<<1|1);
80 }
81 if(len(rx[i],ry[i],lx[j],ly[j])<ll){
82 Insert(i<<1|1,j<<1|1);
83 Insert(j<<1,i<<1);
84 }
85 }
86 }
87 for(int i=0;i<(n<<1);i++){
88 if(dfn[i]==-1) Tarjan(i);
89 }
90 bool flag=true;
91 for(int i=0;i<n;i++){
92 if(blg[i<<1]==blg[i<<1|1]){flag=false;break;}
93 }
94 if(flag){
95 low=(ans=mid);
96 }else{
97 high=mid;
98 }
99 mid=(low+high)/2.0;
100 }
101 printf("%.2lf\n",ans/2);
102 }
103 }
POJ 2296 Map Labeler [中等]
几种情况讨论一下在上边(a),在下边(~a),假设点i在点j的上
如果|xi-xj|>=mid continue;
否则|yi-yj|>=2*mid continue;
否则|yi-yj|>=mid 建边 j->i,~i->~j;
否则|yi-yj|>0 建边 ~i->i,j->~j;
否则|yi-yj|==0 建边 i->~j,~i->j,j->~i,~j->i;
等于0的情况很容易忘记考虑
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<cstdio>
2 #include<cstring>
3 usingnamespace std;
4 #define MIN(x,y) (x>y?y:x)
5 #define N 205
6
7 struct Edge{
8 int vtx,next;
9 }E[999999];
10 int head[N],size;
11 int dfn[N],low[N],blg[N],step,scc;
12 int stk[N],top;
13 bool ins[N];
14
15 void Init(){
16 memset(head,-1,sizeof(head));
17 memset(dfn,-1,sizeof(dfn));
18 memset(ins,false,sizeof(ins));
19 scc=step=size=0; top=-1;
20 }
21
22 void Insert(int u,int v){
23 E[size].vtx=v;
24 E[size].next=head[u];
25 head[u]=size++;
26 }
27
28 void Tarjan(int u){
29 stk[++top]=u; ins[u]=true;
30 dfn[u]=low[u]=step++;
31 for(int i=head[u];i!=-1;i=E[i].next){
32 int v=E[i].vtx;
33 if(dfn[v]==-1){
34 Tarjan(v);
35 low[u]=MIN(low[u],low[v]);
36 }elseif(ins[v]){
37 low[u]=MIN(low[u],dfn[v]);
38 }
39 }
40 if(low[u]==dfn[u]){
41 for(int v=-1;v!=u;top--){
42 v=stk[top];
43 ins[v]=false;
44 blg[v]=scc;
45 }
46 scc++;
47 }
48 }
49
50 inline int ABS(int x){
51 if(x<0) return-x;
52 return x;
53 }
54 int px[N],py[N];
55 int main(){
56 int t;
57 scanf("%d",&t);
58 while(t--){
59 int n;
60 scanf("%d",&n);
61 for(int i=0;i<n;i++){
62 scanf("%d%d",&px[i],&py[i]);
63 }
64 int low=0,high=10000;
65 int mid=(low+high)>>1,ans=mid;
66 while(low<=high){
67 Init();
68 for(int i=0;i<n;i++){
69 for(int j=i+1;j<n;j++){
70 if(ABS(px[i]-px[j])>=mid) continue;
71 if(ABS(py[i]-py[j])>=(mid<<1)) continue;
72 elseif(ABS(py[i]-py[j])>=mid){
73 if(py[i]>py[j]){
74 Insert(i<<1,j<<1);
75 Insert(j<<1|1,i<<1|1);
76 }else{
77 Insert(i<<1|1,j<<1|1);
78 Insert(j<<1,i<<1);
79 }
80 }else{
81 if(py[i]<py[j]){
82 Insert(i<<1|1,i<<1);
83 Insert(j<<1,j<<1|1);
84 }elseif(py[i]>py[j]){
85 Insert(i<<1,i<<1|1);
86 Insert(j<<1|1,j<<1);
87 }else{
88 Insert(i<<1,j<<1|1);
89 Insert(i<<1|1,j<<1);
90 Insert(j<<1,i<<1|1);
91 Insert(j<<1|1,i<<1);
92 }
93 }
94 }
95 }
96 for(int i=0;i<(n<<1);i++){
97 if(dfn[i]==-1) Tarjan(i);
98 }
99 bool flag=true;
100 for(int i=0;i<n;i++){
101 if(blg[i<<1]==blg[i<<1|1]){
102 flag=false; break;
103 }
104 }
105 if(flag){
106 low=(ans=mid)+1;
107 }else{
108 high=mid-1;
109 }
110 mid=(low+high)>>1;
111 }
112 printf("%d\n",ans);
113 }
114 }
HDU 3715 Go Deeper [中等]
x值为1(a)和0(~a)两种情况
c=0时,建边 ~a->b,~b->a
c=1时,建边 a->b,~a->~b,b->a,~b->~a
c=2时,建边 a->~b,b->~a
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<cstdio>
2 #include<cstring>
3 usingnamespace std;
4 #define MIN(x,y) (x>y?y:x)
5 #define N 405
6
7 struct Edge{
8 int vtx,next;
9 }E[999999];
10 int head[N],size;
11 int dfn[N],low[N],blg[N],step,scc;
12 int stk[N],top;
13 bool ins[N];
14
15 void Init(){
16 memset(head,-1,sizeof(head));
17 memset(dfn,-1,sizeof(dfn));
18 memset(ins,false,sizeof(ins));
19 size=step=scc=0; top=-1;
20 }
21
22 void Insert(int u,int v){
23 E[size].vtx=v;
24 E[size].next=head[u];
25 head[u]=size++;
26 }
27
28 void Tarjan(int u){
29 stk[++top]=u; ins[u]=true;
30 dfn[u]=low[u]=step++;
31 for(int i=head[u];i!=-1;i=E[i].next){
32 int v=E[i].vtx;
33 if(dfn[v]==-1){
34 Tarjan(v);
35 low[u]=MIN(low[u],low[v]);
36 }elseif(ins[v]){
37 low[u]=MIN(low[u],dfn[v]);
38 }
39 }
40 if(low[u]==dfn[u]){
41 for(int v=-1;v!=u;top--){
42 v=stk[top];
43 ins[v]=false;
44 blg[v]=scc;
45 }
46 scc++;
47 }
48 }
49
50 int a[10005],b[10005],c[10005];
51 int main(){
52 int t;
53 scanf("%d",&t);
54 while(t--){
55 int n,m;
56 scanf("%d%d",&n,&m);
57 for(int i=0;i<m;i++){
58 scanf("%d%d%d",&a[i],&b[i],&c[i]);
59 }
60 int low=0,high=m;
61 int mid=(low+high)>>1,ans=mid;
62 while(low<=high){
63 Init();
64 for(int i=0;i<mid;i++){
65 if(c[i]==0){
66 Insert(a[i]<<1,b[i]<<1|1);
67 Insert(b[i]<<1,a[i]<<1|1);
68 }elseif(c[i]==1){
69 Insert(a[i]<<1,b[i]<<1);
70 Insert(a[i]<<1|1,b[i]<<1|1);
71 Insert(b[i]<<1,a[i]<<1);
72 Insert(b[i]<<1|1,a[i]<<1|1);
73 }elseif(c[i]==2){
74 Insert(a[i]<<1|1,b[i]<<1);
75 Insert(b[i]<<1|1,a[i]<<1);
76 }
77 }
78 for(int i=0;i<(n<<1);i++){
79 if(dfn[i]==-1) Tarjan(i);
80 }
81 bool flag=true;
82 for(int i=0;i<n;i++){
83 if(blg[i<<1]==blg[i<<1|1]){
84 flag=false; break;
85 }
86 }
87 if(flag){
88 low=(ans=mid)+1;
89 }else{
90 high=mid-1;
91 }
92 mid=(low+high)>>1;
93 }
94 printf("%d\n",ans);
95 }
96 }
POJ 2749 HDU 1815 Building roads [较难]
较难是难在读题上,题目意思是说每个仓库都必须与两个传送点其中的一个且只有一个建路(仓库与传送点以及传送点与传送点之间的路的距离是曼哈顿距离),在满足一些like与hate的关系的情况下要求任意两个仓库的距离的最大值最小(这里指沿道路的距离)
题目理解的话就很简单了,在二分前先把like和hate的关系全部建好(连在s1(a),连在s2(~a))
like: x->y,~x->~y,y->x,y->~x
hate: x->~y,~x->y,y->~x,~y->x;
然后二分距离(d为s1与s2的距离)
dis_s1[x]+dis_s1[y]>mid 建边 x->~y,y->~x
dis_s2[x]+dis_s2[y]>mid 建边 ~x->y,~y->x
dis_s1[x]+dis_s2[y]+d>mid 建边 x->y,~y->~x
dis_s2[x]+dis_s1[y]+d>mid 建边 ~x->~y,y->x
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<cstdio>
2 #include<cstring>
3 usingnamespace std;
4 #define MIN(x,y) (x>y?y:x)
5 #define N 1005
6
7 struct Edge{
8 int vtx,next;
9 }E[999999],tE[999999];
10 int head[N],size;
11 int thead[N],tsize;
12 int dfn[N],low[N],blg[N],step,scc;
13 int stk[N],top;
14 bool ins[N];
15
16 void Init(){
17 memset(head,-1,sizeof(head));
18 memset(thead,-1,sizeof(thead));
19 memset(dfn,-1,sizeof(dfn));
20 memset(ins,false,sizeof(ins));
21 step=scc=0; top=-1;
22 tsize=size=0;
23 }
24
25 void tInit(){
26 memset(thead,-1,sizeof(thead));
27 memset(dfn,-1,sizeof(dfn));
28 memset(ins,false,sizeof(ins));
29 tsize=step=scc=0; top=-1;
30 }
31
32 void Insert(int u,int v){
33 E[size].vtx=v;
34 E[size].next=head[u];
35 head[u]=size++;
36 }
37 void tInsert(int u,int v){
38 tE[tsize].vtx=v;
39 tE[tsize].next=thead[u];
40 thead[u]=tsize++;
41 }
42
43 void Tarjan(int u){
44 stk[++top]=u; ins[u]=true;
45 dfn[u]=low[u]=step++;
46 for(int i=head[u];i!=-1;i=E[i].next){
47 int v=E[i].vtx;
48 if(dfn[v]==-1){
49 Tarjan(v);
50 low[u]=MIN(low[u],low[v]);
51 }elseif(ins[v]){
52 low[u]=MIN(low[u],dfn[v]);
53 }
54 }
55 for(int i=thead[u];i!=-1;i=tE[i].next){
56 int v=tE[i].vtx;
57 if(dfn[v]==-1){
58 Tarjan(v);
59 low[u]=MIN(low[u],low[v]);
60 }elseif(ins[v]){
61 low[u]=MIN(low[u],dfn[v]);
62 }
63 }
64 if(low[u]==dfn[u]){
65 for(int v=-1;v!=u;top--){
66 v=stk[top];
67 ins[v]=false;
68 blg[v]=scc;
69 }
70 scc++;
71 }
72 }
73
74 int s1x,s1y,s2x,s2y;
75 int dis[505][2];
76 int gdis(int x1,int y1,int x2,int y2){
77 int d1=x1-x2; if(d1<0) d1=-d1;
78 int d2=y1-y2; if(d2<0) d2=-d2;
79 return d1+d2;
80 }
81
82 int main(){
83 int n,a,b;
84 while(~scanf("%d%d%d",&n,&a,&b)){
85 scanf("%d%d%d%d",&s1x,&s1y,&s2x,&s2y);
86 int d=gdis(s1x,s1y,s2x,s2y);
87
88 int low=8000000,high=8000000;
89 for(int i=0;i<n;i++){
90 int x,y;
91 scanf("%d%d",&x,&y);
92 dis[i][0]=gdis(x,y,s1x,s1y);
93 if(low>dis[i][0]) low=dis[i][0];
94 dis[i][1]=gdis(x,y,s2x,s2y);
95 if(low>dis[i][1]) low=dis[i][1];
96 }
97 Init();
98 for(int i=0;i<a;i++){
99 int l,r;
100 scanf("%d%d",&l,&r);
101 l--; r--;
102 Insert(l<<1,r<<1|1);
103 Insert(l<<1|1,r<<1);
104 Insert(r<<1,l<<1|1);
105 Insert(r<<1|1,l<<1);
106 }
107 for(int i=0;i<b;i++){
108 int l,r;
109 scanf("%d%d",&l,&r);
110 l--; r--;
111 Insert(l<<1,r<<1);
112 Insert(l<<1|1,r<<1|1);
113 Insert(r<<1,l<<1);
114 Insert(r<<1|1,l<<1|1);
115 }
116 for(int i=0;i<(n<<1);i++){
117 if(dfn[i]==-1) Tarjan(i);
118 }
119 bool flag=true;
120 for(int i=0;i<n;i++){
121 if(blg[i<<1]==blg[i<<1|1]){flag=false; break;}
122 }
123 if(!flag){
124 puts("-1");
125 continue;
126 }
127 int mid=(low+high)>>1,ans=0;
128 while(low<=high){
129 tInit();
130 for(int i=0;i<n;i++){
131 for(int j=i+1;j<n;j++){
132 if(dis[i][0]+dis[j][0]>mid){
133 tInsert(i<<1,j<<1|1);
134 tInsert(j<<1,i<<1|1);
135 }
136 if(dis[i][0]+dis[j][1]+d>mid){
137 tInsert(i<<1,j<<1);
138 tInsert(j<<1|1,i<<1|1);
139 }
140 if(dis[i][1]+dis[j][1]>mid){
141 tInsert(i<<1|1,j<<1);
142 tInsert(j<<1|1,i<<1);
143 }
144 if(dis[i][1]+dis[j][0]+d>mid){
145 tInsert(i<<1|1,j<<1|1);
146 tInsert(j<<1,i<<1);
147 }
148 }
149 }
150 for(int i=0;i<(n<<1);i++){
151 if(dfn[i]==-1) Tarjan(i);
152 }
153 bool flag=true;
154 for(int i=0;i<n;i++){
155 if(blg[i<<1]==blg[i<<1|1]){flag=false; break;}
156 }
157 if(flag){
158 high=(ans=mid)-1;
159 }else{
160 low=mid+1;
161 }
162 mid=(low+high)>>1;
163 }
164 printf("%d\n",ans);
165 }
166 }
POJ 2723 Get Luffy Out [较难]
两把钥匙为一组,以一个为a,另一个就为~a
对于每一扇门的两把锁X和Y,对应钥匙为x和y,如果我们用了与x一组的钥匙~x,那么x钥匙就不能用了,我们就必须去打开Y锁,使用y钥匙,即边 ~x->y,~y->x
二分能打开门的个数,强连通判解
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<cstdio>
2 #include<cstring>
3 usingnamespace std;
4 #define MIN(x,y) (x>y?y:x)
5 #define N 2050
6
7 struct Edge{
8 int vtx,next;
9 }E[999999];
10 int head[N],size;
11 int dfn[N],low[N],blg[N],step,scc;
12 int stk[N],top;
13 bool ins[N];
14
15 void Init(){
16 memset(head,-1,sizeof(head));
17 memset(dfn,-1,sizeof(dfn));
18 memset(ins,false,sizeof(ins));
19 size=step=scc=0; top=-1;
20 }
21
22 void Insert(int u,int v){
23 E[size].vtx=v;
24 E[size].next=head[u];
25 head[u]=size++;
26 }
27
28 void Tarjan(int u){
29 stk[++top]=u; ins[u]=true;
30 dfn[u]=low[u]=step++;
31 for(int i=head[u];i!=-1;i=E[i].next){
32 int v=E[i].vtx;
33 if(dfn[v]==-1){
34 Tarjan(v);
35 low[u]=MIN(low[u],low[v]);
36 }elseif(ins[v]){
37 low[u]=MIN(low[u],dfn[v]);
38 }
39 }
40 if(low[u]==dfn[u]){
41 for(int v=-1;v!=u;top--){
42 v=stk[top];
43 ins[v]=false;
44 blg[v]=scc;
45 }
46 scc++;
47 }
48 }
49
50 int opp[N];
51 int dl[2050],dr[2050];
52 int main(){
53 int n,m;
54 while(scanf("%d%d",&n,&m),n||m){
55 for(int i=0;i<n;i++){
56 int l,r;
57 scanf("%d%d",&l,&r);
58 opp[l]=r;
59 opp[r]=l;
60 }
61 n<<=1;
62 for(int i=0;i<m;i++){
63 scanf("%d%d",&dl[i],&dr[i]);
64 }
65
66 int low=0,high=m;
67 int mid=(low+high)>>1,ans=mid;
68 while(low<=high){
69 Init();
70 for(int i=0;i<mid;i++){
71 Insert(opp[dl[i]],dr[i]);
72 Insert(opp[dr[i]],dl[i]);
73 }
74 for(int i=0;i<n;i++){
75 if(dfn[i]==-1) Tarjan(i);
76 }
77 bool flag=true;
78 for(int i=0;i<n;i++){
79 if(blg[i]==blg[opp[i]]){
80 flag=false;
81 break;
82 }
83 }
84 if(flag) low=(ans=mid)+1;
85 else high=mid-1;
86 mid=(low+high)>>1;
87 }
88 printf("%d\n",ans);
89
90
91 }
92 }
HDU 1816 Get Luffy Out * [较难]
上题的加强版,这题中一把钥匙可以同时属于多组,我们把每把钥匙使用(a)和不使用(~a)为两个状态,先根据组建边,如果一个组x和y两把钥匙,那么有x->~y,y->~x,然后就和上题一样了
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<cstdio>
2 #include<cstring>
3 usingnamespace std;
4 #define MIN(x,y) (x>y?y:x)
5 #define N 4200
6
7 struct Edge{
8 int vtx,next;
9 }E[999999];
10 int head[N],size;
11 int dfn[N],low[N],blg[N],scc,step;
12 int stk[N],top;
13 bool ins[N];
14
15 void Init(){
16 memset(head,-1,sizeof(head));
17 memset(dfn,-1,sizeof(dfn));
18 memset(ins,false,sizeof(ins));
19 size=step=scc=0; top=-1;
20 }
21
22 void Insert(int u,int v){
23 E[size].vtx=v;
24 E[size].next=head[u];
25 head[u]=size++;
26 }
27
28 void Tarjan(int u){
29 stk[++top]=u; ins[u]=true;
30 dfn[u]=low[u]=step++;
31 for(int i=head[u];i!=-1;i=E[i].next){
32 int v=E[i].vtx;
33 if(dfn[v]==-1){
34 Tarjan(v);
35 low[u]=MIN(low[u],low[v]);
36 }elseif(ins[v]){
37 low[u]=MIN(low[u],dfn[v]);
38 }
39 }
40 if(low[u]==dfn[u]){
41 for(int v=-1;v!=u;top--){
42 v=stk[top];
43 ins[v]=false;
44 blg[v]=scc;
45 }
46 scc++;
47 }
48 }
49
50 int kl[1050],kr[1050];
51 int dl[2100],dr[2100];
52 int main(){
53 int n,m;
54 while(scanf("%d%d",&n,&m),n||m){
55 for(int i=0;i<n;i++){
56 scanf("%d%d",&kl[i],&kr[i]);
57 }
58 for(int i=0;i<m;i++){
59 scanf("%d%d",&dl[i],&dr[i]);
60 }
61
62 int low=0,high=m;
63 int mid=(low+high)>>1,ans=mid;
64 while(low<=high){
65 Init();
66 for(int i=0;i<n;i++){
67 Insert(kl[i]<<1,kr[i]<<1|1);
68 Insert(kr[i]<<1,kl[i]<<1|1);
69 }
70 for(int i=0;i<mid;i++){
71 Insert(dl[i]<<1|1,dr[i]<<1);
72 Insert(dr[i]<<1|1,dl[i]<<1);
73 }
74 for(int i=0;i<(n<<1);i++){
75 if(dfn[i]==-1) Tarjan(i);
76 }
77 bool flag=true;
78 for(int i=0;i<n;i++){
79 if(blg[i<<1]==blg[i<<1|1]){
80 flag=false;break;
81 }
82 }
83 if(flag){
84 low=(ans=mid)+1;
85 }else{
86 high=mid-1;
87 }
88 mid=(low+high)>>1;
89 }
90 printf("%d\n",ans);
91 }
92 }
求方案:
赵爽的论文中提到将相对点染为蓝色后要将其子孙全部染为蓝色,今天看代码时发现这个操作完全是多余的,因为染为红色的点是每次找入度为0的点染色的,那么根据他的对称性,染为蓝色的点的出度也是为0的,那么他的子孙也只能是已经染过蓝色的点了,所以这里直接按照拓扑序贪心就行了(11/8/3更新)
HDU 1814 Peaceful Commission [中等]
因为要求最小字典序解,所以就不能用强连通+拓扑排序的方法做了,直接用暴力的方法求出解即可,复杂度有点高啊
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<cstdio>
2 #include<cstring>
3 #include<vector>
4 usingnamespace std;
5
6 vector<int> gra[16005];
7 constint R=1;
8 constint B=2;
9 constint W=0;
10 int col[16005],n;
11 int cnt,ans[16005];
12 bool dfs(int u) {
13 int size=gra[u].size();
14 if(col[u]==B) returnfalse;
15 if(col[u]==R) returntrue;
16 col[u]=R;
17 col[u^1]=B;
18 ans[cnt++]=u;
19 for(int i=0;i<size;i++){
20 int v=gra[u][i];
21 if(!dfs(v)) returnfalse;
22 }
23 returntrue;
24 }
25
26 bool solve() {
27 memset(col,0,sizeof(col));
28 for(int i=0;i<n;i++){
29 if(col[i]) continue;
30 cnt=0;
31 if(!dfs(i)){
32 for(int j=0;j<cnt;j++){
33 col[ans[j]]=W;
34 col[ans[j]^1]=W;
35 }
36 if(!dfs(i^1)) returnfalse;
37 }
38 }
39 returntrue;
40 }
41
42 int main (){
43 int m;
44 while(~scanf("%d%d",&n,&m)) {
45 n<<=1;
46 for(int i=0;i<n;i++) gra[i].clear();
47 for(int i=0;i<m;i++){
48 int a,b;
49 scanf("%d%d",&a,&b);
50 a--; b--;
51 gra[a].push_back(b^1);
52 gra[b].push_back(a^1);
53 }
54 if(solve()){
55 for(int i=0;i<n;i++){
56 if(col[i]==R) printf("%d\n",i+1);
57 }
58 }else puts("NIE");
59 }
60 }
POJ 3683 Priest John’s Busiest Day [中等]
对时间处理之后就是比较裸的2-sat求可行解了
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<cstdio>
2 #include<cstring>
3 usingnamespace std;
4 #define MIN(x,y) (x>y?y:x)
5 #define N 4005
6
7 struct Edge{
8 int vtx,next;
9 }E[2999999];
10 int head[N],size;
11 int dfn[N],low[N],blg[N],step,scc;
12 int stk[N],top;
13 bool ins[N];
14
15 void Init(){
16 memset(head,-1,sizeof(head));
17 memset(dfn,-1,sizeof(dfn));
18 memset(ins,false,sizeof(ins));
19 step=size=0; top=-1;
20 }
21
22 void Insert(int u,int v){
23 E[size].vtx=v;
24 E[size].next=head[u];
25 head[u]=size++;
26 }
27
28 void Tarjan(int u){
29 stk[++top]=u; ins[u]=true;
30 dfn[u]=low[u]=step++;
31 for(int i=head[u];i!=-1;i=E[i].next){
32 int v=E[i].vtx;
33 if(dfn[v]==-1){
34 Tarjan(v);
35 low[u]=MIN(low[u],low[v]);
36 }elseif(ins[v]){
37 low[u]=MIN(low[u],dfn[v]);
38 }
39 }
40 if(low[u]==dfn[u]){
41 for(int v=-1;v!=u;top--){
42 v=stk[top];
43 ins[v]=false;
44 blg[v]=scc;
45 }
46 scc++;
47 }
48 }
49
50 constint W=0;
51 constint R=1;
52 constint B=2;
53 int deg[N],col[N],opp[N];
54 int que[N],ss,ee;
55
56 void topo(){
57 ss=ee=0;
58 memset(col,W,sizeof(col));
59 for(int i=0;i<scc;i++){
60 if(!deg[i]) que[ee++]=i;
61 }
62 while(ss<ee){
63 int u=que[ss++];
64 if(col[u]) continue;
65 col[u]=R;
66 col[opp[u]]=B;
67 for(int i=head[u];i!=-1;i=E[i].next){
68 int v=E[i].vtx; deg[v]--;
69 if(!deg[v]) que[ee++]=v;
70 }
71 }
72 }
73
74 int l[N],r[N],len[N];
75 int main(){
76 int n;
77 while(~scanf("%d",&n)){
78 for(int i=0;i<n;i++){
79 int lh,lm,rh,rm;
80 scanf("%d:%d %d:%d %d",&lh,&lm,&rh,&rm,&len[i]);
81 l[i]=lh*60+lm;
82 r[i]=rh*60+rm;
83 }
84 Init(); scc=0;
85 for(int i=0;i<n;i++){
86 for(int j=i+1;j<n;j++){
87 if(l[i]<l[j]+len[j]&&l[i]+len[i]>l[j]){
88 Insert(i<<1,j<<1|1);
89 Insert(j<<1,i<<1|1);
90 }
91 if(l[i]<r[j]&&l[i]+len[i]>r[j]-len[j]){
92 Insert(i<<1,j<<1);
93 Insert(j<<1|1,i<<1|1);
94 }
95 if(r[i]-len[i]<l[j]+len[j]&&r[i]>l[j]){
96 Insert(i<<1|1,j<<1|1);
97 Insert(j<<1,i<<1);
98 }
99 if(r[i]-len[i]<r[j]&&r[i]>r[j]-len[j]){
100 Insert(i<<1|1,j<<1);
101 Insert(j<<1|1,i<<1);
102 }
103 }
104 }
105 for(int i=0;i<(n<<1);i++){
106 if(dfn[i]==-1) Tarjan(i);
107 }
108 bool flag=true;
109 for(int i=0;i<n;i++){
110 opp[blg[i<<1]]=blg[i<<1|1];
111 opp[blg[i<<1|1]]=blg[i<<1];
112 if(blg[i<<1]==blg[i<<1|1]){
113 flag=false; break;
114 }
115 }
116 if(flag){
117 Init();
118 memset(deg,0,sizeof(deg));
119 for(int i=0;i<n;i++){
120 for(int j=i+1;j<n;j++){
121 int li=blg[i<<1],ri=blg[i<<1|1];
122 int lj=blg[j<<1],rj=blg[j<<1|1];
123 if(li!=rj&&lj!=ri){
124 if(l[i]<l[j]+len[j]&&l[i]+len[i]>l[j]){
125 Insert(ri,lj); deg[lj]++;
126 Insert(rj,li); deg[li]++;
127 }
128 if(r[i]-len[i]<r[j]&&r[i]>r[j]-len[j]){
129 Insert(lj,ri); deg[ri]++;
130 Insert(li,rj); deg[rj]++;
131 }
132 }
133 if(li!=lj&&ri!=rj){
134 if(l[i]<r[j]&&l[i]+len[i]>r[j]-len[j]){
135 Insert(lj,li); deg[li]++;
136 Insert(ri,rj); deg[rj]++;
137 }
138 if(r[i]-len[i]<l[j]+len[j]&&r[i]>l[j]){
139 Insert(rj,ri); deg[ri]++;
140 Insert(li,lj); deg[lj]++;
141 }
142 }
143 }
144 }
145
146 topo();
147 puts("YES");
148 for(int i=0;i<(n<<1);i++){
149 if(col[blg[i]]==R){
150 if(i&1){
151 int ti=i>>1;
152 int lh=(r[ti]-len[ti])/60,lm=(r[ti]-len[ti])%60;
153 int rh=r[ti]/60,rm=r[ti]%60;
154 printf("%0.2d:%0.2d %0.2d:%0.2d\n",lh,lm,rh,rm);
155 }else{
156 int ti=i>>1;
157 int lh=l[ti]/60,lm=l[ti]%60;
158 int rh=(l[ti]+len[ti])/60,rm=(l[ti]+len[ti])%60;
159 printf("%0.2d:%0.2d %0.2d:%0.2d\n",lh,lm,rh,rm);
160 }
161 }
162 }
163 }else puts("NO");
164 }
165 }
POJ 3648 Wedding [中等]
比较裸的2-sat求可行解的题,题意比较绕,就是说新娘在一边,新娘相对的一边不能存在不和谐的人,求新娘一边的人是那些,这样就要求新郎必须在另一边,所以就要加一条新娘到新郎的边以满足情况,最后输出被染成蓝色的点
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<cstdio>
2 #include<cstring>
3 usingnamespace std;
4 #define MIN(x,y) (x>y?y:x)
5 #define N 2005
6
7 struct Edge{
8 int vtx,next;
9 }E[999999];
10 int head[N],size;
11 int dfn[N],low[N],blg[N],step,scc;
12 int stk[N],top;
13 bool ins[N];
14
15 void Init(){
16 memset(head,-1,sizeof(head));
17 memset(dfn,-1,sizeof(dfn));
18 memset(ins,false,sizeof(ins));
19 step=size=0; top=-1;
20 }
21
22 void Insert(int u,int v){
23 E[size].vtx=v;
24 E[size].next=head[u];
25 head[u]=size++;
26 }
27
28 void Tarjan(int u){
29 stk[++top]=u; ins[u]=true;
30 dfn[u]=low[u]=step++;
31 for(int i=head[u];i!=-1;i=E[i].next){
32 int v=E[i].vtx;
33 if(dfn[v]==-1){
34 Tarjan(v);
35 low[u]=MIN(low[u],low[v]);
36 }elseif(ins[v]){
37 low[u]=MIN(low[u],dfn[v]);
38 }
39 }
40 if(low[u]==dfn[u]){
41 for(int v=-1;v!=u;top--){
42 v=stk[top];
43 ins[v]=false;
44 blg[v]=scc;
45 }
46 scc++;
47 }
48 }
49
50 constint W=0;
51 constint R=1;
52 constint B=2;
53 int deg[N],col[N],opp[N];
54 int que[N],ss,ee;
55
56 void topo(){
57 ss=ee=0;
58 memset(col,W,sizeof(col));
59 for(int i=0;i<scc;i++){
60 if(!deg[i]) que[ee++]=i;
61 }
62 while(ss<ee){
63 int u=que[ss++];
64 if(col[u]) continue;
65 col[u]=R; col[opp[u]]=B;
66 for(int i=head[u];i!=-1;i=E[i].next){
67 int v=E[i].vtx; deg[v]--;
68 if(!deg[v]) que[ee++]=v;
69 }
70 }
71 }
72
73 int l[N],r[N];
74 int main(){
75 int n,m;
76 while(scanf("%d%d",&n,&m),n||m){
77 n<<=1;
78 Init(); scc=0;
79 Insert(1,0);
80 for(int i=0;i<m;i++){
81 char cl,cr;
82 scanf("%d%c%d%c",&l[i],&cl,&r[i],&cr);
83 l[i]<<=1; r[i]<<=1;
84 if(cl=='w') l[i]^=1;
85 if(cr=='w') r[i]^=1;
86 Insert(l[i],r[i]^1);
87 Insert(r[i],l[i]^1);
88 }
89 for(int i=0;i<n;i++){
90 if(dfn[i]==-1) Tarjan(i);
91 }
92 bool flag=true;
93 for(int i=0;i<(n>>1);i++){
94 opp[blg[i<<1]]=blg[i<<1|1];
95 opp[blg[i<<1|1]]=blg[i<<1];
96 if(blg[i<<1]==blg[i<<1|1]){
97 flag=false; break;
98 }
99 }
100 if(flag){
101 Init();
102 memset(deg,0,sizeof(deg));
103 Insert(blg[0],blg[1]);
104 deg[blg[1]]++;
105 for(int i=0;i<m;i++){
106 int rl=blg[l[i]];
107 int rr=blg[r[i]];
108 if(rl==opp[rr]||rr==opp[rl]) continue;
109 Insert(opp[rl],rr);
110 Insert(opp[rr],rl);
111 deg[rr]++;
112 deg[rl]++;
113 }
114
115 topo();
116 for(int i=2;i<n;i++){
117 if(col[blg[i]]==B){
118 if(i>3) putchar('');
119 if(i&1) printf("%dw",i>>1);
120 else printf("%dh",i>>1);
121 }
122 }puts("");
123 }else puts("bad luck");
124 }
125 }
更新:
SRM464 div1 550题
在MatRush大神的博客里看到这题,就去做了下,比较裸的2-sat,bomb game那题的圆这里是方形而已,很快敲完并过了test,正在得意时去看各大神代码才发现自己又水了,因为只有100个点,直接floyd做一遍然后判连通就行了,看来离真正懂2-sat还差的远啊~~
POJ 3905 Perfect Election
北大的第八道2-sat,和前面几题差不多,按照题目要求建图就好了
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<cstdio>
2 #include<cstring>
3 usingnamespace std;
4 #define MIN(x,y) (x>y?y:x)
5 #define N 2005
6 #define M 2000005
7
8 struct Edge{
9 int vtx,next;
10 }E[M];
11 int head[N],size;
12 int dfn[N],low[N],blg[N],step,scc;
13 int stk[N],top;
14 bool ins[N];
15
16 void Init(){
17 memset(head,-1,sizeof(head));
18 memset(dfn,-1,sizeof(dfn));
19 memset(ins,false,sizeof(ins));
20 size=step=scc=0; top=-1;
21 }
22
23 void Insert(int u,int v){
24 E[size].vtx=v;
25 E[size].next=head[u];
26 head[u]=size++;
27 }
28
29 void Tarjan(int u){
30 stk[++top]=u; ins[u]=true;
31 dfn[u]=low[u]=step++;
32 for(int i=head[u];i!=-1;i=E[i].next){
33 int v=E[i].vtx;
34 if(dfn[v]==-1){
35 Tarjan(v);
36 low[u]=MIN(low[u],low[v]);
37 }elseif(ins[v]){
38 low[u]=MIN(low[u],dfn[v]);
39 }
40 }
41 if(low[u]==dfn[u]){
42 for(int v=-1;v!=u;top--){
43 v=stk[top];
44 ins[v]=false;
45 blg[v]=scc;
46 }
47 scc++;
48 }
49 }
50
51 int main(){
52 int n,m;
53 while(~scanf("%d%d",&n,&m)){
54 Init();
55 for(int i=0;i<m;i++){
56 int u,v;
57 scanf("%d%d",&u,&v);
58 if(u>0&&v>0){
59 Insert((u-1)<<1|1,(v-1)<<1);
60 Insert((v-1)<<1|1,(u-1)<<1);
61 }elseif(u>0){
62 Insert((u-1)<<1|1,(-v-1)<<1|1);
63 Insert((-v-1)<<1,(u-1)<<1);
64 }elseif(v>0){
65 Insert((-u-1)<<1,(v-1)<<1);
66 Insert((v-1)<<1|1,(-u-1)<<1|1);
67 }else{
68 Insert((-u-1)<<1,(-v-1)<<1|1);
69 Insert((-v-1)<<1,(-u-1)<<1|1);
70 }
71 }
72 for(int i=0;i<(n<<1);i++){
73 if(dfn[i]==-1) Tarjan(i);
74 }
75 bool flag=true;
76 for(int i=0;i<n;i++){
77 if(blg[i<<1]==blg[i<<1|1]){
78 flag=false;
79 break;
80 }
81 }
82 puts(flag?"1":"0");
83 }
84 }
HDU 4115 Eliminate the Conflict
转载请注明:http://www.cnblogs.com/ambition/archive/2011/07/30/2-sat.html