51nod1331 狭窄的通道

题目传送门

这道题 51nod只Ac了十二个人 没有题解可以研究 所以就自己YY了半天 在这里先感谢一波岚清大爷 orz

然后这道题我分了两种情况 一种是左边的往左跑右边的往右跑 中间有一部分直接走不用跑左右

另一种是全部都要跑左右 然后有一部分点需要左右都跑的 但是左跑和右跑的有一部分点会相互冲突 这个时候需要处理一波矛盾关系

讲的不清楚 给一份数据自己模拟试试吧（再次感谢一波 orzllq）

　　1
　　5 11
　　5 3
　　3 4
　　1 10
　　7 9
　　9 1

这个时候 我们第一个想到的是一波网络流找最小割 经过一波混乱以后 代码成功出世并且意外AC了 贴一波代码

当然也可以贪心枚举而且数据大了会快很多 推荐一篇写的很全的博客吧 orzzsn

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#define LL long long
using namespace std;
const int M=,inf=0x3f3f3f3f;
LL read(){
    LL ans=,f=,c=getchar();
    while(c<''||c>''){if(c=='-') f=-; c=getchar();}
    while(c>=''&&c<=''){ans=ans*+(c-''); c=getchar();}
    return ans*f;
}
LL ans,sum,cnt;
LL n,L,m,S,T;
int f[M][M];
struct node{LL l,r;}e[*M];
bool cmp(node a,node b){return a.l<b.l;}
int first[M],cur[M];
struct note{LL to,next,flow;}q[*M*M];
void ins(LL a,LL b,LL flow){cnt++; q[cnt].to=b; q[cnt].flow=flow; q[cnt].next=first[a]; first[a]=cnt;}
void insert(int a,int b,int flow){ins(a,b,flow); ins(b,a,);}
void find1(LL x,LL k){
    bool f=;
    LL now=e[x].r;
    for(int i=k+;i<=n;i++) if(now>e[i].r) insert(x,i,inf),f=;
    if(f) insert(S,x,*(L-e[x].r));
}
void find2(LL x,LL k){
    LL now=e[x].r;
    for(int i=;i<=k;i++) if(now<e[i].r){
        insert(x,T,*e[x].r);
        break;
    }
}
int check(int x){
    LL mx1=,mx2=L,now=e[x].r;
    for(int i=;i<x;i++) mx1=max(mx1,e[i].r);
    for(int i=x+;i<=n;i++) mx2=min(mx2,e[i].r);
    return mx1<now&&mx2>now;
}
int okay(int x){
    LL mx1=,mn1=L;
    for(int i=;i<=x;i++) mx1=max(mx1,e[i].r);
    for(int i=x+;i<=n;i++) mn1=min(mn1,e[i].r);
    return mx1<mn1;
}
LL calc1(int x){
    LL sum=;
    for(int i=;i<=x;i++) sum=sum+e[i].l+e[i].r;
    return sum;
}
LL calc2(int x){
    LL sum=;
    for(int i=x;i<=n;i++) sum=sum+(*L-e[i].l-e[i].r);
    return sum;
}
LL pd(LL x){return x>=?x:-x;}
LL d[M],st[M];
int bfs(){
    memset(d,-,sizeof(d));
    int head=,tail=;
    st[]=S; d[S]=;
    while(head!=tail){
        int x=st[head++]; if(head>M) head=;
        for(int i=first[x];i;i=q[i].next){
            int now=q[i].to;
            if(q[i].flow&&d[now]==-){st[tail++]=now; d[now]=d[x]+; if(tail>M) tail=;}
        }
    }
    return d[T]!=-;
}
LL dfs(LL x,LL a){
    if(x==T||a==) return a;
    int f,flow=;
    for(int &i=cur[x];i;i=q[i].next){
        int now=q[i].to;
        if(d[now]==d[x]+&&(f=dfs(now,min(a,q[i].flow)))>){
            q[i].flow-=f; q[i^].flow+=f;
            flow+=f;
            a-=f; if(a==) break;
        }
    }
    return flow;
}
LL dinic(){
    LL ans=;
    while(bfs()){for(int i=;i<=T;i++) cur[i]=first[i]; ans+=dfs(S,inf);}
    return ans;
}
int main()
{
    for(int v=read();v;v--){
        ans=inf; n=read(); L=read();
        for(int i=;i<=n;i++) e[i].l=read(),e[i].r=read();
        e[].l=; e[].r=; e[n+].l=L; e[n+].r=L;
        sort(e+,e++n,cmp);
        sum=; for(int i=;i<=n;i++) sum=sum+e[i].l+e[i].r; ans=min(ans,sum);
        sum=; for(int i=;i<=n;i++) sum=sum+(*L-e[i].r-e[i].l); ans=min(ans,sum);
        for(int i=;i<=n;i++)if(okay(i)){
            int k=i+; sum=;
            sum+=calc1(i);
            while(check(k)) sum=sum+pd(e[k].r-e[k].l),k++;
            sum+=calc2(k);
            ans=min(ans,sum);
        }
        for(int i=;i<=n;i++){
            sum=; cnt=; S=; T=n+;
            memset(first,,sizeof(first));
            sum+=calc1(i)+calc2(i+);
            for(int k=;k<=i;k++) find1(k,i);
            for(int k=i+;k<=n;k++) find2(k,i);
            sum+=dinic();
            ans=min(ans,sum);
        }
        printf("%lld\n",ans);
    }
    return ;
}