Tarjan/2-SAT

Tags:图论

作业部落

评论地址


Tarjan

用来求割边或者割点,求点双联通分量或者边双联通分量

点双联通分量:两个点之间有两条点不相交的路径

边双联通分量:两个点之间有两条边不相交的路径

Tarjan求LCA还不会

2-SAT

每种物品有选或者不选两种状态,有些限制条件形如

选了\(A\)则必须选\(B\),\(A\)和\(B\)不能同时选,必须选\(A\)等等

把逻辑限制关系变成连边

a->b表示如果\(a\)成立那么\(b\)一定成立

这个要求你理解逆否命题

逆否命题,举个例子,选\(A\)必须选\(B\),那么我选了\(B'\)就不能选\(A\),选\(B'\)就必须选\(A'\)

由于逆否命题产生的对称性使得\(2-SAT\)问题得以在\(O(n)\)时间求解

具体来说要求同时连接x->y y'->x'

这样跑一遍\(Tarjan\)缩点后如果统一物品的两种状态在同一个边双连通分量中就无解

否则可以输出方案,具体来说是每个点选择\(Tarjan\)缩成的超级点中编号最小的那个(也就是反图拓扑序最小的那个)

题单

Tarjan

2-SAT

Code

边双

void Tarjan(int x)
{
vis[x]=1;sta[++top]=x;
dfn[x]=low[x]=++tot;
for(int i=A.head[x],R=A.a[i].to;i;i=A.a[i].next,R=A.a[i].to)
if(!dfn[R]) Tarjan(R),low[x]=min(low[x],low[R]);
else if(vis[R]) low[x]=min(low[x],low[R]);
if(low[x]!=dfn[x]) return;
for(int k=sta[top],lst=0;lst!=x;lst=k,k=sta[--top])
vis[k]=0,bel[k]=x,val[x]+=val[k]*(k!=x);
}

点双

void Tarjan(int x,int f)
{
int s=0;dfn[x]=low[x]=++tot;
for(int i=head[x];i;i=a[i].next)
{
int R=a[i].to;if(R==f) continue;
if(dfn[R]) {low[x]=min(low[x],dfn[R]);continue;}
s++;Tarjan(R,x);tag[x]|=low[R]>=dfn[x];
low[x]=min(low[x],low[R]);
}
if(!f&&s==1) tag[x]=0;
}

!注意点双第7行一定要是dfn[R]

2-sat(和平委员会)

#include<cstdio>
#include<cstdlib>
#include<queue>
using namespace std;
const int N=41000;
struct edge {int next,fr,to;};
struct Map
{
edge a[N]; int head[N],cnt;
void link(int x,int y) {a[++cnt]=(edge){head[x],x,y};head[x]=cnt;}
}A,B;
int dfn[N],top,sta[N],vis[N],low[N];
int bel[N],tot,n,m,p[N],node;
queue<int> Q;
void Min(int &a,int b) {if(b<a) a=b;}
void Tarjan(int x)
{
vis[x]=1;sta[++top]=x; dfn[x]=low[x]=++tot;
for(int i=A.head[x],R=A.a[i].to;i;i=A.a[i].next,R=A.a[i].to)
if(!dfn[R]) Tarjan(R),Min(low[x],low[R]);
else if(vis[R]) Min(low[x],low[R]);
if(low[x]!=dfn[x]) return;++node;
for(int k=sta[top],lst=0;lst!=x;lst=k,k=sta[--top])
vis[k]=0,bel[k]=node;
}
int main()
{
// freopen("spo.in","r",stdin);
// freopen("spo.out","w",stdout);
scanf("%d%d",&n,&m);
for(int i=1;i<=n*2;i++) p[i]=i&1?i+1:i-1;
for(int i=1;i<=m;i++)
{
int x,y;scanf("%d%d",&x,&y);
A.link(x,p[y]);A.link(y,p[x]);
}
for(int i=1;i<=n*2;i++) if(!dfn[i]) Tarjan(i);
for(int i=1;i<=n*2;i+=2) if(bel[i]==bel[p[i]]) {puts("NIE");exit(0);}
for(int i=1;i<=n*2;i+=2) printf("%d\n",bel[i]<bel[p[i]]?i:p[i]);
}

最新文章

  1. Xamarin教程索引页
  2. pathinfo()、dirname()、basename()获得文件的路径,名称等信息说明
  3. iOS 边框显示不完全
  4. 如何解决虚拟机Mac OS X 不支持二进制编译问题()
  5. Eclipse 3.5使用dropins的插件安装方式
  6. codechef Jewels and Stones 题解
  7. shiro错误No SecurityManager accessible to the calling code
  8. java常用集合类:Deque,ArrayList,HashMap,HashSet
  9. ngnix apache tomcat集群负载均衡配置
  10. “Options模式”下各种类型的Options对象是如何绑定的?
  11. leetcode Palindrome Number python
  12. android 学习之RecyclerView
  13. JAVA入门[13]-Spring装配Bean
  14. python的学习和使用
  15. php删除字符串最后一位
  16. Redis自学笔记 --Hash、List、Set类型简述
  17. USB.资料
  18. Android Service 服务(二)—— BroadcastReceiver
  19. ibatis常用的集中判断语句
  20. grep过滤目录或文件方法

热门文章

  1. 3D打印材料的发展现状(1)
  2. 要提高SQL查询效率where语句条件的先后次序应如何写
  3. 显示iOS所有系统字体
  4. RedHat 7 安装PostgreSQL 10.5
  5. &quot;函中函&quot; -------------------- func2(func) -------------- 函数名可以当做函数的参数
  6. 前段js初学总结
  7. Tomcat6的相关配置
  8. Spring各版本源码下载
  9. 取消centOS7虚拟机锁屏
  10. markdownpad 2 pro版本 注册码