1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
#include <bits/stdc++.h>
using namespace std;
constexpr int M = 500'005;
constexpr int MOD = 1e9+7;
pair<int, int> tab[M];
vector<vector<int>> g(M);
int doIlu[M];
int odw[M];
long long silnia[M];

void dfs(int v, int ktory){
	doIlu[ktory]++;
	odw[v] = 1;
	for(auto i:g[v]) if(!odw[i]) dfs(i, ktory);
}
long long pote(long long a, long long b){
	long long w = 1;
	while(b>0){
		if (b%2 == 1) w = (a*w)%MOD;
		a = (a*a)%MOD;
		b/=2; 
	}
	return w%MOD;
} 
long long odwrotnosc(long long a){
	return pote(a, MOD-2);
}

int main() {
	cin.tie(0)->sync_with_stdio(0);
	int n;
	cin>>n;
	silnia[0] = 1;
	for(int i=1; i<=n; i++) silnia[i] = (silnia[i-1]*i) % MOD;
	
	for(int i=1; i<=n; i++) cin>>tab[i].first>>tab[i].second;
	vector<pair<int, int>> v;
	for(int i=n; i>=1; i--){
		for(int j=0; j<v.size(); j++) if(v[j].first >= tab[i].first && v[j].first <= tab[i].second){
			g[n-i+1].push_back(v[j].second);
			v.erase(v.begin()+j);
			j--;
		}	
		v.push_back({tab[i].first, n-i+1});
		v.push_back({tab[i].second, n-i+1});
	}
	for(int i=1; i<=n; i++){
		for(int j=1; j<=n+1; j++) odw[j] = 0;
		doIlu[i]--;
		dfs(i, i);
	}
	long long suma = 0;
	for(int i=1; i<=n; i++){
		suma = suma + (odwrotnosc(doIlu[i] + 1)) %MOD;
	}
	cout<<suma%MOD<<endl;
	return 0;
}