题解:
这道题用二次剩余或者矩阵都可以做。
相当于链剖之后维护一个区间加等比数列的可持久化线段树。 因为比相同,所以很好维护。
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const int RLEN=1<<18|1;
inline char nc() {
static char ibuf[RLEN],*ib,*ob;
(ib==ob) && (ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
return (ib==ob) ? -1 : *ib++;
}
inline int rd() {
char ch=nc(); int i=0,f=1;
while(!isdigit(ch)) {
if(ch=='-')f=-1; ch=nc();}
while(isdigit(ch)) {
i=(i<<1)+(i<<3)+ch-'0'; ch=nc();}
return i*f;
}
const int N=2e5+50, mod=1e9+9, base=383008016;
inline int add(int x,int y) {
return (x+y>=mod) ? (x+y-mod) : (x+y);}
inline int dec(int x,int y) {
return (x-y<0) ? (x-y+mod) : (x-y);}
inline int mul(int x,int y) {
return (LL)x*y%mod;}
inline int power(int a,int b,int rs=1) {
for(;b;b>>=1,a=mul(a,a)) if(b&1) rs=mul(rs,a); return rs;}
inline int inv(int a) {
return power(a,mod-2);}
const int base1=mul(base+1,inv(2)), base2=mul(dec(1,base),inv(2)), iv1=inv(base1), iv2=inv(base2), iv=inv(base);
int n,m,ind,lastans,pw1[N],pw2[N];
int fa[N],dep[N],sze[N],son[N],top[N],dfn[N];
vector <int> edge[N];
struct ST {
int bs,s[N],pw[N];
int rt[N],lc[N*100],rc[N*100],sum[N*100],d[N*100],tot;
inline void build(int &k,int l,int r) {
k=++tot; if(