题目
题意:
给定正整数序列 x1,⋯,xn。
1.计算其最长非递减子序列的长度 s。
2.计算从给定的序列中最多可取出多少个长度为 s 的非递减子序列。(给定序列中的每个元素最多只能被取出使用一次)
3.如果允许在取出的序列中多次使用 x1 和 xn,则从给定序列中最多可取出多少个长度为 s 的非递减子序列。
思路: n<=500.第一问可以用n^2的dp,这样便于下两问的求解。
第2问的话,首先特判一下s为1的特殊情况,就是递减子序列,很好判断。我们可以在两个能够dp转移的点之间建一条流量为1的边,源点向所有dp[i] == 1的点连一条流量为1的边,所有dp[i] == s的点向汇点连一条流量为1的边。除此之外记得要拆点,因为每个点只能用一次,形如1->3->4,2->3->5这样的连边也是有可能的,但是3这个点用了两次,就不符合题意了。所以每个点拆成入点和出点,入点和出点之间连一条流量为1的边,可以限制每个数只用一次。
第3问的话,就是把x1和xn的限制从1变成INF,其余不变,基于网络流的特性,在残留网络里跑仍然可以得到结果,不需要重新建图。具体操作就是入点和出点之间的边改成INF,然后源点到x1、xn的边改成INF,x1、xn到汇点的边改成INF.实际上,其中有些情况不会出现,但是都判断也没有问题。
代码:
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 502*4;
const int M = 2*N;
const int INF = 1e9;
int f[502];
int a[502];
int cnt = 0;
int n,m,k,T; int st,ed;
int cur[N],d[N];
int h[N],e[M],ne[M],w[M],idx = 0;
void add(int a,int b,int c)
{
e[idx] = b,w[idx] = c,ne[idx] = h[a],h[a] = idx++;
e[idx] = a,w[idx] = 0,ne[idx] = h[b],h[b] = idx++;
}
int now = 0;
bool bfs()
{
queue<int> q;
memset(d,-1,sizeof(d));
d[st] = 0,cur[st] = h[st],q.push(st);
while(q.size())
{
int u = q.front(); q.pop();
for(int i=h[u];~i;i=ne[i])
{
int j = e[i];
if(d[j]==-1&&w[i])
{
d[j] = d[u]+1;
cur[j] = h[j];
if(j==ed) return 1;
q.push(j);
}
}
}
return 0;
}
ll find(int u,ll limit)
{
if(u==ed) return limit;
ll flow = 0;
for(int i=h[u];~i&&flow<limit;i=ne[i])
{
int j = e[i];
cur[u] = i;
if(d[j]==d[u]+1&&w[i])
{
ll t = find(j,min(1ll*w[i],limit-flow));
if(!t) d[j] = -1;
flow += t,w[i] -= t,w[i^1] += t;
}
}
return flow;
}
ll dinic()
{
ll ans = 0,flow;
while(bfs()) while(flow = find(st,INF)) ans += flow;
return ans;
}
void solve()
{
memset(h,-1,sizeof(h)); idx = 0;
cin>>n;
// cout<<"??\n";
for(int i=1;i<=n;++i)
{
cin>>a[i];
}
int mx = 0;
for(int i=1;i<=n;++i)
{
f[i] = 1;
for(int j=1;j<i;++j)
{
if(a[j]<=a[i]) f[i] = max(f[i],f[j]+1);
}
mx = max(mx,f[i]);
}
cout<<mx<<"\n";
if(mx==1)
{
cout<<n<<"\n"<<n;
return ;
}
ll t1,t2;
st = 0, ed = N-1;
for(int i=1;i<=n;++i)
{
add(i,i+n,1); //拆点,入点向出点连一条流量为1的边.
if(f[i]==1)
{
add(st,i,1);
continue;
}
for(int j=1;j<i;++j)
{
if(a[j]<=a[i]&&f[i]==f[j]+1)
{
add(j+n,i,1); //出点
}
}
if(f[i]==mx)
{
add(n+i,ed,1); //序列结尾,出点
}
}
t1 = dinic();
for(int i=0;i<idx;i+=2)
{
int l = e[i^1],r = e[i];
if((l==1||l==n)&&r==l+n) w[i] = INF;
else if(l==st&&(r==1||r==n)) w[i] = INF;
else if(r==ed&&(l==n+1||l==n+n)) w[i] = INF;
}
t2 = t1 + dinic();
cout<<t1<<"\n"<<t2<<"\n";
}
signed main(void)
{
ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
solve();
return 0;
}