#include<bits/stdc++.h>
using namespace std;
const int nax = 503;
const int mod = 1e9 + 7;
const int cax = 1e6 + 3;
const long long sux = 1361;
vector<int>ind[nax*2];
int a[nax], b[nax];
int n, m;
int dp[nax*2][nax];
int inv[nax];
vector<int>imp;
int add(int a, int b) {
a += b;
if(a>=mod) a -= mod;
return a;
}
int mul(int a, int b) {
return (a * 1LL * b) % mod;
}
int quickpow(int a, int b) {
int ret = 1;
while(b) {
if(b&1) ret = mul(ret, a);
a = mul(a, a);
b >>= 1;
}
return ret;
}
int get(int x) {
int &ret = inv[x];
if(ret==-1) {
ret = quickpow(x, mod-2);
}
return ret;
}
vector<pair<long long,int>>table[cax];
int C(int n, int r) {
long long H = n * sux + r;
int H2 = H % cax;
for(auto x : table[H2]) {
if(x.first==H) return x.second;
}
int ret = 1;
for(int i=n-r+1; i<=n; ++i) {
ret = mul(ret, i);
}
for(int i=1; i<=r; ++i) {
ret = mul(ret, get(i));
}
table[H2].emplace_back(H, ret);
assert((int)table[H2].size()<=3);
return ret;
}
vector<pair<long long,int>>andp[cax];
int get(int len, int op) {
long long H = len * sux + op;
int H2 = H % cax;
for(auto x : andp[H2]) {
if(x.first==H) return x.second;
}
int ret = 0;
for(int take=0; take<=min(op, len-1); ++take) {
ret = add(ret, mul(C(op, take), C(len, take+1)));
}
andp[H2].emplace_back(H, ret);
assert((int)andp[H2].size()<=3);
return ret;
}
int solve(int i, int prv) {
if(i==m-1) return prv>0;
int &ret = dp[i][prv];
if(ret!=-1) return ret;
ret = 0;
ret = solve(i+1, prv);
int len = imp[i+1]-imp[i];
for(int mx=0, cur_op=0; mx<ind[i].size(); ++mx) {
if(ind[i][mx]<=prv) continue;
int cur = mul(solve(i+1, ind[i][mx]), get(len, cur_op));
cur_op++;
ret = add(ret, cur);
}
return ret;
}
void PlayGround() {
cin>>n;
set<int>pts;
for(int i=1; i<=n; ++i) {
cin>>a[i]>>b[i];
pts.insert(a[i]);
pts.insert(b[i]+1);
}
for(int x : pts) {
imp.push_back(x);
}
m = imp.size();
for(int i=0; i<m; ++i) {
for(int j=1; j<=n; ++j) {
if(imp[i]>=a[j] && imp[i]<=b[j]) {
ind[i].push_back(j);
}
}
}
memset(inv, -1, sizeof inv);
memset(dp, -1, sizeof dp);
cout<<solve(0, 0)<<'\n';
// cerr << "Time elapsed: " << 1.0 * clock() / CLOCKS_PER_SEC << " s.\n";
}
int main() {
ios_base::sync_with_stdio(false);
cin.tie(NULL);
PlayGround();
return 0;
}
Compilation message
boat.cpp: In function 'int solve(int, int)':
boat.cpp:88:29: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
88 | for(int mx=0, cur_op=0; mx<ind[i].size(); ++mx) {
| ~~^~~~~~~~~~~~~~
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
28 ms |
49356 KB |
Output is correct |
2 |
Correct |
28 ms |
49456 KB |
Output is correct |
3 |
Correct |
29 ms |
49428 KB |
Output is correct |
4 |
Correct |
28 ms |
49348 KB |
Output is correct |
5 |
Correct |
29 ms |
49448 KB |
Output is correct |
6 |
Correct |
30 ms |
49520 KB |
Output is correct |
7 |
Correct |
28 ms |
49356 KB |
Output is correct |
8 |
Correct |
28 ms |
49472 KB |
Output is correct |
9 |
Correct |
28 ms |
49388 KB |
Output is correct |
10 |
Correct |
33 ms |
49412 KB |
Output is correct |
11 |
Correct |
28 ms |
49440 KB |
Output is correct |
12 |
Correct |
29 ms |
49360 KB |
Output is correct |
13 |
Correct |
35 ms |
49424 KB |
Output is correct |
14 |
Correct |
30 ms |
49392 KB |
Output is correct |
15 |
Correct |
28 ms |
49364 KB |
Output is correct |
16 |
Correct |
24 ms |
49236 KB |
Output is correct |
17 |
Correct |
26 ms |
49272 KB |
Output is correct |
18 |
Correct |
28 ms |
49232 KB |
Output is correct |
19 |
Correct |
25 ms |
49236 KB |
Output is correct |
20 |
Correct |
27 ms |
49336 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
28 ms |
49356 KB |
Output is correct |
2 |
Correct |
28 ms |
49456 KB |
Output is correct |
3 |
Correct |
29 ms |
49428 KB |
Output is correct |
4 |
Correct |
28 ms |
49348 KB |
Output is correct |
5 |
Correct |
29 ms |
49448 KB |
Output is correct |
6 |
Correct |
30 ms |
49520 KB |
Output is correct |
7 |
Correct |
28 ms |
49356 KB |
Output is correct |
8 |
Correct |
28 ms |
49472 KB |
Output is correct |
9 |
Correct |
28 ms |
49388 KB |
Output is correct |
10 |
Correct |
33 ms |
49412 KB |
Output is correct |
11 |
Correct |
28 ms |
49440 KB |
Output is correct |
12 |
Correct |
29 ms |
49360 KB |
Output is correct |
13 |
Correct |
35 ms |
49424 KB |
Output is correct |
14 |
Correct |
30 ms |
49392 KB |
Output is correct |
15 |
Correct |
28 ms |
49364 KB |
Output is correct |
16 |
Correct |
24 ms |
49236 KB |
Output is correct |
17 |
Correct |
26 ms |
49272 KB |
Output is correct |
18 |
Correct |
28 ms |
49232 KB |
Output is correct |
19 |
Correct |
25 ms |
49236 KB |
Output is correct |
20 |
Correct |
27 ms |
49336 KB |
Output is correct |
21 |
Correct |
327 ms |
50640 KB |
Output is correct |
22 |
Correct |
316 ms |
50928 KB |
Output is correct |
23 |
Correct |
305 ms |
50752 KB |
Output is correct |
24 |
Correct |
289 ms |
50636 KB |
Output is correct |
25 |
Correct |
331 ms |
50664 KB |
Output is correct |
26 |
Correct |
575 ms |
51392 KB |
Output is correct |
27 |
Correct |
585 ms |
51164 KB |
Output is correct |
28 |
Correct |
572 ms |
51160 KB |
Output is correct |
29 |
Correct |
570 ms |
51320 KB |
Output is correct |
30 |
Correct |
32 ms |
49492 KB |
Output is correct |
31 |
Correct |
34 ms |
49480 KB |
Output is correct |
32 |
Correct |
33 ms |
49492 KB |
Output is correct |
33 |
Correct |
34 ms |
49436 KB |
Output is correct |
34 |
Correct |
32 ms |
49384 KB |
Output is correct |
35 |
Correct |
37 ms |
49440 KB |
Output is correct |
36 |
Correct |
38 ms |
49416 KB |
Output is correct |
37 |
Correct |
33 ms |
49492 KB |
Output is correct |
38 |
Correct |
33 ms |
49424 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
36 ms |
49748 KB |
Output is correct |
2 |
Correct |
34 ms |
49740 KB |
Output is correct |
3 |
Correct |
30 ms |
49716 KB |
Output is correct |
4 |
Correct |
29 ms |
49748 KB |
Output is correct |
5 |
Correct |
32 ms |
49816 KB |
Output is correct |
6 |
Correct |
41 ms |
50124 KB |
Output is correct |
7 |
Correct |
46 ms |
50008 KB |
Output is correct |
8 |
Correct |
40 ms |
50124 KB |
Output is correct |
9 |
Correct |
43 ms |
50124 KB |
Output is correct |
10 |
Correct |
40 ms |
50124 KB |
Output is correct |
11 |
Correct |
33 ms |
49828 KB |
Output is correct |
12 |
Correct |
30 ms |
49672 KB |
Output is correct |
13 |
Correct |
32 ms |
49748 KB |
Output is correct |
14 |
Correct |
32 ms |
49672 KB |
Output is correct |
15 |
Correct |
32 ms |
49772 KB |
Output is correct |
16 |
Correct |
28 ms |
49496 KB |
Output is correct |
17 |
Correct |
27 ms |
49416 KB |
Output is correct |
18 |
Correct |
27 ms |
49544 KB |
Output is correct |
19 |
Correct |
26 ms |
49428 KB |
Output is correct |
20 |
Correct |
28 ms |
49484 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
28 ms |
49356 KB |
Output is correct |
2 |
Correct |
28 ms |
49456 KB |
Output is correct |
3 |
Correct |
29 ms |
49428 KB |
Output is correct |
4 |
Correct |
28 ms |
49348 KB |
Output is correct |
5 |
Correct |
29 ms |
49448 KB |
Output is correct |
6 |
Correct |
30 ms |
49520 KB |
Output is correct |
7 |
Correct |
28 ms |
49356 KB |
Output is correct |
8 |
Correct |
28 ms |
49472 KB |
Output is correct |
9 |
Correct |
28 ms |
49388 KB |
Output is correct |
10 |
Correct |
33 ms |
49412 KB |
Output is correct |
11 |
Correct |
28 ms |
49440 KB |
Output is correct |
12 |
Correct |
29 ms |
49360 KB |
Output is correct |
13 |
Correct |
35 ms |
49424 KB |
Output is correct |
14 |
Correct |
30 ms |
49392 KB |
Output is correct |
15 |
Correct |
28 ms |
49364 KB |
Output is correct |
16 |
Correct |
24 ms |
49236 KB |
Output is correct |
17 |
Correct |
26 ms |
49272 KB |
Output is correct |
18 |
Correct |
28 ms |
49232 KB |
Output is correct |
19 |
Correct |
25 ms |
49236 KB |
Output is correct |
20 |
Correct |
27 ms |
49336 KB |
Output is correct |
21 |
Correct |
327 ms |
50640 KB |
Output is correct |
22 |
Correct |
316 ms |
50928 KB |
Output is correct |
23 |
Correct |
305 ms |
50752 KB |
Output is correct |
24 |
Correct |
289 ms |
50636 KB |
Output is correct |
25 |
Correct |
331 ms |
50664 KB |
Output is correct |
26 |
Correct |
575 ms |
51392 KB |
Output is correct |
27 |
Correct |
585 ms |
51164 KB |
Output is correct |
28 |
Correct |
572 ms |
51160 KB |
Output is correct |
29 |
Correct |
570 ms |
51320 KB |
Output is correct |
30 |
Correct |
32 ms |
49492 KB |
Output is correct |
31 |
Correct |
34 ms |
49480 KB |
Output is correct |
32 |
Correct |
33 ms |
49492 KB |
Output is correct |
33 |
Correct |
34 ms |
49436 KB |
Output is correct |
34 |
Correct |
32 ms |
49384 KB |
Output is correct |
35 |
Correct |
37 ms |
49440 KB |
Output is correct |
36 |
Correct |
38 ms |
49416 KB |
Output is correct |
37 |
Correct |
33 ms |
49492 KB |
Output is correct |
38 |
Correct |
33 ms |
49424 KB |
Output is correct |
39 |
Correct |
36 ms |
49748 KB |
Output is correct |
40 |
Correct |
34 ms |
49740 KB |
Output is correct |
41 |
Correct |
30 ms |
49716 KB |
Output is correct |
42 |
Correct |
29 ms |
49748 KB |
Output is correct |
43 |
Correct |
32 ms |
49816 KB |
Output is correct |
44 |
Correct |
41 ms |
50124 KB |
Output is correct |
45 |
Correct |
46 ms |
50008 KB |
Output is correct |
46 |
Correct |
40 ms |
50124 KB |
Output is correct |
47 |
Correct |
43 ms |
50124 KB |
Output is correct |
48 |
Correct |
40 ms |
50124 KB |
Output is correct |
49 |
Correct |
33 ms |
49828 KB |
Output is correct |
50 |
Correct |
30 ms |
49672 KB |
Output is correct |
51 |
Correct |
32 ms |
49748 KB |
Output is correct |
52 |
Correct |
32 ms |
49672 KB |
Output is correct |
53 |
Correct |
32 ms |
49772 KB |
Output is correct |
54 |
Correct |
28 ms |
49496 KB |
Output is correct |
55 |
Correct |
27 ms |
49416 KB |
Output is correct |
56 |
Correct |
27 ms |
49544 KB |
Output is correct |
57 |
Correct |
26 ms |
49428 KB |
Output is correct |
58 |
Correct |
28 ms |
49484 KB |
Output is correct |
59 |
Correct |
962 ms |
61496 KB |
Output is correct |
60 |
Runtime error |
694 ms |
119144 KB |
Execution killed with signal 6 |
61 |
Halted |
0 ms |
0 KB |
- |