BZOJ4564 [Haoi2016]地图

仙人掌转成序列，然后莫队

仙人掌转成序列的方法是先tarjan找环，对于一个环上的点，如果是环根，那么这个环都属于其子仙人掌，如果不是环根，那么除了他自己这个环上都不是其子仙人掌，这样再dfs一次，对于一个点，优先走其出边里是其子仙人掌的部分，同时对于一个环，把不与环根相连的点的siz加到环根的siz里，这样仙人掌就转成了序列，一个子仙人掌对应的区间就是dfn[x]~dfn[x]+siz[x]-1

之后直接上莫队，用数据结构维护的话是n sqrtn logn的，T到爆，用和GTY的二逼妹子序列相同的方法，对权值分块，维护每块和每个点对答案的贡献，修改是单点修改，这样修改其所在块和这个位置，询问的时候变成sqrtn即可

数据好水啊，我在第二次dfs判这条出边和x是不是在一个环的时候没用tarjan的dfn而用了新dfs序都过了-_-

人生中第一次在BZOJ首A-_-

update：NOI拍照现场看到松爷，吓得我赶紧复习一发仙人掌，结果发现之前的代码有点问题……不过这题数据太弱了，写成啥JB样都能过

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<bitset>
using namespace std;
#define MAXN 100010
#define MAXM 300010
#define MAXC 1000010
#define ll long long
#define INF 1000000000
#define MOD 1000000007
#define eps 1e-8
struct vec{
int to;
int fro;
};
struct que{
int cl;
int v;
int l;
int r;
int k;
int num;
friend bool operator <(que x,que y){
return x.k!=y.k?x.k<y.k:x.r<y.r;
}
};
int A[MAXN],a[MAXN];
int c;
que q[MAXN];
vec mp[MAXM];
int tai[MAXN],cnt=1;
int dfn[MAXN],DFN[MAXN],low[MAXN],ndf[MAXN],tim;
int n,m;
int siz[MAXN];
int ANS[2][MAXC];
int CNT[MAXC];
int Asiz;
int ans[MAXN];
inline void be(int x,int y){
mp[++cnt].to=y;
mp[cnt].fro=tai[x];
tai[x]=cnt;
}
inline void bde(int x,int y){
be(x,y);
be(y,x);
}
void dfs(int x,int f){
int i,y;
DFN[x]=low[x]=++tim;
ndf[tim]=x;
for(i=tai[x];i;i=mp[i].fro){
y=mp[i].to;
if(i==(f^1)){
continue ;
}
if(!DFN[y]){
dfs(y,i);
low[x]=min(low[x],low[y]);
}else{
low[x]=min(low[x],DFN[y]);
}
}
}
void dfs2(int x,int f){
int i,y;
siz[x]=1;
dfn[x]=++tim;
for(i=tai[x];i;i=mp[i].fro){
y=mp[i].to;
if(i==(f^1)){
continue ;
}
if(!dfn[y]&&low[y]>=DFN[x]){
dfs2(y,i);
if(low[y]>DFN[x]){
siz[x]+=siz[y];
}
}
}
for(i=tai[x];i;i=mp[i].fro){
y=mp[i].to;
if(i==(f^1)){
continue ;
}
if(!dfn[y]&&low[y]<DFN[x]){
dfs2(y,i);
}
}
if(low[x]!=DFN[x]){
siz[ndf[low[x]]]+=siz[x];
}
}
int main(){
int i,j,o,x,y;
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++){
scanf("%d",&A[i]);
c=max(c,A[i]);
}
Asiz=sqrt(c);
for(i=1;i<=m;i++){
scanf("%d%d",&x,&y);
bde(x,y);
}
dfs(1,0);
tim=0;
dfs2(1,0);
for(i=1;i<=n;i++){
a[dfn[i]]=A[i];
}
scanf("%d",&m);
int SIZ=sqrt(n);
for(i=1;i<=m;i++){
scanf("%d%d%d",&q[i].cl,&q[i].l,&q[i].v);
q[i].r=dfn[q[i].l]+siz[q[i].l]-1;
q[i].l=dfn[q[i].l];
q[i].k=q[i].l/SIZ;
q[i].num=i;
}
sort(q+1,q+m+1);
int L=1,R=0;
for(i=1;i<=m;i++){
while(L>q[i].l){
L--;
x=(a[L]-1)/Asiz+1;
if(CNT[a[L]]){
ANS[CNT[a[L]]&1][x]--;
}
CNT[a[L]]++;
ANS[CNT[a[L]]&1][x]++;
}
while(R<q[i].r){
R++;
x=(a[R]-1)/Asiz+1;
if(CNT[a[R]]){
ANS[CNT[a[R]]&1][x]--;
}
CNT[a[R]]++;
ANS[CNT[a[R]]&1][x]++;
}
while(L<q[i].l){
x=(a[L]-1)/Asiz+1;
ANS[CNT[a[L]]&1][x]--;
CNT[a[L]]--;
if(CNT[a[L]]){
ANS[CNT[a[L]]&1][x]++;
}
L++;
}
while(R>q[i].r){
x=(a[R]-1)/Asiz+1;
ANS[CNT[a[R]]&1][x]--;
CNT[a[R]]--;
if(CNT[a[R]]){
ANS[CNT[a[R]]&1][x]++;
}
R--;
}
x=(q[i].v-1)/Asiz+1;
for(j=1;j<x;j++){
ans[q[i].num]+=ANS[q[i].cl][j];
}
for(j=(x-1)*Asiz+1;j<=q[i].v;j++){
if(CNT[j]&&((CNT[j]&1)==q[i].cl)){
ans[q[i].num]++;
}
}
}
for(i=1;i<=m;i++){
printf("%d\n",ans[i]);
}
return 0;
}

/*
7 8
1 2
2 3
3 4
4 1
3 5
5 6
6 7
7 3

5 6
2 1 6 7 7
1 2
1 3
2 4
4 5
4 5
1 3
3
0 3 2
0 3 1
0 1 7

*/