给出
n
≤
5
e
2
n\leq5e2
n≤5e2个数,然后分别求最长不下降序列长度
s
s
s,选出最多的长度为
s
s
s的不下降子序列的个数,以及
1
1
1和
n
n
n可以选择无限次情形下最多的不下降子序列的个数。
这个题比较先进。对于第一个问题直接dp当前位置开始的最长的长度。
然后对于第二问,我们知道如果dp值为
s
s
s,才有可能是开头,如果dp值为1,则可以成为末尾。建立一个源点向所有的dp值为
s
s
s的点连流量
1
1
1的边,所有dp值为
1
1
1的点向汇点连流量为
1
1
1的边,并且满足
i
<
j
i<j
i<j且
a
i
≤
a
j
a_{i} \leq a_{j}
ai≤aj且
d
p
i
=
d
p
j
+
1
dp_{i}=dp_{j}+1
dpi=dpj+1的点对
i
i
i向
j
j
j连一条流量为
1
1
1的边。每个点只能用一次,所以要拆点,自己向自己连一条流量为
1
1
1的边。
对于第三问,和第二问的区别在于
1
1
1和
n
n
n结点可以使用无限次。判断一下这两个点是否可以作为开头和结尾,如果是的话,源点连无限流量到
1
1
1,
1
1
1连无限流量到自身。结尾同理。
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int N=5e2+7;
int a[N],dp[N];
namespace Dinic {
const int N=1e4+3,M=2e5+3;
int front[N],head[N],dep[N],tot=-1,n,S,T;
struct edge{ int v,nxt;ll w; }e[2*M];
void init() {
tot=-1;
for(int i=1;i<=n;i++) front[i]=-1;
}
void add(int u,int v,ll w) {
e[++tot].v=v,e[tot].w=w,e[tot].nxt=front[u],front[u]=tot;
e[++tot].v=u,e[tot].w=0,e[tot].nxt=front[v],front[v]=tot;
}
bool bfs() {
for(int i=1;i<=n;i++) dep[i]=0;
queue<int> q;
q.push(S);
dep[S]=1;
while(!q.empty()) {
int u=q.front();
q.pop();
for(int i=front[u];i!=-1;i=e[i].nxt) {
int v=e[i].v;
if(dep[v]==0&&e[i].w) {
dep[v]=dep[u]+1;
q.push(v);
}
}
}
return dep[T];
}
ll dfs(int cur,ll dis=1e18) {
if(cur==T) return dis;
ll flow=0,sum=0;
for(int &i=head[cur];i!=-1;i=e[i].nxt) {
int v=e[i].v;
if(dep[v]==dep[cur]+1&&e[i].w) {
flow=dfs(v,min(dis,e[i].w));
dis-=flow,sum+=flow;
e[i].w-=flow,e[i^1].w+=flow;
if(!dis) break;
}
}
return sum;
}
ll solve() {
ll ans=0;
while(bfs()) {
for(int i=1;i<=n;i++) head[i]=front[i];
ans+=dfs(S);
}
return ans;
}
}
int n,mx=0;
void gao1() {
int S=2*n+1;
int T=2*n+2;
Dinic::n=2*n+2;
Dinic::S=S;
Dinic::T=T;
Dinic::init();
for(int i=1;i<=n;i++) {
if(dp[i]==mx) {
Dinic::add(S,i,1);
}
}
for(int i=1;i<=n;i++) {
if(dp[i]==1) {
Dinic::add(i+n,T,1);
}
}
for(int i=1;i<=n;i++)
Dinic::add(i,i+n,1);
for(int i=1;i<=n;i++) {
for(int j=i+1;j<=n;j++) {
if(a[j]>=a[i]&&dp[i]==dp[j]+1) {
Dinic::add(i+n,j,1);
}
}
}
}
void gao2() {
gao1();
int S=2*n+1;
int T=2*n+2;
if(dp[1]==mx) Dinic::add(1,1+n,1e9),Dinic::add(S,1,1e9);
if(dp[n]==1) Dinic::add(n,n+n,1e9),Dinic::add(n+n,T,1e9);
}
int main() {
scanf("%d",&n);
for(int i=1;i<=n;i++)
scanf("%d",&a[i]);
for(int i=n;i>=1;i--) {
dp[i]=1;
for(int j=i+1;j<=n;j++) {
if(a[j]>=a[i]) {
dp[i]=max(dp[i],dp[j]+1);
}
}
mx=max(mx,dp[i]);
}
printf("%d\n",mx);
gao1();
printf("%lld\n",Dinic::solve());
gao2();
printf("%lld\n",Dinic::solve());
return 0;
}