2-SAT问题一般指:P4782 【模板】2-SAT 问题。这样子的题目
解决方法:
1. 先把一个点拆成两个点,分别表示这个点是1或0(一般表示0的点为x,表示1的点为x+n)
2. 根据方程进行推理连边:如通过“ a 为真或 b 为假”可以推至如果 a 是 0 ,那么 b 是 0 ,如果 a 是 1 。那么 b 是 1。通俗理解可以是:
- a 和 b 正在合作写一份代码,a 要求代码1,不包含空格;2,字母全部小写 b 要求代码1,包含空格;2,字母包含大写。那么如果代码打了空格,那么代码里面的字母就全是小写,那如果代码不打空格,那么代码里面的字母包含大写
这样的话,我们就让a与b连边,b+n与a+n连边
3. 判断是否有解,我们只需要跑一遍tarjan寻找强联通分量,如果一个点a和a+n在同一个强联通分量中,那么就吴姐,否则有解。
- 先证在一个强联通分量中吴姐。因为我们的有向边(x为起点,y为终点)表示x可以推导y,如果二者在一个强联通分量,不就表示那个a可以推到a+n,a+n可以推到a么?
- 再证否则有解,因为a和a+n都不在一个强联通分量,那么我们就会发现,一定存在a可以推到a+n而a+n推不到a,或者a+n推到a,而a推不到a+n我们就可以确定a的答案了,因此有解
4. 计算答案,只用枚举每一个数,然后比较强联通分量的编号大小,因为编号更小的会更后被推到(我们tarjan求强联通分量是从叶子节点开始找的)所以只用取a和a+n两个中,编号更小的为答案即可
放上几道典题:
P4782 【模板】2-SAT 问题
#include<bits/stdc++.h>
using namespace std;
const int F=6000000;
int n,m,tot,num,top,cnt;
int head[F],ver[F],nxt[F],dfn[F],low[F],c[F],sack[F],ins[F],out[F];
void add(int x,int y){
ver[++tot]=y;
nxt[tot]=head[x];
head[x]=tot;
}
//x为false;x+n为ture
void fst(){
ios::sync_with_stdio(0);
cin.tie(0); cout.tie(0);
}
void tarjan(int x){
dfn[x]=low[x]=++num;
sack[++top]=x; ins[x]=1;
for(int i=head[x];i;i=nxt[i]){
int y=ver[i];
if(!dfn[y]){
tarjan(y);
low[x]=min(low[x],low[y]);
}
else{
if(ins[y]) low[x]=min(low[x],dfn[y]);
}
}
if(dfn[x]==low[x]){
cnt++;
int y;
do{
y=sack[top--];
ins[y]=0;
c[y]=cnt;
}while(x!=y);
}
}
void badend(){
cout<<"IMPOSSIBLE"<<'\n';
exit(0);
}
void getbian(int x,int i,int y,int j){
if(i==0&&j==0){
add(x+n,y);
add(y+n,x);
}
if(i==0&&j==1){
add(x+n,y+n);
add(y,x);
}
if(i==1&&j==0){
add(x,y);
add(y+n,x+n);
}
if(i==1&&j==1){
add(x,y+n);
add(y,x+n);
}
}
int main(){
//freopen("try.in","r",stdin);
fst();
cin>>n>>m;
while(m--){
int x,y,i,j;
cin>>x>>i>>y>>j;
getbian(x,i,y,j);
}
for(int i=1;i<=2*n;i++){
if(!dfn[i]) tarjan(i);
}
for(int i=1;i<=n;i++){
if(c[i]==c[i+n]) badend();
}
cout<<"POSSIBLE"<<'\n';
for(int i=1;i<=n;i++) cout<<(c[i+n]<c[i])<<" ";
return 0;
}
P4171 [JSOI2010] 满汉全席
把一个点拆成做满菜和做汉菜两种,然后跑2-SAT模版即可
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int F=400;
int t,n,m,tot,cnt,flag,num,top;
int head[F],ver[F*10],nxt[F*10],dfn[F],low[F],c[F],ins[F],sack[F];
string s1,s2;
//x假,x+n真
void add(int x,int y){
ver[++tot]=y;
nxt[tot]=head[x];
head[x]=tot;
}
void init(){
flag=tot=cnt=num=top=0;
for(int i=1;i<=2*n;i++) c[i]=ins[i]=dfn[i]=low[i]=head[i]=0;
}
void fst(){
ios::sync_with_stdio(0);
cin.tie(0); cout.tie(0);
}
void getlen(int x,int i,int y,int j){
if(!i&&!j){
add(x+n,y);
add(y+n,x);
}
if(!i&& j){
add(x+n,y+n);
add(y,x);
}
if( i&&!j){
add(x,y);
add(y+n,x+n);
}
if( i&& j){
add(x,y+n);
add(y,x+n);
}
}
void tarjan(int x){
dfn[x]=low[x]=++num;
sack[++top]=x;
ins[x]=1;
for(int i=head[x];i;i=nxt[i]){
int y=ver[i];
if(!dfn[y]){
tarjan(y);
low[x]=min(low[x],low[y]);
}
else{
if(ins[y]) low[x]=min(low[x],dfn[y]);
}
}
if(dfn[x]==low[x]){
cnt++;
int y;
do{
y=sack[top--];
ins[y]=0;
c[y]=cnt;
}while(x!=y);
}
}
int pow_(int x,int y){
int num=1;
while(y){
if(y%2) num*=x;
x*=x;
y/=2;
}
return num;
}
signed main(){
freopen("try.in","r",stdin);
fst();
cin>>t;
while(t--){
init();
cin>>n>>m;
while(m--){
int x=0,y=0,i,j,len=0,up;
cin>>s1>>s2;
up=s1.length()-1;
i=(s1[0]=='h');
j=(s2[0]=='h');
while(s1[1+len]<='9'&&s1[1+len]>='0'&&1+len<=up) len++;
for(int i=len;i>=1;i--) x+=(s1[i]-'0')*pow_(10,i-len);
up=s2.length()-1;
len=0;
while(s2[1+len]<='9'&&s2[1+len]>='0'&&1+len<=up) len++;
for(int i=len;i>=1;i--) y+=(s2[i]-'0')*pow_(10,i-len);
getlen(x,i,y,j);
}
for(int i=1;i<=n*2;i++){
if(!dfn[i]) tarjan(i);
}
for(int i=1;i<=n;i++){
if(c[i]==c[i+n]&&c[i]&&c[i+n]){
flag=1;
break;
}
}
if(flag) cout<<"BAD"<<endl;
else cout<<"GOOD"<<endl;
}
return 0;
}
UVA11294 Wedding
题解链接:UVA11294 Wedding题解 - 对策 - 洛谷博客
#include<bits/stdc++.h>
using namespace std;
const int F=6000000;
int n,m,tot,cnt,num,top,flag;
int head[F],ver[F],nxt[F],dfn[F],low[F],c[F],ins[F],sack[F];
string s1,s2;
void add(int x,int y){
ver[++tot]=y;
nxt[tot]=head[x];
head[x]=tot;
}
//i表示第i对的丈夫和新娘坐 i+n表示第i对的妻子和新娘坐
void init(){
flag=tot=cnt=num=0;
for(int i=1;i<=4*n;i++) dfn[i]=low[i]=head[i]=c[i]=0;
}
int pow_(int x,int y){
int ans=1;
while(y){
if(y%2) ans*=x;
x*=x;
y/=2;
}
return ans;
}
void tarjan(int x){
low[x]=dfn[x]=++num;
sack[++top]=x; ins[x]=1;
for(int i=head[x];i;i=nxt[i]){
int y=ver[i];
if(!dfn[y]){
tarjan(y);
low[x]=min(low[x],low[y]);
}
else{
if(ins[y]) low[x]=min(low[x],dfn[y]);
}
}
if(dfn[x]==low[x]){
cnt++;
int y;
do{
y=sack[top--];
ins[y]=0;
c[y]=cnt;
}while(x!=y);
}
}
int wk(int x){
if(x>n) return x-n;
return x+n;
}
int main(){
//freopen("try.in","r",stdin);
while(cin>>n>>m){
if(!n&&!m) break;
init();
while(m--){
cin>>s1>>s2;
int len=0,up1=s1.length()-1,up2=s2.length()-1,x=0,y=0;
while(s1[len]<='9'&&s1[len]>='0') len++;
for(int i=len-1;i>=0;i--) x+=(s1[i]-'0')*pow_(10,i+1-len);
len=0;
while(s2[len]<='9'&&s2[len]>='0') len++;
for(int i=len-1;i>=0;i--) y+=(s2[i]-'0')*pow_(10,i+1-len);
x++; y++;
if(s1[up1]=='w'&&s2[up2]=='h'){
add(x,y);
add(y+n,x+n);
}
if(s1[up1]=='h'&&s2[up2]=='w'){
add(x+n,y+n);
add(y,x);
}
if(s1[up1]=='w'&&s2[up2]=='w'){
add(x,y+n);
add(y,x+n);
}
if(s1[up1]=='h'&&s2[up2]=='h'){
add(x+n,y);
add(y+n,x);
}
}
add(1,1+n);
for(int i=1;i<=2*n;i++){
if(!dfn[i]) tarjan(i);
}
for(int i=1;i<=n&&!flag;i++){
if(c[i]==c[i+n]) flag=1;
}
if(flag) cout<<"bad luck"<<'\n';
else{
for(int i=2;i<=n;i++){
if(c[i]>c[i+n]) cout<<i-1<<"w"<<" ";
else cout<<i-1<<"h"<<" ";
}
cout<<'\n';
}
}
return 0;
}
UVA1391 Astronauts
只用注意到每个宇航员只能选择两个任务就很好解决
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int F=6000000;
int n,m,num,tot,jun,top,cnt,flag;
int dfn[F],low[F],head[F],ver[F],nxt[F],c[F],peo[F],sack[F],ins[F];
//i在AB任务中选,i+n选择C任务
void init(){
jun=flag=top=num=tot=cnt=0;
for(int i=1;i<=n*2;i++) head[i]=low[i]=dfn[i]=ins[i]=c[i]=0;
}
void add(int x,int y){
ver[++tot]=y;
nxt[tot]=head[x];
head[x]=tot;
}
void tarjan(int x){
dfn[x]=low[x]=++num;
sack[++top]=x; ins[x]=1;
for(int i=head[x];i;i=nxt[i]){
int y=ver[i];
if(!dfn[y]){
tarjan(y);
low[x]=min(low[x],low[y]);
}
else{
if(ins[y]) low[x]=min(low[x],dfn[y]);
}
}
if(low[x]==dfn[x]){
int y;
cnt++;
do{
y=sack[top--];
ins[y]=0;
c[y]=cnt;
}while(x!=y);
}
}
void fst(){
ios::sync_with_stdio(0);
cin.tie(0); cout.tie(0);
}
signed main(){
fst();
//freopen("try.in","r",stdin);
while(cin>>n>>m){
if(!n&&!m) break;
init();
for(int i=1;i<=n;i++){
cin>>peo[i];
jun+=peo[i];
peo[i]*=n;
}
while(m--){
int x,y;
cin>>x>>y;
if(x!=y){
if((peo[x]>=jun&&peo[y]>=jun)||(peo[x]<jun&&peo[y]<jun)){
add(x+n,y);
add(y+n,x);
add(x,y+n);
add(y,x+n);
}
else{
add(x+n,y);
add(y+n,x);
}
}
}
for(int i=1;i<=2*n;i++){
if(!dfn[i]) tarjan(i);
}
for(int i=1;i<=n;i++){
if(c[i]==c[i+n]){
flag=1;
break;
}
}
if(flag) cout<<"No solution."<<endl;
else{
for(int i=1;i<=n;i++){
if(c[i]>c[i+n]) cout<<"C"<<endl;
else{
if(peo[i]>=jun) cout<<"A"<<endl;
else cout<<"B"<<endl;
}
}
}
}
return 0;
}
CF1215F Radio Stations
#include<bits/stdc++.h>
using namespace std;
const int F=6000000;
struct sl1{
int l,r;
}len[F];
int n,m,p,k,f,up,tot,cnt,num,top,flag;
int head[F],ver[F],nxt[F],dfn[F],low[F],sack[F],ins[F],c[F];
void add(int x,int y){
ver[++tot]=y;
nxt[tot]=head[x];
head[x]=tot;
}
//i不选那个点 i+p选那个点 i+p*2 表示<=i i+p*2+up 表示>i
void tarjan(int x){
dfn[x]=low[x]=++num;
sack[++top]=x; ins[x]=1;
for(int i=head[x];i;i=nxt[i]){
int y=ver[i];
if(!dfn[y]){
tarjan(y);
low[x]=min(low[x],low[y]);
}
else{
if(ins[y]) low[x]=min(low[x],dfn[y]);
}
}
if(low[x]==dfn[x]){
int y;
cnt++;
do{
y=sack[top--];
c[y]=cnt;
ins[y]=0;
}while(x!=y);
}
}
int main(){
//freopen("try.in","r",stdin);
cin>>n>>p>>up>>m;
while(n--){
int x,y;
cin>>x>>y;
add(x+p,y);
add(y+p,x);
}
for(int x,y,i=1;i<=p;i++){
cin>>x>>y;
add(i,2*p+x);
add(i,2*p+y+up+2);
add(2*p+x+up+1,i+p);
add(2*p+y+1,i+p);
}
while(m--){
int x,y;
cin>>x>>y;
add(x,y+p);
add(y,x+p);
}
for(int i=1;i<=up;i++){
add(2*p+i+1,2*p+i);
}
for(int i=up+1;i>=2;i--){
add(2*p+i+up,2*p+i+up+1);
}
add(2*p+1+(up+1),2*p+1);
for(int i=1;i<=2*(p+up+1);i++){
if(!dfn[i]) tarjan(i);
}
for(int i=1;i<=p;i++){
if(c[i]==c[i+p]){
flag=1;
break;
}
if(c[i+p]>c[i]) k++;
}
if(flag) cout<<"-1"<<'\n';
else{
for(int i=1;i<=up;i++){
if(c[2*p+i]<c[2*p+i+up+1]&&c[2*p+i+1]>c[2*p+i+up+2]){
f=i;
break;
}
}
cout<<k<<" "<<f<<'\n';
for(int i=1;i<=p;i++){
if(c[i+p]>c[i]) cout<<i<<" ";
}
}
return 0;
}