【线性规划与网络流24题 6】最长递增子序列
Time Limit:10000MS Memory Limit:65536K
Total Submit:19 Accepted:2
Case Time Limit:1000MS
Description
给定正整数序列x1, .., xn。
(1)计算其最长递增子序列的长度s。
(2)计算从给定的序列中最多可取出多少个长度为s的递增子序列。
(3)如果允许在取出的序列中多次使用x1和 xn,则从给定序列中最多可取出多少个长度为s的递增子序列。
编程任务:
设计有效算法完成(1)(2)(3)提出的计算任务。
Input
第1 行有 1个正整数n(n<500),表示给定序列的长度。
接下来的1 行有 n个正整数x1, ..., xn.
Output
程序运行结束时,将任务(1)(2)(3)的解答输出。
第1 行是最长递增子序列的长度s。第2行是可取出的长度为s的递增子序列个数。第3行是允许在取出的序列中多次使用x1和xn时可取出的长度为s的递增子序列个数。
Sample Input
4
3 6 2 5
Sample Output
2
2
3
Source
感谢 Wo_ai_WangYuan 放上数据
对于第一个问题,用f[i]表示到第i个数能得到的递增子序列的个数,很容易dp求解,ans=max(ans,f[i])
对于第二个问题,就要用网络流解决了
源点连f[i]==1的i,汇点连f[j]==ans的j,这里注意拆点,两个点之间容量为1,保证对于每一个序列每个数最多用一次,求最大流就是答案
对于第三个问题,源点和1号点的容量和n号点和汇点的容量都变为inf,这样对于不同的序列就可以用多次1号点和n号点了
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<vector>
using namespace std;
const int maxn=1005,inf=1e9;
inline void _read(int &x){
char t=getchar();bool sign=true;
while(t<'0'||t>'9')
{if(t=='-')sign=false;t=getchar();}
for(x=0;t>='0'&&t<='9';t=getchar())x=x*10+t-'0';
if(!sign)x=-x;
}
struct edge{
int a,b,c,flow;
edge(){}
edge(int fr,int t,int c,int f)
:a(fr),b(t),c(c),flow(f){}
};
int s[maxn],f[maxn],ans=-inf,op,ed;
int dis[maxn],vd[maxn],g[maxn][maxn],n;
void dp(){
int i,j;
for(i=1;i<=n;i++){
f[i]=1;
for(j=1;j<i;j++)
if(s[i]>=s[j]&&f[i]<f[j]+1)
f[i]=f[j]+1;
ans=max(ans,f[i]);
}
}
int dfs(int u,int flow){
if(u==ed)return flow;
int v,temp,delta=0;
for(v=1;v<=ed;v++)
if(g[u][v]&&dis[u]==dis[v]+1){
temp=dfs(v,min(flow-delta,g[u][v]));
g[u][v]-=temp;
g[v][u]+=temp;
delta+=temp;
if(delta==flow||dis[op]>=ed)return delta;
}
if(dis[op]>=ed)return delta;
vd[dis[u]]--;
if(vd[dis[u]]==0)dis[op]=ed;
dis[u]++;
vd[dis[u]]++;
return delta;
}
int main(){
_read(n);
int i,j;
for(i=1;i<=n;i++)_read(s[i]);
dp();
cout<<ans<<endl;
if(ans==1){
cout<<n<<endl<<n<<endl;
return 0;
}
op=2*n+1,ed=op+1;
for(i=1;i<=n;i++){
g[i][i+n]=1;
if(f[i]==1)g[op][i]=1;
if(f[i]==ans)g[i+n][ed]=1;
}
for(i=2;i<=n;i++)
for(j=1;j<i;j++)
if(s[i]>=s[j]&&f[i]==f[j]+1)
g[j+n][i]=1;
int maxflow=0;
while(dis[op]<ed)maxflow+=dfs(op,inf);
cout<<maxflow<<endl;
maxflow=0;
memset(dis,0,sizeof(dis));
memset(vd,0,sizeof(vd));
memset(g,0,sizeof(g));
for(i=1;i<=n;i++){
int val=1;
if(i==1||i==n)val=inf;
g[i][i+n]=val;
if(f[i]==1)g[op][i]=val;
if(f[i]==ans)g[i+n][ed]=val;
}
for(i=2;i<=n;i++)
for(j=1;j<i;j++)
if(s[i]>=s[j]&&f[i]==f[j]+1)
g[j+n][i]=1;
while(dis[op]<ed)maxflow+=dfs(op,inf);
cout<<maxflow<<endl;
}