#include "circuit.h"
#include <bits/stdc++.h>
using namespace std;
#define int long long
#define FOR(i,s,e) for(int i = s; i <= (int)e; ++i)
#define DEC(i,s,e) for(int i = s; i >= (int)e; --i)
#define IAMSPEED ios_base::sync_with_stdio(false); cin.tie(0);
#ifdef LOCAL
#define db(x) cerr << #x << "=" << x << "\n"
#define db2(x, y) cerr << #x << "=" << x << " , " << #y << "=" << y << "\n"
#define db3(a,b,c) cerr<<#a<<"="<<a<<","<<#b<<"="<<b<<","<<#c<<"="<<c<<"\n"
#define dbv(v) cerr << #v << ":"; for (auto ite : v) cerr << ite << ' '; cerr <<"\n"
#define dbvp(v) cerr << #v << ":"; for (auto ite : v) cerr << "{" << ite.f << ',' << ite.s << "} "; cerr << "\n"
#define dba(a,ss,ee) cerr << #a << ":"; FOR(ite,ss,ee) cerr << a[ite] << ' '; cerr << "\n"
#define reach cerr << "LINE: " << __LINE__ << "\n";
#else
#define reach
#define db(x)
#define db2(x,y)
#define db3(a,b,c)
#define dbv(v)
#define dbvp(v)
#define dba(a,ss,ee)
#endif
mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
#define pb push_back
#define eb emplace_back
#define all(x) (x).begin(), (x).end()
#define f first
#define s second
#define g0(x) get<0>(x)
#define g1(x) get<1>(x)
#define g2(x) get<2>(x)
#define g3(x) get<3>(x)
typedef pair <int, int> pi;
typedef tuple<int,int,int> ti3;
typedef tuple<int,int,int,int> ti4;
int rand(int a, int b) { return a + rng() % (b-a+1); }
const int MOD = 1'000'002'022;
const int inf = (int)1e9 + 500;
const long long oo = (long long)1e18 + 500;
template <typename T> bool chmax(T& a, const T b) { return a<b ? a = b, 1 : 0; }
template <typename T> bool chmin(T& a, const T b) { return a>b ? a = b, 1 : 0; }
const int MAXN = 200005;
vector<int> V[MAXN];
int state[MAXN];
int C[MAXN];
int choices[MAXN];
int n,m;
void pre(int x) {
int S=V[x].size();
choices[x]=max(1ll,S);
for(auto v:V[x]) {
pre(v);
choices[x]=choices[x]*choices[v]%MOD;
}
}
void dfs(int x, int mult) {
int S=V[x].size();
if(S == 0) {
C[x-n]=mult;
return;
}
int pre[S+3];
int suf[S+3];
memset(pre, 0, sizeof pre);
memset(suf, 0, sizeof suf);
pre[0]=1;
suf[S+1]=1;
FOR(i,1,S) {
int node = V[x][i-1];
pre[i] = (pre[i-1] * choices[node]) % MOD;
}
DEC(i,S,1) {
int node = V[x][i-1];
suf[i] = (suf[i+1] * choices[node]) % MOD;
}
FOR(i,1,S) {
int node = V[x][i-1];
int mul = pre[i-1] * suf[i+1] % MOD;
dfs(node, mult * mul % MOD);
}
}
struct node {
int s,e,m,pot,lazy,sum;
node *l, *r;
node (int _s, int _e) {
s=_s;e=_e;m=(s+e)/2;
if(s!=e){
l=new node(s,m);
r=new node(m+1,e);
pot=(l->pot + r->pot)%MOD;
sum = lazy = 0;
} else {
pot=C[s];
}
}
void value() {
if(lazy==0)return;
sum=(MOD+(pot-sum)%MOD)%MOD;
if(s==e){
lazy=0;
return;
}
l->lazy^=lazy;
r->lazy^=lazy;
lazy=0;
}
void update(int x, int y) {
value();
if(s==x&&e==y) {
lazy=1;
return;
}
if(x>m)r->update(x,y);
else if(y<=m)l->update(x,y);
else l->update(x,m), r->update(m+1,y);
l->value(); r->value();
sum = (l->sum + r->sum)%MOD;
}
int query(int x, int y) {
value();
if(s==x&&e==y) return sum;
if(x>m)return r->query(x,y);
else if(y<=m) return l->query(x,y);
else return (l->query(x,m) + r->query(m+1,y))%MOD;
}
} *seg;
void init(int32_t n, int32_t m, vector<int32_t> P, vector<int32_t> A) {
::n=n;
::m=m;
FOR(i,1,n+m-1) {
V[P[i]].pb(i);
}
FOR(i,0,m-1) {
state[i]=A[i];
}
pre(0);
dfs(0,1);
seg=new node(0,m);
dba(C,0,m-1);
FOR(i,0,m-1) {
if(state[i])seg->update(i,i);
}
}
int32_t count_ways(int32_t L, int32_t R) {
L-=n;R-=n;
seg->update(L,R);
return seg->sum;
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
3 ms |
4944 KB |
Output is correct |
2 |
Correct |
3 ms |
4944 KB |
Output is correct |
3 |
Incorrect |
3 ms |
5200 KB |
1st lines differ - on the 1st token, expected: '509', found: '3102' |
4 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
3 ms |
4968 KB |
Output is correct |
2 |
Correct |
3 ms |
5072 KB |
Output is correct |
3 |
Correct |
3 ms |
5072 KB |
Output is correct |
4 |
Correct |
3 ms |
5072 KB |
Output is correct |
5 |
Correct |
3 ms |
5072 KB |
Output is correct |
6 |
Correct |
3 ms |
5072 KB |
Output is correct |
7 |
Correct |
3 ms |
5200 KB |
Output is correct |
8 |
Correct |
3 ms |
5200 KB |
Output is correct |
9 |
Correct |
3 ms |
5200 KB |
Output is correct |
10 |
Correct |
3 ms |
5456 KB |
Output is correct |
11 |
Correct |
4 ms |
5456 KB |
Output is correct |
12 |
Correct |
3 ms |
5220 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
3 ms |
4944 KB |
Output is correct |
2 |
Correct |
3 ms |
4944 KB |
Output is correct |
3 |
Incorrect |
3 ms |
5200 KB |
1st lines differ - on the 1st token, expected: '509', found: '3102' |
4 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
580 ms |
12948 KB |
Output is correct |
2 |
Correct |
974 ms |
20852 KB |
Output is correct |
3 |
Correct |
890 ms |
20864 KB |
Output is correct |
4 |
Correct |
834 ms |
20864 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
580 ms |
12948 KB |
Output is correct |
2 |
Correct |
974 ms |
20852 KB |
Output is correct |
3 |
Correct |
890 ms |
20864 KB |
Output is correct |
4 |
Correct |
834 ms |
20864 KB |
Output is correct |
5 |
Correct |
759 ms |
12944 KB |
Output is correct |
6 |
Correct |
1038 ms |
20820 KB |
Output is correct |
7 |
Correct |
981 ms |
20880 KB |
Output is correct |
8 |
Correct |
973 ms |
20868 KB |
Output is correct |
9 |
Correct |
336 ms |
5456 KB |
Output is correct |
10 |
Correct |
935 ms |
5968 KB |
Output is correct |
11 |
Correct |
816 ms |
5968 KB |
Output is correct |
12 |
Correct |
897 ms |
5968 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
3 ms |
4968 KB |
Output is correct |
2 |
Correct |
3 ms |
5072 KB |
Output is correct |
3 |
Correct |
3 ms |
5072 KB |
Output is correct |
4 |
Correct |
3 ms |
5072 KB |
Output is correct |
5 |
Correct |
3 ms |
5072 KB |
Output is correct |
6 |
Correct |
3 ms |
5072 KB |
Output is correct |
7 |
Correct |
3 ms |
5200 KB |
Output is correct |
8 |
Correct |
3 ms |
5200 KB |
Output is correct |
9 |
Correct |
3 ms |
5200 KB |
Output is correct |
10 |
Correct |
3 ms |
5456 KB |
Output is correct |
11 |
Correct |
4 ms |
5456 KB |
Output is correct |
12 |
Correct |
3 ms |
5220 KB |
Output is correct |
13 |
Correct |
580 ms |
12948 KB |
Output is correct |
14 |
Correct |
974 ms |
20852 KB |
Output is correct |
15 |
Correct |
890 ms |
20864 KB |
Output is correct |
16 |
Correct |
834 ms |
20864 KB |
Output is correct |
17 |
Correct |
759 ms |
12944 KB |
Output is correct |
18 |
Correct |
1038 ms |
20820 KB |
Output is correct |
19 |
Correct |
981 ms |
20880 KB |
Output is correct |
20 |
Correct |
973 ms |
20868 KB |
Output is correct |
21 |
Correct |
336 ms |
5456 KB |
Output is correct |
22 |
Correct |
935 ms |
5968 KB |
Output is correct |
23 |
Correct |
816 ms |
5968 KB |
Output is correct |
24 |
Correct |
897 ms |
5968 KB |
Output is correct |
25 |
Correct |
1054 ms |
28744 KB |
Output is correct |
26 |
Correct |
1158 ms |
29276 KB |
Output is correct |
27 |
Correct |
1146 ms |
29212 KB |
Output is correct |
28 |
Correct |
913 ms |
29232 KB |
Output is correct |
29 |
Correct |
1277 ms |
51068 KB |
Output is correct |
30 |
Correct |
854 ms |
51144 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
3 ms |
4944 KB |
Output is correct |
2 |
Correct |
3 ms |
4944 KB |
Output is correct |
3 |
Incorrect |
3 ms |
5200 KB |
1st lines differ - on the 1st token, expected: '509', found: '3102' |
4 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
3 ms |
4944 KB |
Output is correct |
2 |
Correct |
3 ms |
4944 KB |
Output is correct |
3 |
Incorrect |
3 ms |
5200 KB |
1st lines differ - on the 1st token, expected: '509', found: '3102' |
4 |
Halted |
0 ms |
0 KB |
- |