博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
pku 3207 Ikki's Story IV - Panda's Trick 2-sat判定是否存在可行解
阅读量:5009 次
发布时间:2019-06-12

本文共 2750 字,大约阅读时间需要 9 分钟。

题意:

一个圆盘的边沿上有n个点, 下标从0开始, 有m条线连接2m个互不相同的点, 线可以在圆盘内部,也可以在圆盘外部, 要求任意两条线不能相交. 给出m条线(内外随意), 问是否满足每条线都不相交.

思路:

可以将第i条线看成一对顶点,编号分别为2*i和2*i+1.那么如果线段i与j相交,就在2*i与2*j+1以及2*i+1与2*j之间连一条双向边。然后就转化到2-sat上判断是否存在可行解了。

//#pragma comment(linker,"/STACK:327680000,327680000")#include 
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define CL(arr, val) memset(arr, val, sizeof(arr))#define ll long long#define inf 0x7f7f7f7f#define lc l,m,rt<<1#define rc m + 1,r,rt<<1|1#define pi acos(-1.0)#define ll long long#define L(x) (x) << 1#define R(x) (x) << 1 | 1#define MID(l, r) (l + r) >> 1#define Min(x, y) (x) < (y) ? (x) : (y)#define Max(x, y) (x) < (y) ? (y) : (x)#define E(x) (1 << (x))#define iabs(x) (x) < 0 ? -(x) : (x)#define OUT(x) printf("%I64d\n", x)#define lowbit(x) (x)&(-x)#define Read() freopen("din.txt", "r", stdin)#define Write() freopen("dout.txt", "w", stdout);#define N 1007#define M 1007using namespace std;struct side{ int s,e;}sid[M];int n,m;struct node{ int v; int next;}g[M*M];int head[M],ct;int dfn[M],low[M];int belong[M],stk[M];bool isn[M];int idx,cnt,top;bool isok(int i,int j){ //判断相交,建图是关键 if ((sid[i].s < sid[j].s && sid[i].e > sid[j].s && sid[i].e < sid[j].e) || (sid[i].s > sid[j].s && sid[i].s < sid[j].e && sid[i].e > sid[j].e)) return true; else return false;}void add(int u,int v){ g[ct].v = v; g[ct].next = head[u]; head[u] = ct++; g[ct].v = u; g[ct].next = head[v]; head[v] = ct++;}void tarjan(int u){ int i,j; dfn[u] = low[u] = ++idx; stk[++top] = u; isn[u] = true; for (i = head[u]; i != - 1; i = g[i].next) { int v = g[i].v; if (dfn[v] == -1) { tarjan(v); low[u] = min(low[u],low[v]); } else if (isn[v]) { low[u] = min(low[u],dfn[v]); } } if (dfn[u] == low[u]) { cnt++; do { j = stk[top--]; belong[j] = cnt; isn[j] = false; }while (j != u); }}void solve(){ int i; for (i = 0; i < 2*m; ++i) { dfn[i] = low[i] = -1; isn[i] = false; belong[i] = 0; } idx = cnt = top = 0; for (i = 0; i < 2*m; ++i) { if (dfn[i] == -1) tarjan(i); } bool flag = false; for (i = 0; i < m; ++i) { if (belong[2*i] == belong[2*i + 1]) { flag = true; break; } } if (flag) printf("the evil panda is lying again\n"); else printf("panda is telling the truth...\n");}int main(){ // Read(); int i,j; while (~scanf("%d%d",&n,&m)) { for (i = 0; i < m; ++i) { scanf("%d%d",&sid[i].s,&sid[i].e); if (sid[i].s > sid[i].e) swap(sid[i].s,sid[i].e); } CL(head,-1); ct = 0; for (i = 0; i < m; ++i) { for (j = i + 1; j < m; ++j) { if (isok(i,j)) { add(2*i,2*j + 1); add(2*j,2*i + 1); } } } solve(); } return 0;}

  

 

转载于:https://www.cnblogs.com/E-star/archive/2013/01/29/2882041.html

你可能感兴趣的文章
Elastic Search 语法总结
查看>>
py自动化之环境配置
查看>>
Winodws SNMP服务安装和配置(Windows 2003 & 2008 R2)
查看>>
红黑树-想说爱你不容易
查看>>
【题目】英文字符进行频率的统计,直方图输出
查看>>
LeetCode-Binary Tree Level Order Traversal
查看>>
COM组件开发实践
查看>>
yii2 源码分析1从入口开始
查看>>
浅谈网站推广
查看>>
Away3D基础之摄像机
查看>>
Leetcode 128. Longest Consecutive Sequence
查看>>
程序员必须知道的几个Git代码托管平台
查看>>
导电塑料入梦来
查看>>
C# 线程手册 第五章 扩展多线程应用程序 - 什么是线程池
查看>>
笔记1126ASP.NET面试题(转)
查看>>
考研路茫茫--单词情结 - HDU 2243(AC自动机+矩阵乘法)
查看>>
HTTP运行期与页面执行模型
查看>>
tableView优化方案
查看>>
近期思考(2019.07.20)
查看>>
Apache2.4使用require指令进行访问控制
查看>>