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 不错的题的更多相关文章

  1. ZOJ Monthly, June 2014 月赛BCDEFGH题题解

    比赛链接:点击打开链接 上来先搞了f.c,,然后发现状态不正确,一下午都是脑洞大开,, 无脑wa,无脑ce...一样的错犯2次.. 硬着头皮搞了几发,最后20分钟码了一下G,不知道为什么把1直接当成不 ...

  2. POJ1062不错的题——spfa倒向建图——枚举等级限制

    POJ1062 虽然是中文题目但是还是有一定几率都不准题目意思的:1.所有可能降价的措施不是降价多少钱而是降至多少钱2.等级范围:是你所走的那一条路中所有人中最好最低等级差不允许超过limit限制 思 ...

  3. ZOJ 3992 One-Dimensional Maze(思维题)

    L - One-Dimensional Maze Time Limit:1000MS     Memory Limit:65536KB     64bit IO Format:%lld & % ...

  4. ZOJ 3983 Crusaders Quest(思维题)

    C - Crusaders Quest Time Limit:1000MS     Memory Limit:65536KB     64bit IO Format:%lld & %llu S ...

  5. ZOJ 3778 C - Talented Chef 水题

    LINK:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3778 题意:有n道菜,每道菜需要\(a_i\)道工序,有m个锅可 ...

  6. HDU1160:FatMouse's Speed(最长上升子序列,不错的题)

    题目:http://acm.hdu.edu.cn/showproblem.php?pid=1160 学的东西还是不深入啊,明明会最长上升子序列,可是还是没有A出这题,反而做的一点思路没有,题意就不多说 ...

  7. ZOJ - 3890 Wumpus(BFS基础题)

    Wumpus Time Limit: 2 Seconds      Memory Limit: 65536 KB One day Leon finds a very classic game call ...

  8. ZOJ 1494 Climbing Worm 数学水题

    http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=494 题目大意: 一只蜗牛要从爬上n英寸高的地方,他速度为u每分钟,他爬完u需要 ...

  9. ZOJ 3652 Maze 模拟,bfs,读题 难度:2

    http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=4842 要注意题目中两点: 1.在踏入妖怪控制的区域那一刹那,先减行动力,然后才 ...

随机推荐

  1. 小技巧-Try Catch

    与多线程,业务逻辑等比较复杂的功能打交道时,免不了对部分有可能产生不可预期的代码进行异常捕获. 这种异常可能不处理,比如: try { } catch {} 但是如果一旦发生异常,在程序调试的时候,发 ...

  2. Highchart 改编风力风向图

    基于highchart修改而成的风力风向图 1.替换highchart.js drawpoints方法,替换内容如下: drawPoints: function() { var a, b = this ...

  3. php错误 mysql_query():supplied argument is not a valid MySQL result resource

    出现这种错误,原因是出现该错误的函数的参数出现了问题 参数出现问题有多种情况: 1.sql查询语句有问题,可能多了一个逗号,少了一个逗号,多了括号之类的: 2.与数据库连接的参数有问题,用户名.密码. ...

  4. 【原创】JPEG图像密写研究(二) 哈夫曼树的建立

    [原创]记录自己研究的过程,仅供参考,欢迎讨论... 在根据JPEG图像文件结构读取完文件后,提取出其中DHT段,利用其中内容建立哈夫曼树,便于之后译码工作.这里需要注意的是文件中的哈夫曼表数量不固定 ...

  5. 1分钟搞定超慢SQL

    前几天,一个用户的研发人员找到我了,说他们有个SQL语句非常慢,我说多慢?他们说:半个小时也没出结果.于是问他们要了SQL语句和执行计划,SQL语句就不能再这里贴出来了,下面是调整前的执行计划(略去某 ...

  6. 飘逸的python - 用urlparse从url中抽离出想要的信息

    最近有个需求,要检测配置中的那些url的域名是否都正常,即是否都能ping通. 不过配置中url格式是这样的 http://www.something.com:1234/ . 要ping的是www.s ...

  7. Android 开发UI牛博[转]

    Android 新兴的UI模式——侧边导航栏 侧边导航栏也就是大家熟知的SliddingMenu,英文也叫Fly-In App Menu.Side Navigation等.当然谷歌现在已经推出类似这个 ...

  8. U3d 手游优化概述

    移动平台瓶颈 体积小 芯片要求改 功耗小 影响计算系能 带宽小 传输方面受限 性能优化 资源方面 美术方面 自带地形(地形是非常占用资源的) a.控制地形的分辨率 b.地形高度图尺寸小于257 c.地 ...

  9. result 相关

    1.dispatcher 2.redirect 3.chain 4.redirectAction 5.freemarker 6.httpheader 7.stream 8.velocity 9.xsl ...

  10. 一般处理程序在VS2012中打开问题

    问题:如果你用vs2012建立的一个一般处理程序,运行查看是,出现这样的界面 原因:VS2012默认使用IIS Web服务器,而不是Visual Studio开发服务器,基于安全考虑IIS默认不允许浏 ...