多校第16场 HDU 3998 Sequence(最多不相交路径)
2011-09-02 19:51
399 查看
披着DP外衣的网络流
桂哥哥写的没有拆点的算法居然AC了,结果我给他出了一个数据,结果就挂了,说明这题的数据给的太水了。。。
找到最长上升子序列,在所有的最长上升子序列中找到一种元素标号互不相同的情况,使序列个数做多
构图方法:
将所有点拆点,保证其唯一性(在所有序列中只出现一次)
所有的DP为0,即可以做路径起始点的和源点连接,DP为ans的点和汇点相连,边权正无穷。
满足dp[i]+1=dp[j]且a[i]<d[j]的两点连一条正无穷的边,这样保证了每条可行流都是一条长度为ans+1的上升子序列,其他构图方法有的能跑过数据,但是是错的。
桂哥哥写的没有拆点的算法居然AC了,结果我给他出了一个数据,结果就挂了,说明这题的数据给的太水了。。。
找到最长上升子序列,在所有的最长上升子序列中找到一种元素标号互不相同的情况,使序列个数做多
构图方法:
将所有点拆点,保证其唯一性(在所有序列中只出现一次)
所有的DP为0,即可以做路径起始点的和源点连接,DP为ans的点和汇点相连,边权正无穷。
满足dp[i]+1=dp[j]且a[i]<d[j]的两点连一条正无穷的边,这样保证了每条可行流都是一条长度为ans+1的上升子序列,其他构图方法有的能跑过数据,但是是错的。
#include <cstdio> #include <cstring> #define max(a,b) (a>b?a:b) const int maxn=710; const int inf=5; const int s=0; struct edge{ int v,next,w; }edge[40005]; int head[2*maxn],cnt;//for sap void addedge(int u, int v, int w) { edge[cnt].v=v; edge[cnt].w=w; edge[cnt].next=head[u]; head[u]=cnt++; edge[cnt].v=u; edge[cnt].w=0; edge[cnt].next=head[v]; head[v]=cnt++; } int sap(int t) { int pre[2*maxn],cur[2*maxn],dis[2*maxn],gap[2*maxn]; int flow=0 , aug=inf ,u; bool flag; for (int i=0 ; i<=t ; ++i) { cur[i]=head[i]; gap[i]=dis[i]=0; } gap[s]=t+1; u=pre[s]=s; while (dis[s]<=t) { flag=0 ; for (int &j=cur[u] ; ~j ; j=edge[j].next) { int v=edge[j].v; if (edge[j].w>0 && dis[u]==dis[v]+1) { flag=1; if(edge[j].w<aug)aug=edge[j].w; pre[v]=u; u=v; if (u==t) { flow+=aug; while (u!=s) { u=pre[u]; edge[cur[u]].w-=aug; edge[cur[u]^1].w+=aug; } aug=inf; } break; } } if (flag)continue ; int mindis=t+1; for (int j=head[u]; ~j ; j=edge[j].next) { int v=edge[j].v; if (edge[j].w>0 && dis[v]<mindis) { mindis=dis[v]; cur[u]=j; } } if(--gap[dis[u]]==0)break; gap[dis[u]=mindis+1]++; u=pre[u]; } return flow; } void init () { memset (head , -1 , sizeof(head)); cnt=0; } int n; int a[maxn]; int dp[maxn]; void build_graph(int len) { for (int i=0 ; i<n ; ++i) { addedge(i+1 , i+1+n , 1); if(dp[i]==0)//只有DP为0才能作为路径起点,否则会出现更长的序列,与最长子序列矛盾 { addedge(0 , i+1 , inf); } if(dp[i]==len)//只有DP为len才能作为路径终点 { addedge(i+1+n , 1404 , inf); } for (int j=i+1 ; j<n ; ++j) { if(a[i]<a[j]) { addedge(i+1+n , j+1 , inf); } } } } int main () { while (~scanf("%d",&n)) { init(); for (int i=0 ; i<n ; ++i) { scanf("%d",a+i); } memset (dp , 0 , sizeof(dp)); for (int i=0 ; i<n ; ++i) { for (int j=0 ; j<i ; ++j) { if(a[i]>a[j]) dp[i]=max(dp[i],dp[j]+1); } } int ans=0; for (int i=0 ; i<n ; ++i) ans=max(ans,dp[i]); printf("%d\n",ans+1); build_graph(ans); printf("%d\n",sap(1404)); } return 0; }
相关文章推荐
- hdu 3998 Sequence(DP+最大流,求最多的不相交路径)
- hdu 3998 最大流 (最多不相交路径)
- 【网络流】 HDU 3998 Sequence 最多不重合路径
- HDU-3998 Sequence 最多不相交上升子序列
- 【网络流】HDU 3599 War 最多不相交路径
- HDU 3998 Sequence 最长上升子序列+网络流 求不相交的最长上升子序列个数
- hdu 3998 Sequence 求不相交最长上升自序列个数
- hdu 5288 OO’s Sequence(2015多校第一场第1题)枚举因子
- HDU 5288 OO’s Sequence(数学啊 多校2015)
- Wavel Sequence HDU - 6078 多校#4 dp
- [网络流24题] 06 最长递增子序列(最多不相交路径,最大流)
- HDU 3998 Sequence (最长上升子序列 + 最大流)
- HDU 6166 2017多校 Team09 1006:Dijkstra顶点子集最短路径
- hdu 5288 OO’s Sequence(2015多校第一场第1题)枚举因子
- HDU 4915 多校5 Parenthese sequence
- HDU 5828 多校第八场 1008 Rikka with Sequence(线段树--数据加强版)
- HDU 3998 Sequence
- hdu 1025 上面n个点与下面n个点对应连线 求最多能连有多少条不相交的线 (LIS)
- hdu 4606(线段相交 + 最短路 + 最小路径覆盖)
- [多校2015.01.1010 容斥+迭代] hdu 5297 Y sequence