LOJ6005 网络流24题

再来一道 短短的

#include
#define inf ~0u>>2
using namespace std;
const int N=1e3+5;
int n,ans,ans2,ans3,s,t,tot=1;
int head[N],dis[N],cur[N],a[N],f[N];
struct Egde{
    int v,nxt,w;
}e[N*2];
queue<int>q;
void add(int u,int v,int w){
    e[++tot].v=v,e[tot].nxt=head[u],e[tot].w=w,head[u]=tot;
    e[++tot].v=u,e[tot].nxt=head[v],e[tot].w=0,head[v]=tot;
}
bool bfs(){
    memset(dis,0,sizeof dis);
    while(!q.empty()) q.pop();
    q.push(s),dis[s]=1;
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=head[x];i;i=e[i].nxt){
            int j=e[i].v;
            if(!dis[j]&&e[i].w){
                dis[j]=dis[x]+1;
                q.push(j); 
            }
        }
    }
    return dis[t]!=0;
}
int dfs(int x,int del){
    if(x==t||!del) return del;
    int ret=0;
    for(int& i=cur[x];i&&del;i=e[i].nxt){
        int j=e[i].v;
        if(e[i].w&&dis[j]==dis[x]+1){
            int dd=dfs(j,min(del,e[i].w));
            e[i].w-=dd;
            e[i^1].w+=dd;
            ret+=dd;
            del-=dd;
        }
    }
    return ret;
}
int dinic(int& o){
    while(bfs()){
        for(int i=1;i<=n+2;i++) cur[i]=head[i];
        o+=dfs(s,inf);
    }
}
int main(){
    cin>>n,s=n+1,t=n+2; 
    for(int i=1;i<=n;i++) {
        cin>>a[i],f[i]=1;
        for(int j=1;jif(a[j]<=a[i])
                f[i]=max(f[i],f[j]+1);
        ans=max(ans,f[i]);
    }
    if(ans==1){
        return printf("1\n%d\n%d",n,n),0;
    }
    for(int i=1;i<=n;i++){
        if(f[i]==1) add(s,i,1);
        if(f[i]==ans) add(i,t,1);
        for(int j=1;jif(a[j]<=a[i]&&f[j]+1==f[i]){
                add(j,i,1);
            }
    } 
    dinic(ans2);
    tot=1;memset(head,0,sizeof head);
    for(int i=1;i<=n;i++){
        if(f[i]==1) add(s,i,inf);
        if(f[i]==ans) add(i,t,inf);
        for(int j=1;jif(a[j]<=a[i]&&f[j]+1==f[i])
                add(j,i,1);
    } 
    dinic(ans3);
    printf("%d\n%d\n%d",ans,ans2,ans3);
}