【NOIP 校内模拟】行星通道计划（二维树状数组）

两个管道相交是啥情况？

有环展链 我们把环展成链过后 发现只会出现如下图两种情况

然后我们维护二维BIT1：表示左端点小于等于x 右端点小于等于y的个数

BIT2：左端点大于等于x 右端点小于等于y

BIT3：左端点大于等于x 右端点大于等于y

查询就很简单了

对于操作1相当于一个更新

#include<bits/stdc++.h>
#define N 1005
#define M 500005
using namespace std;
int n,m,tree1[N][N],tree2[N][N],tree3[N][N];
inline int read(){
	int ans=0;
	char ch=getchar();
	while(!isdigit(ch))ch=getchar();
	while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
	return ans;
}
inline int lowbit(int x)
{
	return x&(-x);
}
inline int query1(int x,int y)	//左端点小于等于x 右端点小于等于y
{
	int ans=0;
	for(int i=x;i;i-=lowbit(i))
		for(int j=y;j;j-=lowbit(j))
			ans+=tree1[i][j];
	return ans;
} 
inline int query2(int x,int y)	//左端点大于等于x 右端点小于等于y
{
	int ans=0;
	for(int i=x;i<=n;i+=lowbit(i))
		for(int j=y;j;j-=lowbit(j))
			ans+=tree2[i][j]; 
	return ans;
} 
inline int query3(int x,int y)	//左端点大于x 右端点大于y
{
	int ans=0;
	for(int i=x;i<=n;i+=lowbit(i))
		for(int j=y;j<=n;j+=lowbit(j))
			ans+=tree3[i][j]; 
	return ans;
} 
inline void update1(int x,int y,int v)
{
	for(int i=x;i<=n;i+=lowbit(i))
		for(int j=y;j<=n;j+=lowbit(j))
			tree1[i][j]+=v;
}
inline void update2(int x,int y,int v)
{
	for(int i=x;i;i-=lowbit(i))
		for(int j=y;j<=n;j+=lowbit(j))
			tree2[i][j]+=v;
}
inline void update3(int x,int y,int v)
{
	for(int i=x;i;i-=lowbit(i))
		for(int j=y;j;j-=lowbit(j))
			tree3[i][j]+=v; 
}
struct S
{
	int x,y;
}guan[M];
int main()
{
	cin>>n>>m;
	for(int i=1;i<=m;i++)
	{
		int opt;
		opt=read();
		if(opt==0)
		{
			int x,y;
			x=read();	y=read();
			if(x>y)	swap(x,y);
			int ansl=query1(x,y)-query2(1,x);
			int ansr=query3(x,y)-query2(y,n);
			printf("%d\n",ansl+ansr);
			update1(x,y,1);	update2(x,y,1);	update3(x,y,1);
			guan[i].x=x;	guan[i].y=y;
		}
		if(opt==1)
		{
			int pos;
			pos=read();
			int x=guan[pos].x;
			int y=guan[pos].y;
			update1(x,y,-1);	update2(x,y,-1);	update3(x,y,-1);
		}
	}
}