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[<<],*p1=buf,*p2=buf;
 inline int read(){
     #define num ch-'0'
     char ch;bool flag=;int res;
     while(!isdigit(ch=getc()))
     (ch=='-')&&(flag=true);
     for(res=num;isdigit(ch=getc());res=res*+num);
     (flag)&&(res=-res);
     #undef num
     return res;
 }
 const int N=;
 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+,bx++n),sort(by+,by++n);
     sort(bb1+,bb1++n),sort(bb2+,bb2++n);
     jsx=unique(bx+,bx++n)-bx-;
     jsy=unique(by+,by++n)-by-;
     js1=unique(bb1+,bb1++n)-bb1-;
     js2=unique(bb2+,bb2++n)-bb2-;
     for(int i=;i<=n;++i){
         p[i].x=lower_bound(bx+,bx++jsx,p[i].x)-bx,tx[i]=p[i].x;
         p[i].y=lower_bound(by+,by++jsy,p[i].y)-by,ty[i]=p[i].y;
         p[i].b1=lower_bound(bb1+,bb1++js1,p[i].b1)-bb1;
         p[i].b2=lower_bound(bb2+,bb2++js2,p[i].b2)-bb2;
     }
     sort(p+,p++n,cmpy);
     for(int i=;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+,p++n,cmpx);
     for(int i=;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=;i<sz;++i)
         if(tx[iny[y][i]]>tx[num]) printf("%d ",iny[y][i]);
         for(int i=sz-;~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-;~i;--i)
         if(tx[iny[y][i]]<tx[num]) printf("%d ",iny[y][i]);
         for(int i=;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=,kbj=,sz=iny[y].size();
         for(int i=;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=;
         for(int i=sz-;~i;--i){
             int k=iny[y][i];
             if(kmx>f[k]) f[k]=kmx,bj[k]=kbj;
             if(gup[k]+>f[k]) f[k]=gup[k]+,bj[k]=k;
             if(i++gup[k]>kmx) kmx=i++gup[k],kbj=k;
         }
     }
     printf("%d\n",f[n]-),print(n),puts("");
 }
 int S,T,s,t,ans,tot=;
 int ok[N],du[N],cur[N],head[N],Next[N*],ver[N*],edge[N*],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]=;
 }
 inline void add(int u,int v){
     for(int i=head[u];i;i=Next[i]) if(ver[i]==v) return;
     ok[v]=,++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=;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+;j<sz;++j){
         int k=iny[y][j];
         if(up[k]&&f[up[k]]+j+==f[num]) add(k,up[k]);
         if(lup[k]&&f[lup[k]]+j+==f[num]) add(k,lup[k]);
         if(rup[k]&&f[rup[k]]+j+==f[num]) add(k,rup[k]);
     }
     if(up[num]&&f[up[num]]+==f[num]) add(num,up[num]);
     if(lup[num]&&f[lup[num]]+==f[num]) add(num,lup[num]);
     if(rup[num]&&f[rup[num]]+==f[num]) add(num,rup[num]);
 }
 void build(){
     ok[n]=;
     for(int y=;y<=jsy;++y){
         int sz=iny[y].size();
         for(int i=;i<sz;++i) if(ok[iny[y][i]]) calc(y,i);
     }
 }
 bool bfs(){
     while(!q.empty()) q.pop();
     for(int i=;i<=n+;++i) cur[i]=head[i],dep[i]=-;
     q.push(S),dep[S]=;
     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]<&&edge[i]){
                 dep[v]=dep[u]+,q.push(v);
                 if(v==T) return true;
             }
         }
     }
     return false;
 }
 int dfs(int u,int limit){
     if(u==T||!limit) return limit;
     int flow=,f;
     for(int i=cur[u];i;cur[u]=i=Next[i]){
         int v=ver[i];
         if(dep[v]==dep[u]+&&(f=dfs(v,min(limit,edge[i])))){
             flow+=f,limit-=f;
             edge[i]-=f,edge[i^]+=f;
             if(!limit) break;
         }
     }
     if(!flow) dep[u]=-;
     return flow;
 }
 void work2(){
     build();
     s=n+,t=n+,S=n+,T=n+;
     for(int i=;i<=n;++i)
     du[i]>?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=;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 ;
 }