zoj 3656 2-sat 不错的题
http://acm.zju.edu.cn/onlinejudge/showProblem.do?
problemId=4879
TLE了一下午。然后没办法了 去搜题解 发现思路跟我的差点儿相同 可是就是我的T 后来扩大了数组 然后AC，无语啊
按我的估算 500个点 开到1000+就够了 可是不够 奇怪不懂......
YES or NO的题，一般就是并查集跟2-sat了 目測这道题并查集写起来更easy
http://blog.csdn.net/u011026968/article/details/10823853
看着道题 poj 3678 看完就能发现，仅仅要将每一位都当做poj3678 执行32次 然后就OK
建图方法跟poj3678一样
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <stack>
#include <iostream>
using namespace std; #define CL(a,b) memset(a,b,sizeof(a))
#define IN(s) freopen(s,"r",stdin) const int MAXN = 6000;
int n,b[505][505];
int head[MAXN],dfn[MAXN],low[MAXN],id[MAXN];
int cnt,scnt;
stack<int>st;
struct Node{
int to,nxt;
}edge[1000010];/// inline void addedge(int u,int v,int k)
{
edge[k].to=v;
edge[k].nxt=head[u];
head[u]=k;
//printf("u=%d v=%d k=%d\n",u,v,k);
}
void tarjan(int u)
{
int v,i,min1=dfn[u]=low[u]=cnt++;
st.push(u);
for(i=head[u];i!=-1;i=edge[i].nxt)
{
v=edge[i].to;
if(dfn[v]==-1)tarjan(v);
min1=min(min1,low[v]);
}
if(min1<low[u]){low[u]=min1;return;}
do
{
v=st.top();
id[v]=scnt;
st.pop();
low[v]=n*2;
}while(v!=u);
scnt++;
} int solve(int pos)
{
CL(dfn,0xff);
CL(id,0xff);
CL(head,0xff);
CL(low,0xff); scnt=cnt=0;
while(!st.empty())st.pop();
int num=0;
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
{
if(i==j)continue;///
int f= (b[i][j]>>pos)&1;
//////////
//printf("pos=%d bij=%d f=%d\n",pos,b[i][j],f);
if(i%2 && j%2)
{
if(f)
{
addedge(i,j+n,num++);
addedge(j,i+n,num++); }
else
{ addedge(i+n,i,num++);
addedge(j+n,j,num++);
addedge(i,j,num++);
addedge(j,i,num++);
}
continue;
}
if(i%2==0 && j%2==0)
{
if(f)
{ addedge(i,i+n,num++);
addedge(j,j+n,num++);
addedge(i+n,j+n,num++);
addedge(j+n,i+n,num++);
}
else
{ addedge(j+n,i,num++);
addedge(i+n,j,num++);
}
continue;
}
if(f)
{
addedge(i,j+n,num++);
addedge(j+n,i,num++);
addedge(j,i+n,num++);
addedge(i+n,j,num++);
}
else
{
addedge(i,j,num++);
addedge(j,i,num++);
addedge(i+n,j+n,num++);
addedge(j+n,i+n,num++);
}
}
int flag=1;
for(int i=0;i<n*2;i++)
if(dfn[i] == -1)
{
////////
//printf("tari=%d\n",i);
////////
tarjan(i);
} for(int i=0;i<n;i++)
if(id[i] == id[i+n])//在同一个连通分量
{
flag=0;
break;
}
if(flag)return 1;
else return 0;
} int main()
{
//IN("zoj3656.txt");
while(~scanf("%d",&n))
{
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
scanf("%d",&b[i][j]);
int flag=0;
for(int i=0;i<n;i++)
{
if(b[i][i])
{
flag=2;
puts("NO");
break;
}
}
if(flag==2)continue;
for(int i=0;i<n-1;i++)
for(int j=i+1;j<n;j++)
{
if(b[i][j]!=b[j][i])
{
flag=2;
puts("NO");
break;
}
}
if(flag==2)continue;
for(int i=0;i<32;i++)
{
flag=solve(i);
if(flag==0)break;
}
if(!flag)puts("NO");
else puts("YES");
}
return 0;
}
并查集的做法代码短了非常多 好流弊的样子 http://blog.csdn.net/lasolmi/article/details/38979207
另一种建图方法，网上找的 快了20ms 不是非常理解的说
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <stack>
#include <iostream>
using namespace std; #define CL(a,b) memset(a,b,sizeof(a))
#define IN(s) freopen(s,"r",stdin) const int MAXN = 6000;
int n,b[505][505];
int head[MAXN],dfn[MAXN],low[MAXN],id[MAXN];
int cnt,scnt;
stack<int>st;
struct Node{
int to,nxt;
}edge[1000010];/// inline void addedge(int u,int v,int k)
{
edge[k].to=v;
edge[k].nxt=head[u];
head[u]=k;
//printf("u=%d v=%d k=%d\n",u,v,k);
}
void tarjan(int u)
{
int v,i,min1=dfn[u]=low[u]=cnt++;
st.push(u);
for(i=head[u];i!=-1;i=edge[i].nxt)
{
v=edge[i].to;
if(dfn[v]==-1)tarjan(v);
min1=min(min1,low[v]);
}
if(min1<low[u]){low[u]=min1;return;}
do
{
v=st.top();
id[v]=scnt;
st.pop();
low[v]=n*2;
}while(v!=u);
scnt++;
} int solve(int pos)
{
CL(dfn,0xff);
CL(id,0xff);
CL(head,0xff);
CL(low,0xff); scnt=cnt=0;
while(!st.empty())st.pop();
int num=0;
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
{
if(i==j)continue;///
int f= (b[i][j]>>pos)&1;
//////////
//printf("pos=%d bij=%d f=%d\n",pos,b[i][j],f);
if(i%2 && j%2)
{
if(f)
{
//addedge(i,j+n,num++);
//addedge(j,i+n,num++);
addedge(i+n,j,num++);
addedge(j+n,i,num++);
}
else
{
addedge(i,i+n,num++);
addedge(j,j+n,num++);
//addedge(i+n,i,num++);
//addedge(j+n,j,num++);
//addedge(i,j,num++);
//addedge(j,i,num++);
}
continue;
}
if(i%2==0 && j%2==0)
{
if(f)
{
addedge(i+n,i,num++);
addedge(j+n,j,num++); //addedge(i,i+n,num++);
//addedge(j,j+n,num++);
//addedge(i+n,j+n,num++);
//addedge(j+n,i+n,num++);
}
else
{
addedge(i,j+n,num++);
addedge(j,i+n,num++);
//addedge(j+n,i,num++);
//addedge(i+n,j,num++);
}
continue;
}
if(f)
{
addedge(i,j+n,num++);
addedge(j+n,i,num++);
addedge(j,i+n,num++);
addedge(i+n,j,num++);
}
else
{
addedge(i,j,num++);
addedge(j,i,num++);
addedge(i+n,j+n,num++);
addedge(j+n,i+n,num++);
}
}
int flag=1;
for(int i=0;i<n*2;i++)
if(dfn[i] == -1)
{
////////
//printf("tari=%d\n",i);
////////
tarjan(i);
} for(int i=0;i<n;i++)
if(id[i] == id[i+n])//在同一个连通分量
{
flag=0;
break;
}
if(flag)return 1;
else return 0;
} int main()
{
//IN("zoj3656.txt");
while(~scanf("%d",&n))
{
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
scanf("%d",&b[i][j]);
int flag=0;
for(int i=0;i<n;i++)
{
if(b[i][i])
{
flag=2;
puts("NO");
break;
}
}
if(flag==2)continue;
for(int i=0;i<n-1;i++)
for(int j=i+1;j<n;j++)
{
if(b[i][j]!=b[j][i])
{
flag=2;
puts("NO");
break;
}
}
if(flag==2)continue;
for(int i=0;i<32;i++)
{
flag=solve(i);
if(flag==0)break;
}
if(!flag)puts("NO");
else puts("YES");
}
return 0;
}
zoj 3656 2-sat 不错的题的更多相关文章
- ZOJ Monthly, June 2014 月赛BCDEFGH题题解
比赛链接:点击打开链接 上来先搞了f.c,,然后发现状态不正确,一下午都是脑洞大开,, 无脑wa,无脑ce...一样的错犯2次.. 硬着头皮搞了几发,最后20分钟码了一下G,不知道为什么把1直接当成不 ...
- POJ1062不错的题——spfa倒向建图——枚举等级限制
POJ1062 虽然是中文题目但是还是有一定几率都不准题目意思的:1.所有可能降价的措施不是降价多少钱而是降至多少钱2.等级范围:是你所走的那一条路中所有人中最好最低等级差不允许超过limit限制 思 ...
- ZOJ 3992 One-Dimensional Maze（思维题）
L - One-Dimensional Maze Time Limit:1000MS Memory Limit:65536KB 64bit IO Format:%lld & % ...
- ZOJ 3983 Crusaders Quest（思维题）
C - Crusaders Quest Time Limit:1000MS Memory Limit:65536KB 64bit IO Format:%lld & %llu S ...
- ZOJ 3778 C - Talented Chef 水题
LINK:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3778 题意:有n道菜,每道菜需要\(a_i\)道工序,有m个锅可 ...
- HDU1160:FatMouse's Speed(最长上升子序列，不错的题）
题目:http://acm.hdu.edu.cn/showproblem.php?pid=1160 学的东西还是不深入啊,明明会最长上升子序列,可是还是没有A出这题,反而做的一点思路没有,题意就不多说 ...
- ZOJ - 3890 Wumpus（BFS基础题）
Wumpus Time Limit: 2 Seconds Memory Limit: 65536 KB One day Leon finds a very classic game call ...
- ZOJ 1494 Climbing Worm 数学水题
http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=494 题目大意: 一只蜗牛要从爬上n英寸高的地方,他速度为u每分钟,他爬完u需要 ...
- ZOJ 3652 Maze 模拟,bfs,读题 难度:2
http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=4842 要注意题目中两点: 1.在踏入妖怪控制的区域那一刹那,先减行动力,然后才 ...
随机推荐
- 小技巧-Try Catch
与多线程,业务逻辑等比较复杂的功能打交道时,免不了对部分有可能产生不可预期的代码进行异常捕获. 这种异常可能不处理,比如: try { } catch {} 但是如果一旦发生异常,在程序调试的时候,发 ...
- Highchart 改编风力风向图
基于highchart修改而成的风力风向图 1.替换highchart.js drawpoints方法,替换内容如下: drawPoints: function() { var a, b = this ...
- php错误 mysql_query()：supplied argument is not a valid MySQL result resource
出现这种错误,原因是出现该错误的函数的参数出现了问题 参数出现问题有多种情况: 1.sql查询语句有问题,可能多了一个逗号,少了一个逗号,多了括号之类的: 2.与数据库连接的参数有问题,用户名.密码. ...
- 【原创】JPEG图像密写研究（二） 哈夫曼树的建立
[原创]记录自己研究的过程,仅供参考,欢迎讨论... 在根据JPEG图像文件结构读取完文件后,提取出其中DHT段,利用其中内容建立哈夫曼树,便于之后译码工作.这里需要注意的是文件中的哈夫曼表数量不固定 ...
- 1分钟搞定超慢SQL
前几天,一个用户的研发人员找到我了,说他们有个SQL语句非常慢,我说多慢?他们说:半个小时也没出结果.于是问他们要了SQL语句和执行计划,SQL语句就不能再这里贴出来了,下面是调整前的执行计划(略去某 ...
- 飘逸的python - 用urlparse从url中抽离出想要的信息
最近有个需求,要检测配置中的那些url的域名是否都正常,即是否都能ping通. 不过配置中url格式是这样的 http://www.something.com:1234/ . 要ping的是www.s ...
- Android 开发UI牛博[转]
Android 新兴的UI模式——侧边导航栏 侧边导航栏也就是大家熟知的SliddingMenu,英文也叫Fly-In App Menu.Side Navigation等.当然谷歌现在已经推出类似这个 ...
- U3d 手游优化概述
移动平台瓶颈 体积小 芯片要求改 功耗小 影响计算系能 带宽小 传输方面受限 性能优化 资源方面 美术方面 自带地形(地形是非常占用资源的) a.控制地形的分辨率 b.地形高度图尺寸小于257 c.地 ...
- result 相关
1.dispatcher 2.redirect 3.chain 4.redirectAction 5.freemarker 6.httpheader 7.stream 8.velocity 9.xsl ...
- 一般处理程序在VS2012中打开问题
问题:如果你用vs2012建立的一个一般处理程序,运行查看是,出现这样的界面 原因:VS2012默认使用IIS Web服务器,而不是Visual Studio开发服务器,基于安全考虑IIS默认不允许浏 ...