https://blog.csdn.net/litble/article/details/80463466
https://www.cnblogs.com/cjyyb/p/9288367.html
一上午就弄这道题 犯了很多错误 比如排序排错 没更新啥的。。
大概是把\(y\)从大到小dp \(dp_i\)表示从\(i\)向上或斜向上走最多多少点
\(g_i\)表示假如\(i\)是这一层第一个到达的点从\(i\)开始走最多多少点 可以在这一层走一遍
同一层内从\(i\)开始到\(j\)离开该层的话会先绕到该层端点再到\(j\)
当然也可以直接从\(i\)离开该层
第二问需要上下界最小流
litble给出了一个神奇的针对这个图的上下界最小流方法 太神了
我实现的好麻烦。。
#include
using namespace std;
#define fp(i,l,r) for(register int (i)=(l);i<=(r);++(i))
#define fd(i,l,r) for(register int (i)=(l);i>=(r);--(i))
#define fe(i,u) for(register int (i)=front[(u)];(i);(i)=e[(i)].next)
#define mem(a) memset((a),0,sizeof (a))
#define O(x) cerr<<#x<<':'<=10)wr(x/10);
putchar('0'+x%10);
}
const int MAXN=50020,inf=1e9+1e7;
struct poi{
int x,y,id;
}p[MAXN],sta[MAXN];
inline bool cmp1(poi a,poi b){
if(a.x!=b.x)return a.xb.y;
return a.xmp;
int n,frm[MAXN],g[MAXN],dp[MAXN],to[MAXN][5],top,pos[MAXN],grm[MAXN],gp[MAXN];
void output(int x){
if(!x)return;
if(p[gp[grm[x]]].y!=p[gp[x]].y){printf("%d ",x),output(grm[x]);return;}
top=0;int l=gp[x],r=gp[x],p1=-1,p2=-1;
while(l>1&&p[l-1].y==p[gp[x]].y)--l;while(r+1<=n&&p[r+1].y==p[gp[x]].y)++r;
fp(i,l,r)sta[++top]=p[i];
fp(i,1,top)if(sta[i].id==x)p1=i;
fp(i,1,top)if(sta[i].id==grm[x])p2=i;
assert(p1!=-1);assert(p2!=-1);
if(p2>p1){
fd(i,p1,1)printf("%d ",sta[i].id);
fp(i,p1+1,p2)printf("%d ",sta[i].id);
}
else{
fp(i,p1,top)printf("%d ",sta[i].id);
fd(i,p1-1,p2)printf("%d ",sta[i].id);
}
output(frm[grm[x]]);
}
inline void solve1(){
sort(p+1,p+1+n,cmp1);
fp(i,1,n-1)if(p[i].x==p[i+1].x)to[p[i].id][2]=p[i+1].id;
fp(i,1,n)if(p[i].x==0){to[0][2]=p[i].id;break;}
mp.clear();
fp(i,1,n){
int t=p[i].x+p[i].y;
if(mp[t])to[p[i].id][3]=mp[t];mp[t]=p[i].id;
if(!t)to[0][3]=p[i].id;
}
mp.clear();
fd(i,n,1){
int t=p[i].x-p[i].y;
if(mp[t])to[p[i].id][4]=mp[t];mp[t]=p[i].id;
if(!t)to[0][4]=p[i].id;
}
sort(p+1,p+1+n,cmp2);
for(int l=1,r;l<=n;l=r+1){
r=l;while(rg[sta[i].id])g[sta[i].id]=t,grm[sta[i].id]=sta[p].id;
if(dp[sta[i].id]+top-i>t)p=i;
}
p=top;
fd(i,top-1,1){
int t=dp[sta[p].id]+p-1;
if(t>g[sta[i].id])g[sta[i].id]=t,grm[sta[i].id]=sta[p].id;
if(tbuc;buc.clear();
int t=dp[sta[1].id]+top-1;buc.push_back(sta[1].id);
fp(i,2,top){
if(visg[sta[i].id]&&g[sta[i].id]==t){
for(int t:buc)visf[t]=1;buc.clear();
}
int x=dp[sta[i].id]+top-i;
if(x>t)t=x,buc.clear(),buc.push_back(sta[i].id);
else if(x==t)buc.push_back(sta[i].id);
}
buc.clear();
buc.push_back(sta[top].id);t=dp[sta[top].id]+top-1;
fd(i,top-1,1){
if(visg[sta[i].id]&&g[sta[i].id]==t){
for(int t:buc)visf[t]=1;buc.clear();
}
int x=dp[sta[i].id]+i-1;
if(x>t)t=x,buc.clear(),buc.push_back(sta[i].id);
else if(x==t)buc.push_back(sta[i].id);
}
fp(i,1,top)if(visf[sta[i].id]){
int t=sta[i].id;
fp(j,2,4){
if(g[to[t][j]]+1==dp[t])add(t,to[t][j]),visg[to[t][j]]=1;
}
}
}
S=n+1;T=n+2;SS=n+3;TT=n+4;
fp(i,0,n)if(visg[i]||visf[i]){
adde(S,i,inf);adde(i,T,inf);
}
fp(i,0,n)if(deg[i]){
if(deg[i]>0)adde(SS,i,deg[i]);
else adde(i,TT,-deg[i]);
}
while(bfs()){
memcpy(head,front,sizeof head);dfs(SS,inf);
}
adde(T,S,inf);
while(bfs()){
memcpy(head,front,sizeof head);ans+=dfs(SS,inf);
}
printf("%d\n",ans);
}
main(){
n=read();
fp(i,1,n)p[i].x=read(),p[i].y=read(),p[i].id=i;
solve1();solve2();
return 0;
}