#include <cstdio>
#include <algorithm>
using namespace std;
const int M = 1005;
const int MOD = 1e9+7;
#define int long long
int read()
{
	int x=0,f=1;char c;
	while((c=getchar())<'0' || c>'9') {if(c=='-') f=-1;}
	while(c>='0' && c<='9') {x=(x<<3)+(x<<1)+(c^48);c=getchar();}
	return x*f;
}
int n,m,ans,a[M],b[M],t[M],inv[M],sum[M],dp[M][M];
signed main()
{
	n=read();inv[0]=inv[1]=1;
	for(int i=2;i<=n;i++)
		inv[i]=inv[MOD%i]*(MOD-MOD/i)%MOD;
	for(int i=1;i<=n;i++)
	{
		a[i]=read();b[i]=read()+1;
		t[++m]=a[i];t[++m]=b[i];
	}
	sort(t+1,t+1+m);
	m=unique(t+1,t+1+m)-t-1;
	for(int i=1;i<=n;i++)
	{
		a[i]=lower_bound(t+1,t+1+m,a[i])-t;
		b[i]=lower_bound(t+1,t+1+m,b[i])-t;
	}
	dp[0][0]=1;
	for(int i=0;i<=m;i++) sum[i]=1;
	for(int i=1;i<=n;i++)
	{
		for(int j=a[i];j<b[i];j++)
		{
			int l=t[j+1]-t[j];
			for(int k=i;k>=2;k--)
				dp[j][k]=(dp[j][k]+(l-k+1)
				*inv[k]%MOD*dp[j][k-1])%MOD;
			dp[j][1]=(dp[j][1]+sum[j-1]*l)%MOD;
		}
		for(int j=1;j<m;j++)
		{
			sum[j]=sum[j-1];
			for(int k=1;k<=i;k++)
				sum[j]=(sum[j]+dp[j][k])%MOD;
		}
	}
	printf("%lld\n",sum[m-1]-1);
}

E99 线性DP+前缀和优化 P3643 [APIO2016] 划艇