noip 2018 雅礼模拟赛day1 T2

玄学dp...

看见第一个条件会非常自然地按y排序，然后原地炸飞...

（当然，zwz大佬除外...）

因为第一个条件的可优化点非常少，同时卡时间和空间，所以并不是很好。

所以我们先按x来排序，设计状态dp[i][1]表示从i向右射出一条光线，dp[i][0]表示从i向左射出一条光线

然后考虑转移：

我们从前向后枚举每个点，然后再枚举要射向哪个点，这样时间复杂度控制在O（n^2）

因为我们排过序，所以i以前的点的横坐标均小于i的横坐标，所以我们只需讨论y的大小即可

分两类讨论：

①：yi>yj，说明可以从i向左下方向j射出一条光线，那么更新：dp[i][0]+=dp[j][1]

②：yi<yj，说明可以从j向右下方射出一条光线，那么更新：dp[j][1]+=dp[i][0]

#include <cstdio>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <queue>
#include <stack>
#define ll long long
#define mode 1000000007
using namespace std;
ll dp[6005][2];
int n;
struct node
{
	int x,y;
}p[6005];
bool cmp(node a,node b)
{
	return a.x<b.x;
}
int main()
{
	freopen("refract.in","r",stdin);
	freopen("refract.out","w",stdout);
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	{
		scanf("%d%d",&p[i].x,&p[i].y);
	}
	sort(p+1,p+n+1,cmp);
	for(int i=1;i<=n;i++)
	{
		dp[i][0]=dp[i][1]=1;
		for(int j=i-1;j>=1;j--)
		{
			if(p[i].y>p[j].y)
			{
				dp[i][0]+=dp[j][1];
				dp[i][0]%=mode;
			}else
			{
				dp[j][1]+=dp[i][0];
				dp[j][1]%=mode;
			}
		}
	}
	ll ret=mode-n;
	for(int i=1;i<=n;i++)
	{
		ret+=dp[i][1]+dp[i][0];
		ret%=mode;
	}
	printf("%I64d\n",ret);
	return 0;
}