「BZOJ2718 / 1143」[Violet 4] 毕业旅行
Description
.gif)
Input
.gif)
Output
最多可选多少景点
Sample Input
7 6
1 2
2 3
5 4
4 3
3 6
6 7
1 2
2 3
5 4
4 3
3 6
6 7
Sample Output
2
HINT
.gif)
题解
最长反链=最小路径覆盖。。。
至于证明。。。百度vfk的博客
floyd传递闭包后，用n-二分图最大匹配数即为答案
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 |
#include<cstdio> #include<cmath> #include<ctime> #include<cstring> #include<iostream> #include<algorithm> #include<queue> #define ll long long #define mod 1000000007 #define inf 1000000000 using namespace std; int read() { int x=0;char ch=getchar(); while(ch<'0'||ch>'9')ch=getchar(); while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} return x; } int n,m,T,cnt=1,ans; int q[405],h[405],last[405],cur[405]; int mp[205][205]; struct edge{ int to,next,v; }e[50005]; void insert(int u,int v,int w) { e[++cnt].to=v;e[cnt].next=last[u];last[u]=cnt;e[cnt].v=w; e[++cnt].to=u;e[cnt].next=last[v];last[v]=cnt;e[cnt].v=0; } void floyd() { for(int k=1;k<=n;k++) for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) mp[i][j]|=(mp[i][k]&mp[k][j]); } bool bfs() { int head=0,tail=1; for(int i=0;i<=T;i++)h[i]=-1; q[0]=0;h[0]=0; while(head!=tail) { int now=q[head];head++; for(int i=last[now];i;i=e[i].next) if(h[e[i].to]==-1&&e[i].v) { h[e[i].to]=h[now]+1; q[tail++]=e[i].to; } } return h[T]!=-1; } int dfs(int x,int f) { if(x==T)return f; int w,used=0; for(int i=cur[x];i;i=e[i].next) if(h[e[i].to]==h[x]+1) { w=dfs(e[i].to,min(e[i].v,f-used)); e[i].v-=w;e[i^1].v+=w; if(e[i].v)cur[x]=i; used+=w;if(used==f)return f; } if(!used)h[x]=-1; return used; } void dinic() { while(bfs()) { for(int i=0;i<=T;i++) cur[i]=last[i]; ans+=dfs(0,inf); } } void build() { for(int i=1;i<=n;i++) { insert(0,i,1); insert(i+n,T,1); } for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) if(mp[i][j]) insert(i,n+j,inf); } int main() { n=read();m=read();T=2*n+1; for(int i=1;i<=m;i++) { int u=read(),v=read(); mp[u][v]=1; } floyd(); build(); dinic(); printf("%d",n-ans); return 0; } |
Subscribe
计数器 9786106
OI 课件=>
OI 课件=>