#include <bits/stdc++.h>
#pragma GCC optimize("O3")
#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")
#define int long long
using namespace std;
const int NS = (int)1e6 + 4, mod = (int)1e6 + 3;
int n, k;
int a[NS];
int r[NS], l[NS], chk[NS];
int pw(int x, int y){
if(!y) return 1;
if(y == 1) return x;
int v = pw(x, y / 2);
return v * v % mod * (y % 2 ? x : 1) % mod;
}
int fdr(int x){
return (x == r[x] ? x : r[x] = fdr(r[x]));
}
int fdl(int x){
return (x == l[x] || x == -1 ? x : l[x] = fdl(l[x]));
}
struct Fenwick{
int n;
vector<int> tr;
Fenwick(int m){
n = m + 4;
tr.resize(n);
}
void push(int pos, int val){
pos += 2;
for(int i = pos; i < n; i += (i & -i)){
tr[i] += val;
}
}
int get(int pos){
int rv = 0;
pos += 2;
for(int i = pos; i > 0; i -= (i & -i)){
rv += tr[i];
}
return rv;
}
int get(int l, int r){
return get(r) - get(l - 1);
}
};
signed main(){
ios_base::sync_with_stdio(false);
cin.tie(0);
iota(r, r + NS, 0);
iota(l, l + NS, 0);
memset(chk, -1, sizeof(chk));
int n, k;
cin >> n >> k;
vector<int> a(n - k + 1);
for(int i = 0; i < n - k + 1; ++i){
cin >> a[i];
}
Fenwick tr(n);
auto cant = [&]{
cout << "0\n";
exit(0);
};
auto nor = [&](vector<int> x)->vector<int>{
int nl = fdr(x[0]);
int nr = fdl(x[1]);
if(nl > nr){
int nv = x[2] - tr.get(x[0], x[1]);
if(nv) cant();
return {nl, nr, 0};
}
int nv = x[2] - tr.get(x[0], nl - 1) - tr.get(nr + 1, x[1]);
return {nl, nr, nv};
};
auto upd = [&](int pos, int val){
if(val < 0 || val > 1) cant();
chk[pos] = val;
l[pos] = pos - 1;
r[pos] = pos + 1;
};
vector<vector<int>> stk;
stk.push_back({0, k - 1, a[0]});
for(int i = 1; i < n - k + 1; ++i){
vector<int> now = nor({i, i + k - 1, a[i]});
if(!(int)stk.size()){
upd(now[1], now[2]);
continue;
}
while((int)stk.size() && now[0] <= now[1]){
stk.back() = nor(stk.back());
if(stk.back()[0] > stk.back()[1]){
stk.pop_back();
continue;
}
// cout << i << ' ' << stk.back()[0] << ' ' << stk.back()[1] << ' ' << now[0] << ' ' << now[1] << endl;
assert(stk.back()[0] < now[0] && stk.back()[1] < now[1]);
if(stk.back()[2] == now[2]) break;
if(stk.back()[2] > now[2]){
upd(stk.back()[0], stk.back()[2] - now[2]);
upd(now[1], 0);
}
else{
upd(now[1], now[2] - stk.back()[2]);
upd(stk.back()[0], 0);
}
stk.pop_back();
now = nor(now);
}
if((int)stk.size()){
stk.push_back(now);
}
}
if(!(int)stk.size()){
cout << "1\n";
}
int ac = 0, oc = a[0];
for(int i = 0; i < k; ++i){
ac += (chk[i] == -1);
oc -= (chk[i] == 1);
}
int ans = 1;
for(int i = oc + 1; i <= ac; ++i){
(ans *= i) %= mod;
}
for(int i = 2; i <= ac - oc; ++i){
(ans *= pw(i, mod - 2)) %= mod;
}
cout << ans << '\n';
return 0;
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
9 ms |
23764 KB |
Output is correct |
2 |
Correct |
9 ms |
23776 KB |
Output is correct |
3 |
Incorrect |
9 ms |
23764 KB |
Output isn't correct |
4 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
9 ms |
23764 KB |
Output is correct |
2 |
Correct |
9 ms |
23776 KB |
Output is correct |
3 |
Incorrect |
9 ms |
23764 KB |
Output isn't correct |
4 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
9 ms |
23764 KB |
Output is correct |
2 |
Correct |
9 ms |
23776 KB |
Output is correct |
3 |
Incorrect |
9 ms |
23764 KB |
Output isn't correct |
4 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
9 ms |
23764 KB |
Output is correct |
2 |
Correct |
9 ms |
23776 KB |
Output is correct |
3 |
Incorrect |
9 ms |
23764 KB |
Output isn't correct |
4 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
9 ms |
23764 KB |
Output is correct |
2 |
Correct |
9 ms |
23776 KB |
Output is correct |
3 |
Incorrect |
9 ms |
23764 KB |
Output isn't correct |
4 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
9 ms |
23764 KB |
Output is correct |
2 |
Correct |
9 ms |
23776 KB |
Output is correct |
3 |
Incorrect |
9 ms |
23764 KB |
Output isn't correct |
4 |
Halted |
0 ms |
0 KB |
- |