【Luogu】P2766最长不下降子序列问题(暴力网络流)

  题目链接

  水题qwq,数据都那么水。

  我要是出数据的人我就卡$n^3$建图。

  qwq。

  然而这么水的题我!居!然!没!有!1!A!!还!提!交!了!五!遍!!!

  md从现在开始要锻炼1A率了

  看我从今往后做完一道题之后至少检查TM十分钟

  可恶qwq。

  第一问$n^2$sbDP可解。然而你们知道我提交五遍TM是错在哪里了吗?????

  我TM就错在这个pj-,sb到不能再sb的sb暴力DP上!!!

  气死我了!!!

  关于第二问和第三问,先拆点再拆点qwq。

  先把每个点拆成入点和出点用来限制流量,然后把出点拆成s个点为的是跑最大流的时候统计比较方便。

  能跑到汇点就是一种方式,跑不到就说明没有那么长的最长不下降子序列。

  第三问暴力修改边容量再跑一遍即可。

  qwqqqqqqqqq我网络流都没写错栽在一个入门DP上

  qwqqqqqq

  

#include<cstdio>
#include<cctype>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<queue>
#define maxm 200010
using namespace std;
inline long long read(){
    long long num=0,f=1;
    char ch=getchar();
    while(!isdigit(ch)){
        if(ch=='-')    f=-1;
        ch=getchar();
    }
    while(isdigit(ch)){
        num=num*10+ch-'0';
        ch=getchar();
    }
    return num*f;
}

int n,m;
int q[100010];
int s[100010];
int ne[510][510],tot[510];
int Start,End;
inline int count(int i){    return i&1?i+1:i-1;    }

struct Edge{
    int next,to,val;
}edge[maxm*4];
int head[maxm],num;
inline void addedge(int from,int to,int val){
    edge[++num]=(Edge){head[from],to,val};
    head[from]=num;
}
inline void add(int from,int to,int val){
    addedge(from,to,val);
    addedge(to,from,0);
}

bool vis[maxm];
int dfn[maxm];
int list[maxm];

bool bfs(){
    memset(vis,0,sizeof(vis));    vis[Start]=1;    dfn[Start]=1;
    queue<int>q;    q.push(Start);
    while(!q.empty()){
        int from=q.front();q.pop();
        for(int i=head[from];i;i=edge[i].next){
            int to=edge[i].to;
            if(edge[i].val==0||vis[to])    continue;
            vis[to]=1;    dfn[to]=dfn[from]+1;
            q.push(to);
        }
    }
    return vis[End];
}

int dfs(int x,int val){
    if(x==End||val==0)    return val;
    int flow=0;    vis[x]=1;
    for(int &i=list[x];i;i=edge[i].next){
        int to=edge[i].to;
        if(vis[to]||edge[i].val==0||dfn[to]!=dfn[x]+1)    continue;
        int now=dfs(to,min(val,edge[i].val));
        val-=now;    flow+=now;    edge[i].val-=now;    edge[count(i)].val+=now;
        if(val<=0)    break;
    }
    if(val!=flow)    dfn[x]=-1;
    return flow;
}

int maxflow(){
    int ans=0;
    while(bfs()){
        memset(vis,0,sizeof(vis));
        for(int i=Start;i<=End;++i)    list[i]=head[i];
        int now=dfs(Start,0x7fffffff);
        if(!now)    break;
        ans+=now;
    }
    return ans;
}

int main(){
    int n=read();
    for(int i=1;i<=n;++i){
        q[i]=read();s[i]=1;
    }
    m=1;
    for(int i=n-1;i>=1;--i){
        for(int j=i+1;j<=n;++j){
            if(q[i]>q[j])    continue;
            if(s[i]<s[j]+1)    s[i]=s[j]+1;
            if(m<s[i])    m=s[i];
        }
    }
    printf("%d
",m);
    End=(m+1)*n+n+1;
    for(int i=1;i<=n;++i){
        add(Start,i,0x7fffffff);
        add(i,i+n,1);
        add(m*n+i,End,1);
        for(int j=i+1;j<=n;++j){
            if(q[i]>q[j])    continue;
            for(register int k=1;k<m;++k)    add(k*n+i,(k+1)*n+j,1);
        }
    }
    int ans=maxflow();
    printf("%d
",ans);
    int ret=1;
    for(int j=head[ret];j;j=edge[j].next){
        int to=edge[j].to;
        if(to<ret)    continue;
        edge[j].val+=1234567;
    }
    /*for(int j=head[1];j;j=edge[j].next){
        int to=edge[j].to;
        edge[j].val+=1234567;
    }*/
    ret=m*n+n;
    for(int j=head[ret];j;j=edge[j].next){
        int to=edge[j].to;
        if(to!=End)    continue;
        edge[j].val+=1234567;
    }
    printf("%d",ans+maxflow());
    return 0;
}
原文地址:https://www.cnblogs.com/cellular-automaton/p/8270123.html