bzoj D. Captain America（TLE）

题意：

平面上有 n 个点, 第 i 个点的坐标为 $(X_i, Y_i)$, 你需要把每个点染成红色或者蓝色, 染成红色的花费为 r, 染成蓝色的花费为 b。有 m 个限制条件, 有两种类型, 第一种类型为$x = l_i$ 上的红点与蓝点个数差的绝对值不超过 d, 第二种类型为 $y = l_i$ 上的红点与蓝点个数差的绝对值不超过 d。

题意：

肯定是将尽量多的点染成代价较小的颜色。
将坐标离散化，对每一行，每一列开一个点。将对应的点连上，流量为1。
其它限制其实就是上下界，直接建上下界网络流即可。
code：（TLE）

#include<queue>
#include<map>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define LL long long
using namespace std;
const int inf=1<<28;
struct node{
    int X,Y,num,x,y;
}p[100010],P[100010];
int n,m,R,B,numx[100010],numy[100010];
int st,ed,ST,ED,to[100010];
map<int,int> X,Y;
bool cmpx(node a,node b) {return a.X<b.X;}
bool cmpy(node a,node b) {return a.Y<b.Y;}
struct Node{
    int y,c,next;
}a[200010<<3];int len=1,last[200010];
int dx[100010],dy[100010];
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;
}
void ins(int x,int y,int c)
{
    a[++len].y=y;a[len].c=c;
    a[len].next=last[x];last[x]=len;
    a[++len].y=x;a[len].c=0;
    a[len].next=last[y];last[y]=len;
}
int xnum,ynum;
void end() {puts("-1");exit(0);}
void build()
{
    ins(ed,st,inf);
    for(int i=1;i<=n;i++) ins(P[i].x,P[i].y+xnum,1),to[i]=len-1;
    for(int i=1;i<=xnum;i++)
    {
        if(dx[i]==inf) ins(st,i,numx[i]);
        else
        {
            int l=max((numx[i]-dx[i]+1)>>1,0),r=min((numx[i]+dx[i])>>1,numx[i]);
            if(l>r) end();
            ins(st,i,r-l);ins(ST,i,l);ins(st,ED,l);
        }
    }
    for(int i=1;i<=ynum;i++)
    {
        if(dy[i]==inf) ins(i+xnum,ed,numy[i]);
        else
        {
            int l=max((numy[i]-dy[i]+1)>>1,0),r=min((numy[i]+dy[i])>>1,numy[i]);
            if(l>r) end();
            ins(i+xnum,ed,r-l);ins(ST,ed,l);ins(i+xnum,ED,l);
        }
    }
}
int h[200010];
queue<int> q;
int findflow(int x,int ed,int f)
{
    if(x==ed) return f;
    int ans=0,t;
    for(int i=last[x];i;i=a[i].next)
    {
        int y=a[i].y;
        if(ans<f&&h[y]==h[x]+1&&a[i].c)
        {
            ans+=(t=findflow(y,ed,min(a[i].c,f-ans)));
            a[i].c-=t;a[i^1].c+=t;
        }
    }
    if(ans==0) h[x]=0;
    return ans;
}
bool bt_h(int st,int ed)
{
    memset(h,0,sizeof(h));h[st]=1;
    q.push(st);
    while(!q.empty())
    {
        int x=q.front();q.pop();
        for(int i=last[x];i;i=a[i].next)
        {
            int y=a[i].y;
            if(h[y]==0&&a[i].c) h[y]=h[x]+1,q.push(y);
        }
    }
    return h[ed]!=0;
}
int dinic(int st,int ed)
{
    int ans=0;
    while(bt_h(st,ed)) ans+=findflow(st,ed,inf);
    return ans;
}
int flag=false;
int main()
{
    n=read();m=read();R=read();B=read();
    for(int i=1;i<=n;i++) p[i].X=read(),p[i].Y=read(),p[i].num=i;
    sort(p+1,p+n+1,cmpx);xnum=0;
    memset(numx,0,sizeof(numx));
    for(int i=1;i<=n;i++)
    {
        if(i==1||p[i].X!=p[i-1].X) xnum++;
        X[p[i].X]=xnum;P[p[i].num].x=xnum;
        numx[xnum]++;
    }
    sort(p+1,p+n+1,cmpy);ynum=0;
    memset(numy,0,sizeof(numy));
    for(int i=1;i<=n;i++)
    {
        if(i==1||p[i].Y!=p[i-1].Y) ynum++;
        Y[p[i].Y]=ynum;P[p[i].num].y=ynum;
        numy[ynum]++;
    }
    st=xnum+ynum+1;ed=st+1;ST=ed+1;ED=ST+1;
    for(int i=1;i<=xnum;i++) dx[i]=inf;
    for(int i=1;i<=ynum;i++) dy[i]=inf;
    for(int i=1;i<=m;i++)
    {
        int op,l,d;op=read(),l=read(),d=read();
        if(op==1) {int x=X[l];if(x) dx[x]=min(dx[x],d);}
        if(op==2) {int y=Y[l];if(y) dy[y]=min(dy[y],d);}
    }
    build();
    dinic(ST,ED);int ans,tot=0;
    ans=a[3].c;
    for(int i=last[ST];i;i=a[i].next) tot+=a[i].c;
    a[2].c=a[3].c=0;
    if(tot) end();
    ans+=dinic(st,ed);
    if(R>B) swap(R,B),flag=true;
    printf("%lld\n",(LL)R*(LL)ans+(LL)B*(LL)(n-ans));
    for(int i=1;i<=n;i++) putchar((flag&&a[to[i]].c)||(!flag&&!a[to[i]].c)?'r':'b');
}