#include "circuit.h"
#include<bits/stdc++.h>
#include <vector>
using namespace std;
#define pb push_back
#define fi first
#define se second
#define ll long long
const int N=1e5+500;
const ll mod=1000002022;
ll Puta(ll a,ll b){a%=mod,b%=mod;ll x=a*b;x%=mod;return x;}
ll Plus(ll a,ll b){a%=mod,b%=mod;ll x=a+b;x%=mod;while(x<0) x+=mod;return x;}
vector<int>E[2*N+500];
int a[N+500],n,m,P[2*N+500];
ll res[N+500],dp[2*N+500][2],F[2*N+500],F1[2*N+500];
bool stepen2;
void DFSsetup(int u){
if(E[u].empty()) return;
F[u]=E[u].size();
for(auto i:E[u]){
DFSsetup(i);
F[u]=Puta(F[u],F[i]);
}
for(auto i:E[u]){
F1[i]=1;
for(auto j:E[u]){
if(j!=i) F1[i]=Puta(F1[i],F[j]);
}
}
}
void DFS(int u){
if(E[u].empty()) return;
for(auto i:E[u]) DFS(i);
dp[u][1]=0;
for(auto i:E[u]) dp[u][1]=Plus(dp[u][1],Puta(dp[i][1],F1[i]));
//printf("%i: %lld %lld\n",u,dp[u][1],F[u]);
/*int mask=0;
while(mask< (1<<(E[u].size()))){
ll temp1=1,sum=0;
for(int i=0,temp=mask;i<E[u].size();i++,temp/=2){
temp1=Puta(temp1,dp[E[u][i]][temp%2]);sum+=temp%2;
}
if(mask!=0) dp[u][1]=Plus(dp[u][1],Puta(sum,temp1));
dp[u][0]=Plus(dp[u][0],Puta(E[u].size()-sum,temp1));
mask++;
}*/
//dp[u][0]=F[u]-dp[u][1];
}
void Update1(int v){
dp[v][1]=1-dp[v][1];
dp[v][0]=1-dp[v][0];
int u=P[v];
while(u!=-1){
dp[u][1]=0;
for(auto i:E[u]) dp[u][1]=Plus(dp[u][1],Puta(dp[i][1],F1[i]));
/*dp[u][1]=dp[u][0]=0;
int mask=0;
while(mask< (1<<(E[u].size()))){
ll temp1=1,sum=0;
for(int i=0,temp=mask;i<E[u].size();i++,temp/=2){
temp1=Puta(temp1,dp[E[u][i]][temp%2]);sum+=temp%2;
}
if(mask!=0) dp[u][1]=Plus(dp[u][1],Puta(sum,temp1));
dp[u][0]=Plus(dp[u][0],Puta(E[u].size()-sum,temp1));
mask++;
}*/
u=P[u];
}
}
int nc,root,lc[2*N],rc[2*N],sum[2*N],lazy[2*N];
void update(int &c,int ss,int se,int qval){
if(!c) c=++nc;
if(qval==0) return;
sum[c]=se-ss+1-sum[c];
lazy[c]++;lazy[c]%=2;
}
void push(int &c,int ss,int se){
int mid=(ss+se)/2;
update(lc[c],ss,mid,lazy[c]);
update(rc[c],mid+1,se,lazy[c]);
lazy[c]=0;
}
void Update(int &c,int ss,int se,int qs,int qe){
if(!c) c=++nc;
if(qs<=ss && se<=qe){
update(c,ss,se,1);
//printf("%i %i %i: %i\n",c,ss,se,sum[c]);
return;
}
if(qe<ss || se<qs) return;
int mid=(ss+se)/2;
push(c,ss,se);
Update(lc[c],ss,mid,qs,qe);Update(rc[c],mid+1,se,qs,qe);
sum[c]=sum[lc[c]]+sum[rc[c]];
//printf("%i %i %i: %i\n",c,ss,se,sum[c]);
}
int Get(int c,int ss,int se,int qs,int qe){
//printf("%i %i %i: %i\n",c,ss,se,sum[c]);
if(qs<=ss && se<=qe) return sum[c];
if(qe<ss || se<qs) return 0;
int mid=(ss+se)/2;
push(c,ss,se);
return Get(lc[c],ss,mid,qs,qe)+Get(rc[c],mid+1,se,qs,qe);
}
int tajm,in[2*N],out[2*N],LL[2*N],RR[2*N];
void DFSsetup2(int u){
in[u]=++tajm;
for(auto i:E[u]) DFSsetup2(i);
out[u]=tajm;
}
int lazy2[2*N];
void update2(int u,int qval){
if(qval==0) return;
dp[u][1]=Plus(F[u],-dp[u][1]);
lazy2[u]++;lazy2[u]%=2;
}
void push2(int u){
for(auto i:E[u]) update2(i,lazy2[u]);
lazy2[u]=0;
}
void Update2(int u,int qs,int qe){
if(qs<=LL[u] && RR[u]<=qe) {update2(u,1);return;}
if(qe<LL[u] || RR[u]<qs) return;
push2(u);
for(auto i:E[u]) Update2(i,qs,qe);
dp[u][1]=0;
for(auto i:E[u]) dp[u][1]=Plus(dp[u][1],Puta(dp[i][1],F1[i]));
//printf("%i %i %i: %i\n",u,LL[u],RR[u],dp[u][1]);
}
void init(int N1, int M1, std::vector<int> P1, std::vector<int> A) {
n=N1,m=M1;
for(int i=0;i<P1.size();i++) P[i]=P1[i];
for(int i=0;i<P1.size();i++){
E[P[i]].pb(i);
}
for(int i=0;i<m;i++){
dp[n+i][1]=A[i],dp[n+i][0]=1-A[i];
a[i]=A[i];
F[n+i]=1;
}
DFSsetup2(0);
for(int i=0;i<n+m;i++){
int lb=lower_bound(in+n,in+n+m-1,in[i])-in,ub=lower_bound(out+n,out+n+m-1,out[i])-out;
LL[i]=lb,RR[i]=ub;
//printf("%i: %i %i %i %i\n",i,in[i],out[i],LL[i],RR[i]);
}
DFSsetup(0);
DFS(0);
int temp=m,ct=0;while(temp>0) ct+=temp%2,temp/=2;if(ct==1) stepen2=true;
if(!stepen2 || n<=5000) return;
for(int i=n;i<n+m;i++) dp[i][1]=0;
DFS(0);
res[0]=dp[0][1];
for(int i=0;i<m;i++){
Update1(n+i);
res[i+1]=dp[0][1];
}
for(int i=0;i<m;i++) if(a[i]==1) Update(root,0,m-1,i,i);
}
int count_ways(int L, int R) {
if(n<=5000){
for(int i=L;i<=R;i++) dp[i][1]=1-dp[i][1];
DFS(0);
//Update2(0,L,R);
return dp[0][1];
}
if(stepen2){
Update(root,0,m-1,L-n,R-n);
return res[Get(root,0,m-1,0,m-1)];
}
for(int i=L;i<=R;i++) Update1(i);
return dp[0][1];
}
Compilation message
circuit.cpp: In function 'void init(int, int, std::vector<int>, std::vector<int>)':
circuit.cpp:134:15: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
134 | for(int i=0;i<P1.size();i++) P[i]=P1[i];
| ~^~~~~~~~~~
circuit.cpp:135:15: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
135 | for(int i=0;i<P1.size();i++){
| ~^~~~~~~~~~
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
14168 KB |
Output is correct |
2 |
Correct |
3 ms |
14168 KB |
Output is correct |
3 |
Correct |
10 ms |
14244 KB |
Output is correct |
4 |
Correct |
11 ms |
14168 KB |
Output is correct |
5 |
Correct |
11 ms |
14240 KB |
Output is correct |
6 |
Correct |
10 ms |
14168 KB |
Output is correct |
7 |
Correct |
10 ms |
14164 KB |
Output is correct |
8 |
Correct |
11 ms |
14168 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
14168 KB |
Output is correct |
2 |
Correct |
3 ms |
14168 KB |
Output is correct |
3 |
Correct |
3 ms |
14168 KB |
Output is correct |
4 |
Correct |
4 ms |
14168 KB |
Output is correct |
5 |
Correct |
4 ms |
14004 KB |
Output is correct |
6 |
Correct |
3 ms |
14424 KB |
Output is correct |
7 |
Correct |
3 ms |
14168 KB |
Output is correct |
8 |
Correct |
3 ms |
14168 KB |
Output is correct |
9 |
Correct |
3 ms |
14168 KB |
Output is correct |
10 |
Correct |
3 ms |
14168 KB |
Output is correct |
11 |
Correct |
4 ms |
14168 KB |
Output is correct |
12 |
Correct |
3 ms |
14168 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
14168 KB |
Output is correct |
2 |
Correct |
3 ms |
14168 KB |
Output is correct |
3 |
Correct |
10 ms |
14244 KB |
Output is correct |
4 |
Correct |
11 ms |
14168 KB |
Output is correct |
5 |
Correct |
11 ms |
14240 KB |
Output is correct |
6 |
Correct |
10 ms |
14168 KB |
Output is correct |
7 |
Correct |
10 ms |
14164 KB |
Output is correct |
8 |
Correct |
11 ms |
14168 KB |
Output is correct |
9 |
Correct |
3 ms |
14168 KB |
Output is correct |
10 |
Correct |
3 ms |
14168 KB |
Output is correct |
11 |
Correct |
3 ms |
14168 KB |
Output is correct |
12 |
Correct |
4 ms |
14168 KB |
Output is correct |
13 |
Correct |
4 ms |
14004 KB |
Output is correct |
14 |
Correct |
3 ms |
14424 KB |
Output is correct |
15 |
Correct |
3 ms |
14168 KB |
Output is correct |
16 |
Correct |
3 ms |
14168 KB |
Output is correct |
17 |
Correct |
3 ms |
14168 KB |
Output is correct |
18 |
Correct |
3 ms |
14168 KB |
Output is correct |
19 |
Correct |
4 ms |
14168 KB |
Output is correct |
20 |
Correct |
3 ms |
14168 KB |
Output is correct |
21 |
Correct |
3 ms |
14168 KB |
Output is correct |
22 |
Correct |
3 ms |
14168 KB |
Output is correct |
23 |
Correct |
3 ms |
14168 KB |
Output is correct |
24 |
Correct |
3 ms |
14168 KB |
Output is correct |
25 |
Correct |
4 ms |
14168 KB |
Output is correct |
26 |
Correct |
3 ms |
14168 KB |
Output is correct |
27 |
Correct |
3 ms |
14168 KB |
Output is correct |
28 |
Correct |
3 ms |
14420 KB |
Output is correct |
29 |
Correct |
10 ms |
14168 KB |
Output is correct |
30 |
Correct |
10 ms |
14168 KB |
Output is correct |
31 |
Correct |
3 ms |
14168 KB |
Output is correct |
32 |
Correct |
4 ms |
14168 KB |
Output is correct |
33 |
Correct |
4 ms |
14168 KB |
Output is correct |
34 |
Correct |
3 ms |
14168 KB |
Output is correct |
35 |
Correct |
4 ms |
14168 KB |
Output is correct |
36 |
Correct |
4 ms |
14168 KB |
Output is correct |
37 |
Correct |
10 ms |
14368 KB |
Output is correct |
38 |
Correct |
10 ms |
14368 KB |
Output is correct |
39 |
Correct |
4 ms |
14168 KB |
Output is correct |
40 |
Correct |
3 ms |
14168 KB |
Output is correct |
41 |
Correct |
4 ms |
14188 KB |
Output is correct |
42 |
Correct |
3 ms |
14168 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
459 ms |
19480 KB |
Output is correct |
2 |
Correct |
671 ms |
22412 KB |
Output is correct |
3 |
Correct |
703 ms |
22684 KB |
Output is correct |
4 |
Correct |
658 ms |
22680 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
459 ms |
19480 KB |
Output is correct |
2 |
Correct |
671 ms |
22412 KB |
Output is correct |
3 |
Correct |
703 ms |
22684 KB |
Output is correct |
4 |
Correct |
658 ms |
22680 KB |
Output is correct |
5 |
Correct |
577 ms |
19496 KB |
Output is correct |
6 |
Correct |
708 ms |
22392 KB |
Output is correct |
7 |
Correct |
745 ms |
22480 KB |
Output is correct |
8 |
Correct |
656 ms |
22180 KB |
Output is correct |
9 |
Correct |
2256 ms |
14168 KB |
Output is correct |
10 |
Execution timed out |
3083 ms |
14424 KB |
Time limit exceeded |
11 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
14168 KB |
Output is correct |
2 |
Correct |
3 ms |
14168 KB |
Output is correct |
3 |
Correct |
3 ms |
14168 KB |
Output is correct |
4 |
Correct |
4 ms |
14168 KB |
Output is correct |
5 |
Correct |
4 ms |
14004 KB |
Output is correct |
6 |
Correct |
3 ms |
14424 KB |
Output is correct |
7 |
Correct |
3 ms |
14168 KB |
Output is correct |
8 |
Correct |
3 ms |
14168 KB |
Output is correct |
9 |
Correct |
3 ms |
14168 KB |
Output is correct |
10 |
Correct |
3 ms |
14168 KB |
Output is correct |
11 |
Correct |
4 ms |
14168 KB |
Output is correct |
12 |
Correct |
3 ms |
14168 KB |
Output is correct |
13 |
Correct |
459 ms |
19480 KB |
Output is correct |
14 |
Correct |
671 ms |
22412 KB |
Output is correct |
15 |
Correct |
703 ms |
22684 KB |
Output is correct |
16 |
Correct |
658 ms |
22680 KB |
Output is correct |
17 |
Correct |
577 ms |
19496 KB |
Output is correct |
18 |
Correct |
708 ms |
22392 KB |
Output is correct |
19 |
Correct |
745 ms |
22480 KB |
Output is correct |
20 |
Correct |
656 ms |
22180 KB |
Output is correct |
21 |
Correct |
2256 ms |
14168 KB |
Output is correct |
22 |
Execution timed out |
3083 ms |
14424 KB |
Time limit exceeded |
23 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
14168 KB |
Output is correct |
2 |
Correct |
3 ms |
14168 KB |
Output is correct |
3 |
Correct |
10 ms |
14244 KB |
Output is correct |
4 |
Correct |
11 ms |
14168 KB |
Output is correct |
5 |
Correct |
11 ms |
14240 KB |
Output is correct |
6 |
Correct |
10 ms |
14168 KB |
Output is correct |
7 |
Correct |
10 ms |
14164 KB |
Output is correct |
8 |
Correct |
11 ms |
14168 KB |
Output is correct |
9 |
Correct |
3 ms |
14168 KB |
Output is correct |
10 |
Correct |
3 ms |
14168 KB |
Output is correct |
11 |
Correct |
3 ms |
14168 KB |
Output is correct |
12 |
Correct |
4 ms |
14168 KB |
Output is correct |
13 |
Correct |
4 ms |
14004 KB |
Output is correct |
14 |
Correct |
3 ms |
14424 KB |
Output is correct |
15 |
Correct |
3 ms |
14168 KB |
Output is correct |
16 |
Correct |
3 ms |
14168 KB |
Output is correct |
17 |
Correct |
3 ms |
14168 KB |
Output is correct |
18 |
Correct |
3 ms |
14168 KB |
Output is correct |
19 |
Correct |
4 ms |
14168 KB |
Output is correct |
20 |
Correct |
3 ms |
14168 KB |
Output is correct |
21 |
Correct |
3 ms |
14168 KB |
Output is correct |
22 |
Correct |
3 ms |
14168 KB |
Output is correct |
23 |
Correct |
3 ms |
14168 KB |
Output is correct |
24 |
Correct |
3 ms |
14168 KB |
Output is correct |
25 |
Correct |
4 ms |
14168 KB |
Output is correct |
26 |
Correct |
3 ms |
14168 KB |
Output is correct |
27 |
Correct |
3 ms |
14168 KB |
Output is correct |
28 |
Correct |
3 ms |
14420 KB |
Output is correct |
29 |
Correct |
10 ms |
14168 KB |
Output is correct |
30 |
Correct |
10 ms |
14168 KB |
Output is correct |
31 |
Correct |
3 ms |
14168 KB |
Output is correct |
32 |
Correct |
4 ms |
14168 KB |
Output is correct |
33 |
Correct |
4 ms |
14168 KB |
Output is correct |
34 |
Correct |
3 ms |
14168 KB |
Output is correct |
35 |
Correct |
4 ms |
14168 KB |
Output is correct |
36 |
Correct |
4 ms |
14168 KB |
Output is correct |
37 |
Correct |
10 ms |
14368 KB |
Output is correct |
38 |
Correct |
10 ms |
14368 KB |
Output is correct |
39 |
Correct |
4 ms |
14168 KB |
Output is correct |
40 |
Correct |
3 ms |
14168 KB |
Output is correct |
41 |
Correct |
4 ms |
14188 KB |
Output is correct |
42 |
Correct |
3 ms |
14168 KB |
Output is correct |
43 |
Execution timed out |
3035 ms |
14168 KB |
Time limit exceeded |
44 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
14168 KB |
Output is correct |
2 |
Correct |
3 ms |
14168 KB |
Output is correct |
3 |
Correct |
10 ms |
14244 KB |
Output is correct |
4 |
Correct |
11 ms |
14168 KB |
Output is correct |
5 |
Correct |
11 ms |
14240 KB |
Output is correct |
6 |
Correct |
10 ms |
14168 KB |
Output is correct |
7 |
Correct |
10 ms |
14164 KB |
Output is correct |
8 |
Correct |
11 ms |
14168 KB |
Output is correct |
9 |
Correct |
3 ms |
14168 KB |
Output is correct |
10 |
Correct |
3 ms |
14168 KB |
Output is correct |
11 |
Correct |
3 ms |
14168 KB |
Output is correct |
12 |
Correct |
4 ms |
14168 KB |
Output is correct |
13 |
Correct |
4 ms |
14004 KB |
Output is correct |
14 |
Correct |
3 ms |
14424 KB |
Output is correct |
15 |
Correct |
3 ms |
14168 KB |
Output is correct |
16 |
Correct |
3 ms |
14168 KB |
Output is correct |
17 |
Correct |
3 ms |
14168 KB |
Output is correct |
18 |
Correct |
3 ms |
14168 KB |
Output is correct |
19 |
Correct |
4 ms |
14168 KB |
Output is correct |
20 |
Correct |
3 ms |
14168 KB |
Output is correct |
21 |
Correct |
3 ms |
14168 KB |
Output is correct |
22 |
Correct |
3 ms |
14168 KB |
Output is correct |
23 |
Correct |
3 ms |
14168 KB |
Output is correct |
24 |
Correct |
3 ms |
14168 KB |
Output is correct |
25 |
Correct |
4 ms |
14168 KB |
Output is correct |
26 |
Correct |
3 ms |
14168 KB |
Output is correct |
27 |
Correct |
3 ms |
14168 KB |
Output is correct |
28 |
Correct |
3 ms |
14420 KB |
Output is correct |
29 |
Correct |
10 ms |
14168 KB |
Output is correct |
30 |
Correct |
10 ms |
14168 KB |
Output is correct |
31 |
Correct |
3 ms |
14168 KB |
Output is correct |
32 |
Correct |
4 ms |
14168 KB |
Output is correct |
33 |
Correct |
4 ms |
14168 KB |
Output is correct |
34 |
Correct |
3 ms |
14168 KB |
Output is correct |
35 |
Correct |
4 ms |
14168 KB |
Output is correct |
36 |
Correct |
4 ms |
14168 KB |
Output is correct |
37 |
Correct |
10 ms |
14368 KB |
Output is correct |
38 |
Correct |
10 ms |
14368 KB |
Output is correct |
39 |
Correct |
4 ms |
14168 KB |
Output is correct |
40 |
Correct |
3 ms |
14168 KB |
Output is correct |
41 |
Correct |
4 ms |
14188 KB |
Output is correct |
42 |
Correct |
3 ms |
14168 KB |
Output is correct |
43 |
Correct |
459 ms |
19480 KB |
Output is correct |
44 |
Correct |
671 ms |
22412 KB |
Output is correct |
45 |
Correct |
703 ms |
22684 KB |
Output is correct |
46 |
Correct |
658 ms |
22680 KB |
Output is correct |
47 |
Correct |
577 ms |
19496 KB |
Output is correct |
48 |
Correct |
708 ms |
22392 KB |
Output is correct |
49 |
Correct |
745 ms |
22480 KB |
Output is correct |
50 |
Correct |
656 ms |
22180 KB |
Output is correct |
51 |
Correct |
2256 ms |
14168 KB |
Output is correct |
52 |
Execution timed out |
3083 ms |
14424 KB |
Time limit exceeded |
53 |
Halted |
0 ms |
0 KB |
- |