随笔 - 531  文章 - 0  评论 - 3  阅读 - 10215 

给定正整数序列 x1…,xn

1.计算其最长不下降子序列的长度 ss。

2.如果每个元素只允许使用一次,计算从给定的序列中最多可取出多少个长度为 ss 的不下降子序列。

3.如果允许在取出的序列中多次使用 x1  xn(其他元素仍然只允许使用一次),

  则从给定序列中最多可取出多少个不同的长度为 ss 的不下降子序列

 

建模:

 求出Fi ,

Fi ==1 Add(S,i,1)

Fi==len Add(i+n,T,1)

这题每个点只能选一次,相当于点权,所以点i拆开, Add(i ,i+n,1)

对于F[j]+1=F[i] , Add(j,i+n,1)

 

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
#include <iostream>
#include<cmath>
#include <queue>
#include <cstring>
using namespace std;
 const int N =2e6+2,M=5e4+100;
  
  
 int g[N],a[N],L;
   
 const int inf =1e9+7;
  
 int all=1,hd[N],go[M],w[M],nxt[M];
   
 int S,T,n,m;
 int dis[M],ans=0,now[M];
   
 void add_(int x,int y,int z){
     nxt[++all]=hd[x]; hd[x]=all; go[all]=y;
     w[all]=z;
     swap(x,y);
     nxt[++all]=hd[x]; hd[x]=all; go[all]=y;
     w[all]=0;
 }
  bool bfs(){
     for(int i=0;i<M;i++)dis[i]=inf;
     queue<int> q;
     q.push(S);
     now[S]=hd[S];
     dis[S]=0;
       
     while(q.empty()==0){
         int x=q.front();
         q.pop();
         for(int i=hd[x];i;i=nxt[i]){
             int y=go[i];
             if(w[i]>0&&dis[y]==inf){
                 dis[y]=dis[x]+1;
                 now[y]=hd[y];
                 q.push(y);
                 if(y==T) return 1;
             }
         }
     }
     return 0;
 }
 int dfs(int x,int sum){
     if(x==T) return sum;
     int k,res=0;
       
     for(int i=now[x];i&& sum ;i=nxt[i]){
         now[x]=i;
         int y=go[i];
         if(w[i]>0&&(dis[y]==dis[x]+1)){
             k=dfs(y,min(sum,w[i]));
             if(k==0) dis[y]=inf;
             w[i]-=k;
             w[i^1]+=k;
             res+=k;
             sum-=k;
         }
     }
     return res;
 }
   
   
 signed main(){
    int i,j,k;
    cin>>n;
    for(i=1;i<=n;i++) cin>>a[i],g[i]=1;
     
    for(i=1;i<=n;i++)
     for(k=1;k<i;k++)
      if(a[i]>=a[k])
        g[i]=max(g[i],1+g[k]);
    for(i=1;i<=n;i++) L=max(L,g[i]);
    if(L==1){
        printf("1\n%d\n%d\n",n,n); return 0 ;
    }
    cout<<L<<endl;
    S=0,T=2*n+1;
    for(i=1;i<=n;++i) if(g[i]==1) add_(S,i,1) ;
    for(i=1;i<=n;i++) if(g[i]==L) add_(i+n,T,1);
    for(i=1;i<=n;++i) add_(i,i+n,1);
     
    for(i=1;i<=n;i++)
     for(k=1;k<i;k++)
      if(g[i]==g[k]+1&&a[i]>=a[k]) add_(k+n,i,1);
    int ans=0;
    while(bfs()) ans+=dfs(S,inf) ;
    cout<<ans<<endl;
     
    add_(1,1+n,inf); add_(S,1,inf);
    if(g[n]==L) add_(n,2*n,inf),add_(2*n,T,inf);
     
    while(bfs()) ans+=dfs(S,inf) ;
    cout<<ans<<endl;
 }

 

 

 

 

posted on   towboat  阅读(29)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 阿里最新开源QwQ-32B,效果媲美deepseek-r1满血版,部署成本又又又降低了!
· 单线程的Redis速度为什么快?
· SQL Server 2025 AI相关能力初探
· AI编程工具终极对决:字节Trae VS Cursor,谁才是开发者新宠?
· 展开说说关于C#中ORM框架的用法!
点击右上角即可分享
微信分享提示