bzoj1930: [Shoi2003]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

题解

费用流大水题。。
构图：

    init(st,S,2,0);
    init(T,ed,2,0);
    init(S,u,1,0);
    if(i可以连到u)
    {
    init(u+n,T,1,0);
    init(u,u+n,1,1);
    init(i+n,u,MAX,0);
    }

但是要是你直接连边向可以到的点的话会MLE。。
那怎么办呢，我们发现其实有一条路i—-j——k
那么i到k是不用连边的
于是边数就少了很多。。
但其实你仔细思考一下，发现仅仅这样答案会错~~其实连样例也过不了~
原因就是“中转站”只能用一次
那怎么办呢，我们可以再连一条边

init(i,u,1,0);

表示不获取代价，但是可以用他连出去的边就好了

当然还有另外一种构图方法是长这样子的:
我用这个也A了一次了

for (int u=1;u<=n;u++)
    {
        init(S,u,1,0);
        init(u+n,T,1,0);
        init(u,u+n,1,1);
        init(u,u+n,MAX,0);
        int maxx=0;
        for (int i=u+1;i<=n;i++)
        {
            if (a[i].y<=maxx) continue;
            if (a[i].y<=a[u].y)
            {
                maxx=a[i].y;
                init(i+n,u,MAX,0);
            }
        }
    }

CODE：（空间随便开的，不要来hack我。。。我也没算，反正A了）

#include<cstdio>
#include<cstdlib>
#include<queue>
#include<algorithm>
#include<cstring>
#include<iostream>
#include<cstring>
using namespace std;
const int N=2000*2+6;
const int MAX=1<<30;
int n;
struct qq
{
    int x,y;
}a[N];
int st,ed;
struct qt
{
    int x,y,z,z1,last;//流量    费用 
}s[2000*500];int num,last[N];
void init (int x,int y,int z,int z1)
{
    num++;
    s[num].x=x;s[num].y=y;s[num].z=z;s[num].z1=z1;
    s[num].last=last[x];
    last[x]=num;
    swap(x,y);z=0;z1=-z1;
    num++;
    s[num].x=x;s[num].y=y;s[num].z=z;s[num].z1=z1;
    s[num].last=last[x];
    last[x]=num;
}

int from[N];
int f[N];
bool in[N];
bool SPFA ()
{
    memset(in,false,sizeof(in));
    memset(from,-1,sizeof(from));
    memset(f,-1,sizeof(f));
    queue<int> q;
    q.push(st);
    f[st]=0;in[st]=true;
    while (!q.empty())
    {
        int x=q.front();q.pop();
        for (int u=last[x];u!=-1;u=s[u].last)
        {
            int y=s[u].y;
            if (s[u].z>0&&f[y]<f[x]+s[u].z1)
            {
                f[y]=f[x]+s[u].z1;
                from[y]=u;
                if (in[y]==false)
                {
                    in[y]=true;
                    q.push(y);
                }
            }
        }
        in[x]=false;
    }
    return from[ed]!=-1;
}
int mymin (int x,int y)
{
    return x<y?x:y;
}
int ans=0;
void lalal ()
{
    int x=from[ed],t=MAX;
    while (x!=-1)
    {
        t=mymin(t,s[x].z);
        x=from[s[x].x];
    }

    x=from[ed];
    while (x!=-1)
    {
        ans=ans+t*s[x].z1;
        s[x].z-=t;
        s[x^1].z+=t;
        x=from[s[x].x];
    }
    return ;
}
bool cmp (qq A,qq B){return A.x==B.x?A.y>B.y:A.x>B.x;}
int main()
{
    num=1;memset(last,-1,sizeof(last));
    scanf("%d",&n);st=2*n+1;ed=st+1;int S=ed+1,T=S+1;
    init(st,S,2,0);
    init(T,ed,2,0);
    for (int u=1;u<=n;u++)  scanf("%d%d",&a[u].x,&a[u].y);
    sort(a+1,a+1+n,cmp);
    for (int u=1;u<=n;u++)
    {
        init(S,u,1,0);
        init(u+n,T,1,0);
        init(u,u+n,1,1);
        int maxx=0;
        for (int i=u+1;i<=n;i++)
        {
            if (a[i].y<=maxx) continue;
            if (a[i].y<=a[u].y)
            {
                maxx=a[i].y;
                init(i+n,u,1,0);
                init(i,u,MAX,0);
            }
        }
    }
    while (SPFA()==true) 
        lalal();
    printf("%d\n",ans);
    return 0;
}