[bzoj1930] [Shoi2003]pacman 吃豆豆

　　一开始想DP做法。。发现不会QAQ

　　暴力费用流挺好想的。。就是拆点限制经过次数，然后每个点的出点往能到达的点的入点连边。。跑个最大费用最大流。

　　但显然边数能达到n^2。。显而易见的优化就是，如果存在路径点a->b->c,那么a就没必要连往c了。（虽然优化完极限情况下边数也差不多是n^2了。。。）

　　但是有可能出现借道的情况。就从入点往出点连一条流量1，费用0的边。（不会被借道两次）

　　有一个神坑。。。这题用spfa跑最长路的时候，队列数组要开到n^2。。。。。因为这个原因卡了一个晚上

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#define ll long long
using namespace std;
const int maxn=4023,inf=1002333333;
struct poi{int x,y;}a[maxn];
struct zs{int pre;short too,flow,dis;}e[4233333];int tot,last[maxn];
short dl[4233333];
int dis[maxn];
bool u[maxn],ins[maxn];
int i,j,k,n,m,s,t,ss;
int ans;

int ra;char rx;
inline int read(){
    rx=getchar(),ra=0;
    while(rx<'0'||rx>'9')rx=getchar();
    while(rx>='0'&&rx<='9')ra*=10,ra+=rx-48,rx=getchar();return ra;
}

bool spfa(){
    memset(dis,255,(t+1)<<2),dis[ss]=0,dl[1]=ss,u[ss]=1;
    int l=0,r=1,i,now;
    while(l<r&&!dis[ss])
        for(i=last[now=dl[++l]],u[now]=0;i;i=e[i].pre)if(e[i].flow&&dis[e[i].too]<dis[now]+e[i].dis){
            dis[e[i].too]=dis[now]+e[i].dis;
            if(!u[e[i].too])u[e[i].too]=1,dl[++r]=e[i].too;
        }
    return dis[t]>-1;
}
inline int min(int a,int b){return a<b?a:b;}
int dfs(int x,int mx){
    if(x==t)return mx;
    int used=0,i,w;ins[x]=1;
    for(i=last[x];i;i=e[i].pre)if(e[i].flow&&dis[e[i].too]==dis[x]+e[i].dis&&!ins[e[i].too])
        if((w=dfs(e[i].too,min(mx-used,(int)e[i].flow)))){
            e[i].flow-=w,e[i^1].flow+=w,ans+=w*e[i].dis,used+=w;
            if(used==mx){ins[x]=0;return mx;}
        }
    dis[x]=-1,ins[x]=0;return used;
}

inline void insert(int a,int b,int c,int d){
    e[++tot].too=b,e[tot].flow=c,e[tot].dis= d,e[tot].pre=last[a],last[a]=tot,
    e[++tot].too=a,e[tot].flow=0,e[tot].dis=-d,e[tot].pre=last[b],last[b]=tot;
}
bool cmp(poi a,poi b){return a.x<b.x||(a.x==b.x&&a.y<b.y);}
int main(){
    n=read();
    for(i=1;i<=n;i++)
        a[i].x=read(),a[i].y=read();
    sort(a+1,a+1+n,cmp);
    ss=0,s=n<<1|1,t=s+1,tot=1;
    insert(ss,s,2,0);
    for(i=1;i<=n;i++){
        insert(i,i+n,1,0),insert(i,i+n,1,1),insert(s,i,1,0),insert(i+n,t,1,0);
        int mn=inf;
        for(j=i+1;j<=n;j++)
            if(a[j].y<mn&&a[j].y>=a[i].y)mn=a[j].y,insert(i+n,j,2,0);
    }
    while(spfa())dfs(ss,inf);
    printf("%d
",ans);
}