BZOJ 3218(a + b Problem-二分图套值域线段树)




出这题的人是怎么想出来的……

言归正传,这题是二分图套值域线段树。

首先经过 @Vfleaking的神奇建图后,把图拆成二分图,

不妨利用有向图最小割的性质建图(以前我一直以为最小割和边的方向无关,可这样的话很奇怪哦……)

理解悲剧……

我们可以利用边有向的性质解决黑白色块……

然后发现线段树很多……主席树闪亮登场

然后·就这麽A了?…………

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<functional>
#include<cmath>
#include<cctype>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define RepD(i,n) for(int i=n;i>=0;i--)
#define MEM(a) memset(a,0,sizeof(a))
#define MEMI(a) memset(a,127,sizeof(a))
#define MEMi(a) memset(a,128,sizeof(a))
#define INF (2139062143)
#define MAXN (500000 +10)
#define MAXM (500000 +10)
int n,m,s,t;
int edge[MAXM],next[MAXM]={0},pre[MAXN]={0},weight[MAXM],size=1;
void addedge(int u,int v,int w)
{
	edge[++size]=v;
	weight[size]=w;
	next[size]=pre[u];
	pre[u]=size;
}
void addedge2(int u,int v,int w){/*cout<<u<<' '<<v<<' '<<w<<endl;*/addedge(u,v,w),addedge(v,u,0);}
int cnt[MAXN]={0},d[MAXN]={0};
long long totflow=0;
long long sap(int x,int flow)
{
	if (x==t) return flow;
	int nowflow=0;
	Forp(x)
	{
		int &v=edge[p];
		if (d[v]==d[x]-1&&weight[p])
		{
			long long fl=sap(v,min(weight[p],flow));
			weight[p]-=fl,weight[p^1]+=fl,nowflow+=fl,flow-=fl;
			if (!flow) return nowflow;
		}
	}	
	if (!(--cnt[d[x]++])) d[s]=t+1;
	cnt[d[x]]++;
	return nowflow;
}
struct node
{
	node *ch[2];
	node(){ch[0]=ch[1]=NULL;}
}q[MAXN],*root[MAXN]={NULL};
int tail=0;
void ins(node *&p,node *x,int l,int r,int c)
{
	int m=(l+r) >>1;
	p=&q[++tail];
	if (x) *p=*x;
	if (l==r) return;
	if (c<=m) ins(p->ch[0],p->ch[0],l,m,c);
	else ins(p->ch[1],p->ch[1],m+1,r,c);		
}
void print(node *p)
{
	//cout<<
}
node *ans[MAXN];
int ans_siz=0;
void qur(node *p,int l,int r,int L,int R)
{
	if (!p||L>R) return;
	int m=l+r >>1;
	if (L<=l&&r<=R) {ans[++ans_siz]=p;return;}
	if (l==r) return;
	if (L<=m) qur(p->ch[0],l,m,L,R);
	if (m<R) qur(p->ch[1],m+1,r,L,R);
}
int a2[MAXN],a[MAXN],b[MAXN],w[MAXN],l[MAXN],r[MAXN],p[MAXN];
int main()
{
	freopen("bzoj3218.in","r",stdin);
//	freopen(".out","w",stdout);
	scanf("%d",&n);
	For(i,n) scanf("%d%d%d%d%d%d",&a[i],&b[i],&w[i],&l[i],&r[i],&p[i]);
	memcpy(a2,a,sizeof(a));
	sort(a2+1,a2+n+1);
	int size=unique(a2+1,a2+n+1)-(a2+1);
	For(i,n) 
	{
		a[i]=lower_bound(a2+1,a2+size+1,a[i])-(a2);
		l[i]=lower_bound(a2+1,a2+size+1,l[i])-(a2);
		r[i]=upper_bound(a2+1,a2+size+1,r[i])-(a2)-1;
	}
	For(i,n) ins(root[i],root[i-1],1,size,a[i]);
	
	s=2*n+1; t=2*n+2;
	For(i,n) addedge2(s,i,b[i]),addedge2(i,t,w[i]),addedge2(i,n+i,p[i]);
	For(i,tail)
	{
		if (q[i].ch[0]) addedge2(t+i,t+((q[i].ch[0])-q),INF);
		if (q[i].ch[1]) addedge2(t+i,t+((q[i].ch[1])-q),INF);
	}
	For(i,n)
	{
		ans_siz=0;
		qur(root[i],1,size,a[i],a[i]);
		qur(root[i-1],1,size,a[i],a[i]);
		node *cur=ans[1],*last=NULL;
		if (ans_siz>1) last=ans[2];
		//cout<<t+((cur)-(q))<<'*';
		addedge2(t+((cur)-(q)),i,INF);
		
		if (ans_siz>1) addedge2(t+(cur-(q)),t+((last)-(q)),INF);
	}
	For(i,n)
	{
		ans_siz=0;
		qur(root[i],1,size,l[i],r[i]);
		For(j,ans_siz) addedge2(n+i,t+((ans[j])-(q)),INF);
	}	
	//cout<<2*n+2+tail<<endl;	
	cnt[0]=2*n+2+tail;
	while (d[s]<=t) {totflow+=sap(s,INF);} 
	//For(i,2*n+tail) cout<<d[i]<<' ';
	long long ans=-totflow;
	For(i,n) ans+=b[i]+w[i];
	cout/*<<totflow<<' '*/<<ans<<endl;	
	
	return 0;
}






原文地址:https://www.cnblogs.com/snake-hand/p/3144807.html