bzoj1930: [Shoi2003]pacman 吃豆豆 费用流

70分的想法很好弄，直接按坐标关系建图，可是边会有n*n条，所以必超时。

满分算法是把每个点先按x排序，然后按类似于凸包的做法从每个点取最优值。

贪心可以证明，假如两个点x所标相同y坐标小的能到达的个数不少于y坐标大的个数；y坐标相同x同理。

所以两个路径应该都是按接近凸包上的点来走，所以路径流量为2，这样就缩小了走的范围。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <iostream>
#include <queue>
using namespace std;
#define MAXN 4100
#define MAXM 4100000
#define INF 0x3f3f3f3f
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
struct node
{
    int v,f,c,next;
}e[MAXM];
int k,head[MAXN],dist[MAXN],from[MAXN],p[MAXN];
bool vis[MAXN];
int en=1,s,_s,t,_t,maxflow,mincost,m,n,num;
struct node2
{
    int x,y;
}a[4100];
void add(int u,int v,int f,int c)
{
    en++;
    e[en].v=v;
    e[en].c=c;
    e[en].f=f;
    e[en].next=head[u];
    head[u]=en;

    en++;
    e[en].v=u;
    e[en].c=-c;
    e[en].f=0;
    e[en].next=head[v];
    head[v]=en;
}
int spfa()
{
        static queue<int> q;
        while(!q.empty())   q.pop();
        memset(dist,0x3f,sizeof(dist));
        int kk=dist[t];
        memset(vis,false,sizeof(vis));
        dist[s] = 0;
        q.push(s);
        while(!q.empty())
        {
            int x = q.front(); q.pop();
            vis[x] = false;
            for(int i = head[x]; i; i = e[i].next)
                if(e[i].f&& dist[e[i].v] > dist[x] + e[i].c)
                {
                    dist[e[i].v] = dist[x] + e[i].c;
                    if(!vis[e[i].v])
                        vis[e[i].v] = true,q.push(e[i].v);
                    from[e[i].v] = x;
                    p[e[i].v] = i;
                }
        }
        return dist[t] != kk;
}
void add2()
{
    int v;
    int maxf=INF;
    for(v=t;v!=s;v=from[v])
        maxf=min(maxf,e[p[v]].f);
    for(v=t;v!=s;v=from[v])
    {
        e[p[v]].f-=maxf;
        e[p[v]^1].f+=maxf;
    }
    maxflow+=maxf;
    mincost+=maxf*dist[t];
}
void init()
{
    maxflow=0;
    s=n*2+1;
    _s=s+1;
    t=_s+1;
    _t=t+1;
    //memset(head,-1,sizeof(head));
}
bool cmp(node2 aa,node2 bb)
{
    if(aa.x!=bb.x)
    return aa.x<bb.x;
    else return aa.y<bb.y;
}
int main()
{
    int i,j;
    n=read();
    init();
    for(int i=1;i<=n;i++)
    {
        a[i].x=read();
        a[i].y=read();
    }
    sort(a+1,a+1+n,cmp);
    for(int i=1;i<=n;i++)
    {
        int minn = INF;
        for(int j=i+1;j<=n;j++)
        {
            if(a[j].y<minn &&a[j].y>=a[i].y)
            {
                add(i+n,j,2,0);
            }
            if(a[j].y >= a[i].y)
                minn = min(minn,a[j].y);
        }
    }
    for(int i=1;i<=n;i++)
    {
        add(_s,i,1,0);
        add(i,i+n,1,-1);
        add(i,i+n,1,0);
        add(i+n,_t,1,0);
    }
    add(s,_s,2,0);
    add(_t,t,2,0);
    int times=0;
    while(spfa())
    {
        add2();
    }
    printf("%d",-1*mincost);
    return 0;
}
/*
8
8 1
1 5
5 7
2 2
7 8
4 6
3 3
6 4
*/