24最長のネットワークフロー問題は、サブシーケンスの問題をドロップしません。

トピックポータル

図の構造は、ああ、非常に巧妙な、私は本当に(私もキンマAでした)それを見ていません

まず、あなたが再び実行する必要が最長のサブシーケンス落ちない\(O(N ^ 2) \) に対する解決策、前処理\(I \)最大長減少サブシーケンスの最後の長さ\(DP [I]を\ )、最大検索\(DP [I] \)で示さ、\(CNT \)出力\(CNTを\) 最初の質問は終わりました。

ソース点から長さ\(1 \)容量に接続された点\(1 \)側、長さが\(CNT \)シンク容量に接続された点\(1 \)側。各点は一度だけ使用することができますので、私たちも、それらの間で、二つに容量ポイントスプリットを配置する必要があります\(1 \)側を。後に(\ 2 ^ N)\場合、2点毎に2点を列挙\(I、J \)を満足\(jは{<} I \ ) とDP [I] = DP(\ [J] +1 \) 次いで\(I、J \)もの容量との\(1 \)側。側を実行している最大のフロー図を作成した後、2番目の質問への答えです。

第三の問題は、ソースポイント容易になります(\ 1)\容量の側にもに\(\ RM INF} {\)、\ (1 \)\(1 \)除去後の点\(1「\)もの間で変化容量(\ RM INF} {\)\同様に、\(N- \)\(N '\)との間のエッジ、(\ \ N-')シンクとの間のエッジ(存在する場合)に容量有する(\ \ RMを{} INF \)

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define INF 2147483647
using namespace std;
struct zzz{
    int t,len,nex;
}e[20010<<1]; int head[10010],tot=1;
void add(int x,int y,int z){
    e[++tot].t=y;
    e[tot].len=z;
    e[tot].nex=head[x];
    head[x]=tot;
}
int dis[10010],s,t;
bool bfs(){
    memset(dis,0,sizeof(dis));
    queue <int> q; q.push(s);
    dis[s]=1;
    while(!q.empty()){
        int k=q.front(); q.pop();
        for(int i=head[k];i;i=e[i].nex){
            int to=e[i].t;
            if(dis[to]||!e[i].len) continue;
            dis[to]=dis[k]+1;
            if(to==t) return 1;
            q.push(to);
        }
    }
    return dis[t];
}
int dfs(int now,int flow){
    if(now==t||!flow) return flow;
    int rest=0,fl;
    for(int i=head[now];i;i=e[i].nex){
        int to=e[i].t;
        if(dis[to]==dis[now]+1&&(fl=dfs(to,min(e[i].len,flow)))){
            e[i].len-=fl, e[i^1].len+=fl, rest+=fl;
            if(rest==flow) return rest;
        }
    }
    if(rest<flow) dis[now]=0;
    return rest;
}
int ans;
int dinic(){
    if(bfs()) ans+=dfs(s,INF);
    return ans;
}
int read(){
    int k=0; char c=getchar();
    for(;c<'0'||c>'9';) c=getchar();
    for(;c>='0'&&c<='9';c=getchar())
      k=k*10+c-48;
    return k;
}
int a[510],dp[510],cnt;
int main(){
    int n=read(); t=(n<<1)+1;
    for(int i=1;i<=n;i++) a[i]=read(),dp[i]=1;
    for(int i=1;i<=n;i++)
      for(int j=i+1;j<=n;j++)
        if(a[j]>=a[i]) dp[j]=max(dp[j],dp[i]+1);
    for(int i=1;i<=n;i++) cnt=max(cnt,dp[i]);
    printf("%d\n",cnt);
    if(cnt==1){
        printf("%d\n%d",n,n);
        return 0;
    }
    for(int i=1;i<=n;i++){
        add(i,i+n,1); add(i+n,i,0);
        if(dp[i]==cnt) add(i+n,t,1),add(t,i+n,0);
        if(dp[i]==1) add(s,i,1),add(i,s,0);
    }
    for(int i=1;i<=n;i++)
      for(int j=1;j<i;j++){
          if(a[j]<=a[i]&&dp[j]+1==dp[i])
              add(j+n,i,1),add(i,n+j,0);
      }
    printf("%d\n",dinic());
    add(1,1+n,INF),add(n+1,1,0);
    add(n,n<<1,INF),add(n<<1,n,0);
    if(dp[1]==1) add(s,1,INF), add(1,s,0);
    if(dp[n]==cnt) add(n<<1,t,INF), add(t,n<<1,0);
    printf("%d\n",dinic());
    return 0;
}

おすすめ

転載: www.cnblogs.com/morslin/p/11854756.html