玄學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;
}