//Sylwia Sapkowska
#include <bits/stdc++.h>
#pragma GCC optimize("O3", "unroll-loops")
using namespace std;
void __print(int x) {cerr << x;}
void __print(long long x) {cerr << x;}
void __print(long double x) {cerr << x;}
void __print(char x) {cerr << "'" << x << "'";}
void __print(const char *x) {cerr << '"' << x << '"';}
void __print(const string &x) {cerr << '"' << x << '"';}
void __print(bool x) {cerr << (x ? "true" : "false");}
template<typename T, typename V>
void __print(const pair<T, V> &x) {cerr << '{'; __print(x.first); cerr << ", "; __print(x.second); cerr << '}';}
template<typename T>
void __print(const T &x) {int f = 0; cerr << '{'; for (auto &i: x) cerr << (f++ ? ", " : ""), __print(i); cerr << "}";}
void _print() {cerr << "]\n";}
template <typename T, typename... V>
void _print(T t, V... v) {__print(t); if (sizeof...(v)) cerr << ", "; _print(v...);}
#ifdef LOCAL
#define debug(x...) cerr << "[" << #x << "] = ["; _print(x)
#else
#define debug(x...)
#endif
#define int long long
typedef pair<int, int> T;
const int oo = 1e18, oo2 = 1e9+7;
const int mod = 1e9+7;
const int I = (mod+1)/2;
int mul(int a, int b){
return (a*b)%mod;
}
void add(int &a, int b){
a += b;
if (a >= mod) a-=mod;
}
int add2(int a, int b){
a += b;
if (a >= mod) a-=mod;
return a;
}
int power(int a, int b){
if (!b) return 1ll;
int k = power(a, b/2);
k = mul(k, k);
if (b&1) k = mul(k, a);
return k;
}
void solve(){
int n; cin >> n;
vector<int>a(n+1), b(n+1);
vector<int>s = {0, 1};
for (int i = 1; i<=n; i++) {
cin >> a[i] >> b[i];
s.emplace_back(a[i]);
s.emplace_back(b[i]);
}
sort(s.begin(), s.end());
s.erase(unique(s.begin(), s.end()), s.end());
debug(s);
vector<T>seg;
for (int i = 0; i<(int)s.size()-1; i++){
seg.emplace_back(s[i], s[i]);
if (s[i]+1 <= s[i+1]-1) seg.emplace_back(s[i]+1, s[i+1]-1);
}
seg.emplace_back(s.back(), s.back());
debug(seg);
int t = 2*n+2;
vector<int>f(t, 1), inv(t, 1);
for (int i = 1; i<t; i++){
f[i] = mul(f[i-1], i);
inv[i] = mul(inv[i-1], power(i, mod-2));
}
// auto nck = [&](int N, int K){
// if (N < 0 || K < 0 || N < K) return 0ll;
// return mul(f[N], mul(inv[K], inv[N-K]));
// };
int m = (int)seg.size();
//len choose k moze byc bardzo dlugie, ale na szczescie dlugosci len jest O(n), wartosci k jest O(n)
//i mozna uzyc lepszego wzorku by to spreprocessowac
vector nck(m, vector<int>(n+1, 1));
for (int j = 0; j<m; j++){
int len = seg[j].second - seg[j].first+1;
for (int k = 0; k<=n; k++){
for (int rep = 0, val = len; rep < k; rep++, val--){
nck[j][k] = mul(nck[j][k], val);
}
nck[j][k] = mul(nck[j][k], inv[k]);
}
}
//dp[i][j][k] = na ile sposobow da sie wybrac dzialajace rozwiazanie na prefiksie i, gdzie
//ostatnie k elementow nalezy do przedzialu seg[j]
vector dp(m, vector<int>(n+1, 0));
dp[0][1] = 1; //zerowy przedzial
for (int i = 1; i<=n; i++){
auto new_dp = dp;
vector<int>pref(m);
for (int s = 0; s < m; s++){
if (s) pref[s] += pref[s-1];
for (int k = 0; k <= n; k++) add(pref[s], dp[s][k]);
}
// dp[i] = dp[i-1]; //nic nie dokladamy
for (int j = 1; j<m; j++){
if (seg[j].second < a[i] || seg[j].first > b[i]) continue;
int len = seg[j].second - seg[j].first + 1;
// dp[i][j][1] = suma po dp[i-1][<j][cokolwiek]
add(new_dp[j][1], mul(pref[j-1], len)); //dokladamy pierwszy element tego przedzialu
// for (int s = 0; s < j; s++){
// for (int k = 0; k <= n; k++){
// if (dp[i-1][s][k]){
// debug(dp[i-1][s][k], len);
// }
// add(dp[i][j][1], mul(dp[i-1][s][k], len));
// }
// }
for (int k = 2; k <= min(i, len); k++){
//dokladamy juz nie pierwszy element tego przedzialu
// add(new_dp[j][k], mul(dp[j][k-1], mul(power(nck(len, k-1), mod-2), nck(len, k))));
add(new_dp[j][k], mul(dp[j][k-1], mul(power(nck[j][k-1], mod-2), nck[j][k])));
}
}
dp.swap(new_dp);
}
int ans = mod-1;
for (int i = 0; i<m; i++){
for (int j = 0; j<=n; j++){
add(ans, dp[i][j]);
}
}
cout << ans << "\n";
}
int32_t main(){
ios_base::sync_with_stdio(0);
cin.tie(0);
cout.tie(0);
int t = 1;
//cin >> t;
while (t--) solve();
return 0;
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
975 ms |
12596 KB |
Output is correct |
2 |
Correct |
984 ms |
12508 KB |
Output is correct |
3 |
Correct |
1006 ms |
12452 KB |
Output is correct |
4 |
Correct |
967 ms |
12532 KB |
Output is correct |
5 |
Correct |
969 ms |
12556 KB |
Output is correct |
6 |
Correct |
973 ms |
12796 KB |
Output is correct |
7 |
Correct |
969 ms |
12792 KB |
Output is correct |
8 |
Correct |
978 ms |
12880 KB |
Output is correct |
9 |
Correct |
984 ms |
12652 KB |
Output is correct |
10 |
Correct |
991 ms |
12440 KB |
Output is correct |
11 |
Correct |
957 ms |
12520 KB |
Output is correct |
12 |
Correct |
995 ms |
12528 KB |
Output is correct |
13 |
Correct |
977 ms |
12444 KB |
Output is correct |
14 |
Correct |
978 ms |
12696 KB |
Output is correct |
15 |
Correct |
981 ms |
12516 KB |
Output is correct |
16 |
Correct |
167 ms |
3060 KB |
Output is correct |
17 |
Correct |
178 ms |
2736 KB |
Output is correct |
18 |
Correct |
174 ms |
2648 KB |
Output is correct |
19 |
Correct |
182 ms |
3340 KB |
Output is correct |
20 |
Correct |
171 ms |
2668 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
975 ms |
12596 KB |
Output is correct |
2 |
Correct |
984 ms |
12508 KB |
Output is correct |
3 |
Correct |
1006 ms |
12452 KB |
Output is correct |
4 |
Correct |
967 ms |
12532 KB |
Output is correct |
5 |
Correct |
969 ms |
12556 KB |
Output is correct |
6 |
Correct |
973 ms |
12796 KB |
Output is correct |
7 |
Correct |
969 ms |
12792 KB |
Output is correct |
8 |
Correct |
978 ms |
12880 KB |
Output is correct |
9 |
Correct |
984 ms |
12652 KB |
Output is correct |
10 |
Correct |
991 ms |
12440 KB |
Output is correct |
11 |
Correct |
957 ms |
12520 KB |
Output is correct |
12 |
Correct |
995 ms |
12528 KB |
Output is correct |
13 |
Correct |
977 ms |
12444 KB |
Output is correct |
14 |
Correct |
978 ms |
12696 KB |
Output is correct |
15 |
Correct |
981 ms |
12516 KB |
Output is correct |
16 |
Correct |
167 ms |
3060 KB |
Output is correct |
17 |
Correct |
178 ms |
2736 KB |
Output is correct |
18 |
Correct |
174 ms |
2648 KB |
Output is correct |
19 |
Correct |
182 ms |
3340 KB |
Output is correct |
20 |
Correct |
171 ms |
2668 KB |
Output is correct |
21 |
Correct |
1917 ms |
20648 KB |
Output is correct |
22 |
Correct |
1956 ms |
20604 KB |
Output is correct |
23 |
Correct |
1857 ms |
20360 KB |
Output is correct |
24 |
Correct |
1986 ms |
20760 KB |
Output is correct |
25 |
Correct |
1890 ms |
20308 KB |
Output is correct |
26 |
Correct |
1827 ms |
19336 KB |
Output is correct |
27 |
Correct |
1749 ms |
19724 KB |
Output is correct |
28 |
Correct |
1786 ms |
19420 KB |
Output is correct |
29 |
Correct |
1831 ms |
19512 KB |
Output is correct |
30 |
Execution timed out |
2065 ms |
22792 KB |
Time limit exceeded |
31 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
83 ms |
1436 KB |
Output is correct |
2 |
Correct |
67 ms |
1452 KB |
Output is correct |
3 |
Correct |
72 ms |
1456 KB |
Output is correct |
4 |
Correct |
79 ms |
1436 KB |
Output is correct |
5 |
Correct |
86 ms |
1456 KB |
Output is correct |
6 |
Correct |
108 ms |
1452 KB |
Output is correct |
7 |
Correct |
109 ms |
1452 KB |
Output is correct |
8 |
Correct |
108 ms |
1456 KB |
Output is correct |
9 |
Correct |
110 ms |
1672 KB |
Output is correct |
10 |
Correct |
108 ms |
1456 KB |
Output is correct |
11 |
Correct |
73 ms |
1456 KB |
Output is correct |
12 |
Correct |
60 ms |
1456 KB |
Output is correct |
13 |
Correct |
65 ms |
1452 KB |
Output is correct |
14 |
Correct |
67 ms |
1484 KB |
Output is correct |
15 |
Correct |
72 ms |
1456 KB |
Output is correct |
16 |
Correct |
33 ms |
900 KB |
Output is correct |
17 |
Correct |
30 ms |
892 KB |
Output is correct |
18 |
Correct |
33 ms |
892 KB |
Output is correct |
19 |
Correct |
30 ms |
900 KB |
Output is correct |
20 |
Correct |
36 ms |
908 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
975 ms |
12596 KB |
Output is correct |
2 |
Correct |
984 ms |
12508 KB |
Output is correct |
3 |
Correct |
1006 ms |
12452 KB |
Output is correct |
4 |
Correct |
967 ms |
12532 KB |
Output is correct |
5 |
Correct |
969 ms |
12556 KB |
Output is correct |
6 |
Correct |
973 ms |
12796 KB |
Output is correct |
7 |
Correct |
969 ms |
12792 KB |
Output is correct |
8 |
Correct |
978 ms |
12880 KB |
Output is correct |
9 |
Correct |
984 ms |
12652 KB |
Output is correct |
10 |
Correct |
991 ms |
12440 KB |
Output is correct |
11 |
Correct |
957 ms |
12520 KB |
Output is correct |
12 |
Correct |
995 ms |
12528 KB |
Output is correct |
13 |
Correct |
977 ms |
12444 KB |
Output is correct |
14 |
Correct |
978 ms |
12696 KB |
Output is correct |
15 |
Correct |
981 ms |
12516 KB |
Output is correct |
16 |
Correct |
167 ms |
3060 KB |
Output is correct |
17 |
Correct |
178 ms |
2736 KB |
Output is correct |
18 |
Correct |
174 ms |
2648 KB |
Output is correct |
19 |
Correct |
182 ms |
3340 KB |
Output is correct |
20 |
Correct |
171 ms |
2668 KB |
Output is correct |
21 |
Correct |
1917 ms |
20648 KB |
Output is correct |
22 |
Correct |
1956 ms |
20604 KB |
Output is correct |
23 |
Correct |
1857 ms |
20360 KB |
Output is correct |
24 |
Correct |
1986 ms |
20760 KB |
Output is correct |
25 |
Correct |
1890 ms |
20308 KB |
Output is correct |
26 |
Correct |
1827 ms |
19336 KB |
Output is correct |
27 |
Correct |
1749 ms |
19724 KB |
Output is correct |
28 |
Correct |
1786 ms |
19420 KB |
Output is correct |
29 |
Correct |
1831 ms |
19512 KB |
Output is correct |
30 |
Execution timed out |
2065 ms |
22792 KB |
Time limit exceeded |
31 |
Halted |
0 ms |
0 KB |
- |