bzoj4200: [Noi2015]小园丁与老司机（可行流+dp）

传送门

 

这该死的码农题……

题解在这儿->这里

//minamoto
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define inf 0x3f3f3f3f
using namespace std;
#define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
char buf[1<<21],*p1=buf,*p2=buf;
inline int read(){
    #define num ch-'0'
    char ch;bool flag=0;int res;
    while(!isdigit(ch=getc()))
    (ch=='-')&&(flag=true);
    for(res=num;isdigit(ch=getc());res=res*10+num);
    (flag)&&(res=-res);
    #undef num
    return res;
}
const int N=100005;
int n,jsx,jsy,js1,js2;
struct node{int x,y,b1,b2,id;}p[N];
int bx[N],by[N],bb1[N],bb2[N];
int up[N],lup[N],rup[N],f[N],Tx[N],Tb1[N],Tb2[N],bj[N],tx[N],ty[N],gup[N],bup[N];
bool cmpy(node a,node b){return a.y>b.y;}
bool cmpx(node a,node b){return a.x<b.x;}
vector<int> iny[N];
void prework(){
    sort(bx+1,bx+1+n),sort(by+1,by+1+n);
    sort(bb1+1,bb1+1+n),sort(bb2+1,bb2+1+n);
    jsx=unique(bx+1,bx+1+n)-bx-1;
    jsy=unique(by+1,by+1+n)-by-1;
    js1=unique(bb1+1,bb1+1+n)-bb1-1;
    js2=unique(bb2+1,bb2+1+n)-bb2-1;
    for(int i=1;i<=n;++i){
        p[i].x=lower_bound(bx+1,bx+1+jsx,p[i].x)-bx,tx[i]=p[i].x;
        p[i].y=lower_bound(by+1,by+1+jsy,p[i].y)-by,ty[i]=p[i].y;
        p[i].b1=lower_bound(bb1+1,bb1+1+js1,p[i].b1)-bb1;
        p[i].b2=lower_bound(bb2+1,bb2+1+js2,p[i].b2)-bb2;
    }
    sort(p+1,p+1+n,cmpy);
    for(int i=1;i<=n;++i){
        int k=p[i].id;
        up[k]=Tx[p[i].x],lup[k]=Tb1[p[i].b1],rup[k]=Tb2[p[i].b2];
        Tx[p[i].x]=Tb1[p[i].b1]=Tb2[p[i].b2]=k;
    }
    sort(p+1,p+1+n,cmpx);
    for(int i=1;i<=n;++i)
    iny[p[i].y].push_back(p[i].id);
}
void print(int num){
    int y=ty[num],sz=iny[y].size(),nxt=bj[num];
    if(num!=n) printf("%d ",num);
    if(tx[nxt]<tx[num]){
        for(int i=0;i<sz;++i)
        if(tx[iny[y][i]]>tx[num]) printf("%d ",iny[y][i]);
        for(int i=sz-1;~i;--i)
        if(tx[iny[y][i]]<tx[num]&&tx[iny[y][i]]>=tx[nxt]) printf("%d ",iny[y][i]);
    }
    else if(tx[nxt]>tx[num]){
        for(int i=sz-1;~i;--i)
        if(tx[iny[y][i]]<tx[num]) printf("%d ",iny[y][i]);
        for(int i=0;i<sz;++i)
        if(tx[iny[y][i]]>tx[num]&&tx[iny[y][i]]<=tx[nxt]) printf("%d ",iny[y][i]);
    }
    if(bup[nxt]) print(bup[nxt]);
}
void work1(){
    prework();
    for(int y=jsy;y;--y){
        int kmx=0,kbj=0,sz=iny[y].size();
        for(int i=0;i<sz;++i){
            int k=iny[y][i];
            if(up[k]&&f[up[k]]>gup[k]) gup[k]=f[up[k]],bup[k]=up[k];
            if(lup[k]&&f[lup[k]]>gup[k]) gup[k]=f[lup[k]],bup[k]=lup[k];
            if(rup[k]&&f[rup[k]]>gup[k]) gup[k]=f[rup[k]],bup[k]=rup[k];
            f[k]=kmx,bj[k]=kbj;
            if(sz-i+gup[k]>kmx) kmx=sz-i+gup[k],kbj=k;
        }
        kmx=kbj=0;
        for(int i=sz-1;~i;--i){
            int k=iny[y][i];
            if(kmx>f[k]) f[k]=kmx,bj[k]=kbj;
            if(gup[k]+1>f[k]) f[k]=gup[k]+1,bj[k]=k;
            if(i+1+gup[k]>kmx) kmx=i+1+gup[k],kbj=k;
        }
    }
    printf("%d\n",f[n]-1),print(n),puts("");
}
int S,T,s,t,ans,tot=1;
int ok[N],du[N],cur[N],head[N],Next[N*10],ver[N*10],edge[N*10],dep[N];
queue<int> q;
inline void adde(int u,int v,int e){
    ver[++tot]=v,Next[tot]=head[u],head[u]=tot,edge[tot]=e;
    ver[++tot]=u,Next[tot]=head[v],head[v]=tot,edge[tot]=0;
}
inline void add(int u,int v){
    for(int i=head[u];i;i=Next[i]) if(ver[i]==v) return;
    ok[v]=1,++du[v],--du[u],adde(u,v,inf);
}
void calc(int y,int i){
    int sz=iny[y].size(),num=iny[y][i];
    for(int j=0;j<i;++j){
        int k=iny[y][j];
        if(up[k]&&f[up[k]]+sz-j==f[num]) add(k,up[k]);
        if(lup[k]&&f[lup[k]]+sz-j==f[num]) add(k,lup[k]);
        if(rup[k]&&f[rup[k]]+sz-j==f[num]) add(k,rup[k]);
    }
    for(int j=i+1;j<sz;++j){
        int k=iny[y][j];
        if(up[k]&&f[up[k]]+j+1==f[num]) add(k,up[k]);
        if(lup[k]&&f[lup[k]]+j+1==f[num]) add(k,lup[k]);
        if(rup[k]&&f[rup[k]]+j+1==f[num]) add(k,rup[k]);
    }
    if(up[num]&&f[up[num]]+1==f[num]) add(num,up[num]);
    if(lup[num]&&f[lup[num]]+1==f[num]) add(num,lup[num]);
    if(rup[num]&&f[rup[num]]+1==f[num]) add(num,rup[num]);
}
void build(){
    ok[n]=1;
    for(int y=1;y<=jsy;++y){
        int sz=iny[y].size();
        for(int i=0;i<sz;++i) if(ok[iny[y][i]]) calc(y,i);
    }
}
bool bfs(){
    while(!q.empty()) q.pop();
    for(int i=1;i<=n+4;++i) cur[i]=head[i],dep[i]=-1;
    q.push(S),dep[S]=0;
    while(!q.empty()){
        int u=q.front();q.pop();
        for(int i=head[u];i;i=Next[i]){
            int v=ver[i];
            if(dep[v]<0&&edge[i]){
                dep[v]=dep[u]+1,q.push(v);
                if(v==T) return true;
            }
        }
    }
    return false;
}
int dfs(int u,int limit){
    if(u==T||!limit) return limit;
    int flow=0,f;
    for(int i=cur[u];i;cur[u]=i=Next[i]){
        int v=ver[i];
        if(dep[v]==dep[u]+1&&(f=dfs(v,min(limit,edge[i])))){
            flow+=f,limit-=f;
            edge[i]-=f,edge[i^1]+=f;
            if(!limit) break;
        }
    }
    if(!flow) dep[u]=-1;
    return flow;
}
void work2(){
    build();
    s=n+1,t=n+2,S=n+3,T=n+4;
    for(int i=1;i<=n;++i)
    du[i]>0?ans+=du[i],adde(S,i,du[i]):adde(i,T,-du[i]);
    while(bfs()) ans-=dfs(S,inf);
    printf("%d\n",ans);
}
int main(){
    //freopen("testdata.in","r",stdin);
    n=read();
    for(int i=1;i<=n;++i){
        bx[i]=p[i].x=read(),by[i]=p[i].y=read();
        bb1[i]=p[i].b1=p[i].y-p[i].x,bb2[i]=p[i].b2=p[i].x+p[i].y;
        p[i].id=i;
    }
    ++n,p[n].id=n;
    work1(),work2();
    return 0;
}