首页 > 代码库 > poj 2749 Building roads 2-sat
poj 2749 Building roads 2-sat
题意:
给n个村庄的坐标和两个特殊点s1,s2的坐标,现在要将每个村庄连到s1或s2上,使n个村庄间的最大距离最小。
分析:
二分任意两村庄间的最大距离,用2-sat判断该最大距离是否可行,二分的时候可以顺便记录答案,不用等最后区间为空时再输出l或l-1。
代码:
//poj 2749 //sep9 #include <iostream> #include <vector> using namespace std; const int maxN=1024; const int maxL=4000000; vector<int> g[maxN],ng[maxN]; int n,m,cnt,scc,vis[maxN],dfn[maxN]; int x1,y1,x2,y2; int x[maxN],y[maxN]; int a,b; int ap[maxN],aq[maxN],bp[maxN],bq[maxN]; int d[maxN],D; void addegde(int u,int v) { g[u].push_back(v); ng[v].push_back(u); } void dfs(int k) { vis[k]=1; for(int i=g[k].size()-1;i>=0;--i) if(!vis[g[k][i]]) dfs(g[k][i]); dfn[++cnt]=k; } void ndfs(int k) { vis[k]=scc; for(int i=ng[k].size()-1;i>=0;--i) if(!vis[ng[k][i]]) ndfs(ng[k][i]); } void kosaraju() { memset(vis,0,sizeof(vis)); cnt=0; for(int i=1;i<=2*n;++i) if(!vis[i]) dfs(i); memset(vis,0,sizeof(vis)); scc=0; for(int i=2*n;i>=1;--i) if(!vis[dfn[i]]){ ++scc; ndfs(dfn[i]); } } int two_sat(int dis) { int i,j; for(i=1;i<=2*n;++i) g[i].clear(),ng[i].clear(); for(i=1;i<=a;++i){ int p=ap[i],q=aq[i]; addegde(p,q+n); addegde(p+n,q); addegde(q,p+n); addegde(q+n,p); } for(i=1;i<=b;++i){ int p=bp[i],q=bq[i]; addegde(p,q); addegde(q,p); addegde(p+n,q+n); addegde(q+n,p+n); } for(i=1;i<=n;++i) for(j=i+1;j<=n;++j) if(i!=j){ if(d[i]+d[j]>dis){ addegde(i,j+n); addegde(j,i+n); } if(d[i]+d[j+n]+D>dis){ addegde(j+n,i+n); addegde(i,j); } if(d[i+n]+d[j]+D>dis){ addegde(i+n,j+n); addegde(j,i); } if(d[i+n]+d[j+n]>dis){ addegde(i+n,j); addegde(j+n,i); } } kosaraju(); for(i=1;i<=n;++i) if(vis[i]==vis[i+n]) return 0; return 1; } int main() { int i,j; scanf("%d%d%d",&n,&a,&b); scanf("%d%d%d%d",&x1,&y1,&x2,&y2); for(i=1;i<=n;++i) scanf("%d%d",&x[i],&y[i]); for(i=1;i<=n;++i){ d[i]=abs(x[i]-x1)+abs(y[i]-y1); d[i+n]=abs(x[i]-x2)+abs(y[i]-y2); } D=abs(x1-x2)+abs(y1-y2); int p,q; for(i=1;i<=a;++i) scanf("%d%d",&ap[i],&aq[i]); for(i=1;i<=b;++i) scanf("%d%d",&bp[i],&bq[i]); int l=0,r=maxL,ans=-1; while(l<r){ int mid=(l+r)/2; if(two_sat(mid)==0) l=mid+1; else r=mid,ans=mid; } printf("%d",ans); return 0; }
poj 2749 Building roads 2-sat
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。