[bzoj]1930 pacman吃豆豆

Description

两个PACMAN吃豆豆。一开始的时候，PACMAN都在坐标原点的左下方，豆豆都在右上方。PACMAN走到豆豆处就会吃掉它。PACMAN行走的路线很奇怪，只能向右走或者向上走，他们行走的路线不可以相交。 请你帮这两个PACMAN计算一下，他们俩加起来最多能吃掉多少豆豆。

Input

第一行为一个整数N，表示豆豆的数目。 接下来 N 行，每行一对正整数，表示第i个豆豆的坐标。任意两个豆豆的坐标都不会重合。

Output

仅有一行包含一个整数，即两个PACMAN加起来最多能吃掉的豆豆数量

Sample Input

8

8 1

1 5

5 7

2 2

7 8

4 6

3 3

6 4

Sample Output

7

HINT

 

N < = 2000

Source

拼尽全力优化，最后还是超时1个点。。。

裸的费用流

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
using namespace std;

const int INF=0x7fffffff;

struct Edge
{
    int from,to,v,c,next;
}E[6000000];
int node=1;
int head[5000],from[5000],dis[5000],vis[5000];

int n,ans,S,T;
struct point
{
    int x,y;
}P[2001];

int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-f;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

void ins(int from,int to,int v,int c)
{
    node++;
    E[node]=(Edge){from,to,v,c,head[from]};
    head[from]=node;
}

void insert(int from,int to,int v,int c)
{
    ins(from,to,v,c);ins(to,from,0,-c);
}

bool spfa()
{
    queue<int> Q;
    for(int i=0;i<=T;i++) dis[i]=-INF;
    Q.push(0);dis[0]=0;vis[0]=1;
    while(!Q.empty())
    {
        int q=Q.front();Q.pop();
        for(int i=head[q];i;i=E[i].next)
            if(E[i].v>0&&dis[q]+E[i].c>dis[E[i].to])
            {
                dis[E[i].to]=dis[q]+E[i].c;
                from[E[i].to]=i;
                if(!vis[E[i].to])
                {
                    Q.push(E[i].to);
                    vis[E[i].to]=1;
                }
            }
        vis[q]=0;
    }
    return dis[T]!=-INF;
}

void mcf()
{
    int x=INF;
    for(int i=from[T];i;i=from[E[i].from])
        x=min(E[i].v,x);
    for(int i=from[T];i;i=from[E[i].from])
    {
        ans+=x*E[i].c;
        E[i].v-=x;E[i^1].v+=x;
    }
}

bool cmp(point a,point b)
{
    return a.x<b.x||(a.x==b.x&&a.y<b.y);
}

int main()
{
    n=read();T=2*n+2;S=T-1;
    for(int i=1;i<=n;i++)
        P[i].x=read(),P[i].y=read();
    sort(P+1,P+n+1,cmp);
    insert(0,S,2,0);
    for(int i=1;i<=n;i++)
    {
        insert(S,i,1,0);
        insert(i,i+n,1,1);
        insert(i,i+n,1,0);
        insert(i+n,T,1,0);
    }
    for(int i=1;i<=n;i++)
    {
        int inf=INF;
        for(int j=i+1;j<=n;j++)
        {
            if(P[i].y<=P[j].y)
            {
                if(P[j].y<inf)
                    insert(i+n,j,2,0);
                inf=min(inf,P[j].y);
            }
        }
    }
    while(spfa())
        mcf();
    printf("%d",ans);
    return 0;
}