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