二维偏序 poj 2352 luogu 3431

这两题都是二维偏序，二维偏序的题怎么去做呢？

回忆一下二维偏序，可以化为这样一个问题：设一个点为（x，y），问（1 ~ x-1，1 ~ y-1）（即一个点的左下方）有多少个点？

那我们一般怎么做呢，就是固定一维（x或y）然后排序，在另一维统计答案。这么做的正确性看看题，手动模拟一下就知道了（严格证明我不知道）。

poj 2352：

很明显题目已经固定了y这一维，然后我们统计x就可以了。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#include<vector>
#include<cmath>
#include<map> 
#include<string>
#include<queue>
#include<stack> 
#include<bitset>
#include<list>
#include<set>
#define IO ios::sync_with_stdio(false)
#define int long long
using namespace std;
string s;
int a[50005],n,c[50005],t,ax[50005],ay[50005],ans[50005];
void add(int x,int y)
{
	for(;x<=42000;x+=x&-x)c[x]+=y;
}
int ask(int x)
{
	int ans=0;
	for(;x;x-=x&-x)ans+=c[x];
	return ans;
}
signed main()
{
	IO;
	cin>>n;
	for(int i=1;i<=n;i++)
	{
		cin>>ax[i]>>ay[i];
		ax[i]++;
		int k=ask(ax[i]);
		//cout<<"i:"<<i<<" "<<"zhi:"<<k<<endl;
		ans[k]++;
		add(ax[i],1);
	}
	for(int i=0;i<=n-1;i++)
	{
		cout<<ans[i]<<endl;
	}
}

luogu 3431

暴力dp很容易，然后我们发现其实就是二维偏序，就用二维偏序去优化dp行了。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#include<vector>
#include<cmath>
#include<map> 
#include<string>
#include<queue>
#include<stack> 
#include<bitset>
#include<list>
#include<set>
#define IO ios::sync_with_stdio(false)
//#define int long long
using namespace std;
struct node
{
	int l,r,sum;
}t[100005*4];
struct dian
{
	int x,y,num;
}a[100005];
int n,m,k,ay[100005],ans=-1;
bool cmp(dian x,dian y)
{
	if(x.x==y.x)
	{
		return x.y<y.y;
	}
	else
	{
		return x.x<y.x;
	}
}
void build(int rt,int l,int r)
{
	t[rt].l=l;
	t[rt].r=r;
	if(t[rt].l==t[rt].r)
	{
		t[rt].sum=0;
		return;
	}
	int mid=l+r>>1;
	build(rt*2,l,mid);
	build(rt*2+1,mid+1,r);
	t[rt].sum=max(t[rt*2].sum,t[rt*2+1].sum);
}
void upd(int rt,int x,int v)
{
	if(t[rt].l==t[rt].r)
	{
		t[rt].sum=v;
		return;
	}
	int mid=t[rt].l+t[rt].r>>1;
	if(x<=mid)
	{
		upd(rt*2,x,v);
	}
	else
	{
		upd(rt*2+1,x,v);
	}
	t[rt].sum=max(t[rt*2].sum,t[rt*2+1].sum);
}
int ask(int rt,int l,int r)
{
	//if(l>r)return 0;
	if(l<=t[rt].l&&r>=t[rt].r)
	{
		return t[rt].sum;
	}
	int mid=t[rt].l+t[rt].r>>1;
	if(l>mid)
	{
		return ask(rt*2+1,l,r);
	}
	else if(r<=mid)
	{
		return ask(rt*2,l,r);
	}
	else
	{
		return max(ask(rt*2,l,r),ask(rt*2+1,l,r));
	}
}
int main()
{
	cin>>n>>m>>k;
	for(int i=1;i<=k;i++)
	{
		scanf("%d %d %d",&a[i].x,&a[i].y,&a[i].num);
	}
	sort(a+1,a+k+1,cmp);
	build(1,1,k);
	for(int i=1;i<=k;i++)
	{
		ay[i]=a[i].y;
	}
	sort(ay+1,ay+k+1);
	int len=unique(ay+1,ay+k+1)-(ay+1);
	for(int i=1;i<=k;i++)
	{
		int pos=lower_bound(ay+1,ay+len+1,a[i].y)-ay;
		int zhi=ask(1,1,pos);
		upd(1,pos,zhi+a[i].num);
		ans=max(ans,zhi+a[i].num);
	}
	cout<<ans;
}