BZOJ.2095.[POI2010]Bridges(最大流ISAP 二分 欧拉回路)
题目链接
最小化最大的一条边,二分答案。然后就变成了给一张无向图定向使其为欧拉回路
二分答案后对于一个位置的两条边可能都保留,即双向边,需要给它定向;可能只保留小的一条,即单向边,不需考虑
如何给它定向呢,或者说怎么形成欧拉回路呢
形成欧拉回路的充要条件:弱连通图;每个点出度=入度
记点i的度数 dgr[i]=out[i]-in[i]
对于一条方向确定的边(u,v),如果改变它的方向,会使dgr[u]-=2,dgr[v]+=2,也就是说u的2度数分配给了v
对于in[x]!=out[x]的点都要通过改变某些边方向满足in[x]==out[x]
最大流,构图就不说了
刚开始时对一条无向边保留一次它反向的机会,随便定个向,就权值小的吧
注意度数是每次u,v都改变1的,建边时dgr/2
如果有某个点出入度之差为奇数,则不可能形成欧拉回路
ISAP初始度数的最大值要是des+1...
Dinic无限TLE原来是M小了...mmp..
ISAP无限WA还是找不出原因...
Dinic:
#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
#define gc() getchar()
const int N=1005,M=2005+2*N,INF=1e9;
int n,m,src,des,dgr[N],u[M],v[M],A[M],B[M];
int Enum,cur[N],H[N],nxt[M<<1],fr[M<<1],to[M<<1],cap[M<<1],lev[N],q[N],path[N],num[N];
inline int read()
{
int now=0,f=1;register char c=gc();
for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
inline void AddEdge(int u,int v,int w)
{
to[++Enum]=v, fr[Enum]=u, nxt[Enum]=H[u], cap[Enum]=w, H[u]=Enum;
to[++Enum]=u, fr[Enum]=v, nxt[Enum]=H[v], cap[Enum]=0, H[v]=Enum;
}
int Rebuild(int x)
{
Enum=1, memset(H,0,sizeof H), memset(dgr,0,sizeof dgr);
for(int i=1;i<=m;++i)
{
if(A[i]<=x) ++dgr[u[i]],--dgr[v[i]];//选择作为单向边或是选中这个方向 需改变的度数就是1!
if(B[i]<=x) AddEdge(u[i],v[i],1);//由u->v可以给出1度数
}
int tot=0;
for(int i=1;i<=n;++i)
if(dgr[i]>0) AddEdge(src,i,dgr[i]>>1),tot+=dgr[i]>>1;
else if(dgr[i]<0) AddEdge(i,des,-dgr[i]>>1);
return tot;//返回不平衡的度数
}
bool BFS2()
{
for(int i=0;i<=des;++i) lev[i]=0,cur[i]=H[i];
lev[src]=1, q[0]=src;
int h=0,t=1;
while(h<t)
{
int x=q[h++];
for(int i=H[x];i;i=nxt[i])
if(!lev[to[i]] && cap[i])
{
lev[to[i]]=lev[x]+1, q[t++]=to[i];
if(to[i]==des) return 1;
}
}
return 0;
}
int Dinic(int u,int flow)
{
if(u==des) return flow;
int used=0;
for(int &i=cur[u];i;i=nxt[i])
if(lev[to[i]]==lev[u]+1 && cap[i])
{
int delta=Dinic(to[i],std::min(cap[i],flow-used));
if(delta)
{
cap[i]-=delta, cap[i^1]+=delta, used+=delta;
if(used==flow) return flow;
}
}
lev[u]=0;
return used;
}
int MF()
{
for(int i=1;i<=n;++i)
if(dgr[i]&1) return -1;//如果有度数为1的点,则不可能形成欧拉回路
int tot=0;
while(BFS2()) tot+=Dinic(src,INF);
return tot;
// return ISAP();
}
int main()
{
n=read(),m=read(),src=0,des=n+1;
int mx=0,mn=INF;
for(int i=1;i<=m;++i)
{
u[i]=read(),v[i]=read(),A[i]=read(),B[i]=read();
if(A[i]>B[i]) std::swap(A[i],B[i]),std::swap(u[i],v[i]);//用A做较小的边,单向边也肯定是较小的边啊
mn=std::min(mn,A[i]), mx=std::max(mx,B[i]);
}
int l=mn,r=mx+1,mid;
while(l<r)
{
mid=l+r>>1;
if(Rebuild(mid)==MF()) r=mid;
else l=mid+1;
}
if(l==mx+1) printf("NIE");
else printf("%d",l);
return 0;
}
ISAP(WA):
#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
#define gc() getchar()
const int N=1005,M=2005+2*N,INF=1e9;
int n,m,src,des,dgr[N],u[M],v[M],A[M],B[M];
int Enum,cur[N],H[N],nxt[M<<1],fr[M<<1],to[M<<1],cap[M<<1],lev[N],q[N],path[N],num[N];
inline int read()
{
int now=0,f=1;register char c=gc();
for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
inline void AddEdge(int u,int v,int w)
{
to[++Enum]=v, fr[Enum]=u, nxt[Enum]=H[u], cap[Enum]=w, H[u]=Enum;
to[++Enum]=u, fr[Enum]=v, nxt[Enum]=H[v], cap[Enum]=0, H[v]=Enum;
}
int Rebuild(int x)
{
Enum=1, memset(H,0,sizeof H), memset(dgr,0,sizeof dgr);
for(int i=1;i<=m;++i)
{
if(A[i]<=x) ++dgr[u[i]],--dgr[v[i]];//选择作为单向边或是选中这个方向 需改变的度数就是1
if(B[i]<=x) AddEdge(u[i],v[i],1);//由u->v可以给出1度数
}
int tot=0;
for(int i=1;i<=n;++i)
if(dgr[i]>0) AddEdge(src,i,dgr[i]>>1),tot+=dgr[i]>>1;
else if(dgr[i]<0) AddEdge(i,des,-dgr[i]>>1);
return tot;//返回不平衡的度数
}
bool BFS()
{
for(int i=src;i<=des;++i) lev[i]=des+1,cur[i]=H[i];
q[0]=des, lev[des]=0;
int h=0,t=1;
while(h<t)
{
int x=q[h++];
for(int i=H[x];i;i=nxt[i])
if(lev[to[i]]==des+1 && cap[i^1])//i^1
lev[to[i]]=lev[x]+1, q[t++]=to[i];
}
return lev[src]<=des;
}
int Augment()
{
int mn=INF;
for(int i=des;i!=src;i=fr[path[i]])
mn=std::min(mn,cap[path[i]]);
for(int i=des;i!=src;i=fr[path[i]])
cap[path[i]]-=mn, cap[path[i]^1]+=mn;
return mn;
}
int ISAP()
{
if(!BFS()) return -1;
memset(num,0,sizeof num);//mmp
for(int i=src;i<=des;++i) ++num[lev[i]];
int x=src,flow=0;
while(lev[src]<=des)
{
if(x==des) flow+=Augment(),x=src;
bool can=0;
for(int i=cur[x];i;i=nxt[i])
if(lev[to[i]]==lev[x]-1 && cap[i])
{
can=1, cur[x]=i, path[x=to[i]]=i;
break;
}
if(!can)
{
int mn=des;
for(int i=H[x];i;i=nxt[i])
if(cap[i]) mn=std::min(mn,lev[to[i]]);
if(!--num[lev[x]]) break;
++num[lev[x]=mn+1], cur[x]=H[x];
if(x!=src) x=fr[path[x]];
}
}
return flow;
}
int MF()
{
for(int i=1;i<=n;++i)
if(dgr[i]&1) return -1;//如果有度数为1的点,则不可能形成欧拉回路
return ISAP();
}
int main()
{
n=read(),m=read(),src=0,des=n+1;
int mx=0,mn=INF;
for(int i=1;i<=m;++i)
{
u[i]=read(),v[i]=read(),A[i]=read(),B[i]=read();
if(A[i]>B[i]) std::swap(A[i],B[i]),std::swap(u[i],v[i]);//用A做较小的边,单向边也肯定是较小的边啊
mn=std::min(mn,A[i]), mx=std::max(mx,B[i]);
}
int l=mn,r=mx+1,mid;
while(l<r)
{
mid=l+r>>1;
if(Rebuild(mid)==MF()) r=mid;
else l=mid+1;
}
if(l==mx+1) printf("NIE");
else printf("%d",l);
return 0;
}
------------------------------------------------------------------------------------------------------------------------
很久以前的奇怪但现在依旧成立的签名
attack is our red sun $$\color{red}{\boxed{\color{red}{attack\ is\ our\ red\ sun}}}$$ ------------------------------------------------------------------------------------------------------------------------
很久以前的奇怪但现在依旧成立的签名
attack is our red sun $$\color{red}{\boxed{\color{red}{attack\ is\ our\ red\ sun}}}$$ ------------------------------------------------------------------------------------------------------------------------