首页 > 代码库 > BZOJ 1312 Hard Life
BZOJ 1312 Hard Life
就是裸题啦。
注意最后还要跑一遍网络流才是正确的图。
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<queue> #include<cmath> #define maxv 105 #define maxn 1005 #define maxe 1000050 #define inf 1000000007 #define eps 1e-5 using namespace std; int n,m,x[maxn],y[maxn],nume=1,g[maxv],dis[maxv],ans=0,s,t,d[maxv]; bool vis[maxv]; struct edge { int v,nxt; double f; }e[maxe]; queue <int> q; bool comp(double x,double y) { if (fabs(x-y)<=eps) return 0; else if (x<y) return -1; else return 1; } void addedge(int u,int v,double f) { e[++nume].v=v;e[nume].nxt=g[u]; e[nume].f=f;g[u]=nume; e[++nume].v=u;e[nume].nxt=g[v]; e[nume].f=0;g[v]=nume; } void build(double g) { nume=1; for (int i=1;i<=n;i++) addedge(s,i,(double)m),addedge(i,t,(double)m+2*g-d[i]); for (int i=1;i<=m;i++) addedge(x[i],y[i],1),addedge(y[i],x[i],1); } bool bfs() { for (int i=s;i<=t;i++) dis[i]=inf; dis[s]=0;q.push(s); while (!q.empty()) { int head=q.front();q.pop(); for (int i=g[head];i;i=e[i].nxt) { int v=e[i].v; if ((comp(e[i].f,0)>0) && (dis[v]>dis[head]+1)) { dis[v]=dis[head]+1; q.push(v); } } } return dis[t]!=inf; } double dinic(int x,double low) { double ret=0; if (x==t) return low; for (int i=g[x];i && (comp(low,0)>0);i=e[i].nxt) { int v=e[i].v; if ((comp(e[i].f,0)>0) && (dis[v]==dis[x]+1)) { double dd=dinic(v,min(e[i].f,low)); ret+=dd;low-=dd; e[i].f-=dd;e[i^1].f+=dd; } } if (comp(ret,0)<=0) dis[x]=inf; return ret; } double check(double g) { build(g);double max_flow=0; while (bfs()) max_flow+=dinic(s,inf); return (double)m*n-max_flow; } void binary_search() { double l=0,r=m,min_gap=(double)1/n/n,mid; while (r-l>min_gap) { memset(g,0,sizeof(g)); mid=(l+r)/2; double cs=check(mid); if (comp(cs,0)>0) l=mid;else r=mid; } memset(g,0,sizeof(g));double cs=check(l);cs=0; return; } void find(int x) { ans++;vis[x]=true; for (int i=g[x];i;i=e[i].nxt) { if ((comp(e[i].f,0)>0) && (!vis[e[i].v])) find(e[i].v); } } int main() { scanf("%d%d",&n,&m);s=0;t=n+1; if (!m) {printf("1\n");return 0;} for (int i=1;i<=m;i++) { scanf("%d%d",&x[i],&y[i]); d[x[i]]++;d[y[i]]++; } binary_search(); find(s); printf("%d\n",ans-1); return 0; }
BZOJ 1312 Hard Life
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。