30天网站建设实录视频,素材羊设计师服务平台,外链百科,做学术用的网站首先连通块#xff0c;所以点分治肯定是 Trick1 钦定选根的连通块dp 对于钦定选根的连通块dp#xff0c;有一种常见思路 先对原树求其dfn序#xff0c;按dfn序倒序求解 具体的#xff0c;对于当前点 i i i#xff08;注意这里都是指dfn序#xff09;#xff0c;我们…首先连通块所以点分治肯定是 Trick1 钦定选根的连通块dp 对于钦定选根的连通块dp有一种常见思路 先对原树求其dfn序按dfn序倒序求解 具体的对于当前点 i i i注意这里都是指dfn序我们可以钦定 i i i 是否选 如果 i i i 选就由 i 1 i1 i1也就是 i i i 的第一个儿子转移过来因为只有他选他子树才可能被选 如果 i i i 不选就由 i w i iw_i iwi​ 转移过来因为他的儿子必然不会被选 至于 i i i 和 i w i iw_i iwi​ 同时选的情况我们在 i 1 i1 i1 那里已经算了 对于 i i i 和 i w i iw_i iwi​ 是否连通的问题当他们的lca都被选时则他们必然也被选这里一定会在他们祖先那里被算到 Trick 2 对于乘积类dp的根号优化方法 考虑直接 d p [ x ] [ i ] dp[x][i] dp[x][i] i i i 值域过大。 但我们可以拆分 f ( x , i ) , g ( x , i ) f(x,i),g(x,i) f(x,i),g(x,i)代表已选乘积为 i i i / 还可以选乘积为 i i i 的方案数 这样状态直接压成 O ( m ) O(\sqrt m) O(m ​) 其实也可以用整除分块的证明进行预处理 #includebits/stdc.h using namespace std; #define int long long inline int read(){int x0,f1;char chgetchar(); while(ch0|| ch9){if(ch-)f-1;chgetchar();}while(ch0ch9){ x(x1)(x3)(ch^48);chgetchar();}return x*f;} #define Z(x) (x)*(x) #define pb push_back //mt19937 rand(time(0)); //mt19937_64 rand(time(0)); //srand(time(0)); #define N 4010 #define M 1510 #define mo (int)(1e97) int n, m, i, j, k, T; int Rt, rt, f[N][M], g[N][M]; int mx[N], w[N], dfn[N], tot, sum, u, v; int sq, p[N], v1, v2, a[N], ans; vectorintG[N]; void dfs(int x, int fa) {w[x]mx[x]1; for(int y : G[x]) {if(yfa || p[y]) continue; dfs(y, x); w[x]w[y]; mx[x]max(mx[x], w[y]); }mx[x]max(mx[x], sum-w[x]); if(mx[x]mx[rt]) rtx; }void dfs2(int x, int fa) {dfn[tot]x; for(int y: G[x]) if(y!fa !p[y]) dfs2(y, x); }void Add(int a, int b) {a(ab)%mo; }void dfz(int x) { // printf( %lld\n, x); int i, j, u; tot0; dfs(x, 0); dfs2(x, 0); // for(i1; itot; i) printf(%lld , dfn[i]); printf(\n); for(i0; itot5; i)for(j0; jsq5; j) f[i][j]g[i][j]0; // f[tot1][1]1; for(itot; i1; --i) {udfn[i]; if(a[u]sq) Add(g[i][m/a[u]], 1); else Add(f[i][a[u]], 1); for(j1; jsq; j) {v1i1; v2iw[u]; if(j*a[u]sq j*a[u]m) Add(g[i][m/(j*a[u])], f[v1][j]); else if(j*a[u]m) Add(f[i][j*a[u]], f[v1][j]); if(ja[u]) Add(g[i][j/a[u]], g[v1][j]); // // // Add(f[i][j], f[v2][j]); Add(g[i][j], g[v2][j]); }}for(i1; isq; i) Add(ans, f[1][i]g[1][i]); // printf(# %lld : %lld\n, x, ans); dfs(x, 0); p[x]1; for(int y : G[x]) if(!p[y]) {dfs(y, x); sumw[y]; mx[rt0]1e9; dfs(y, x); dfz(rt); } }signed main() { // freopen(in.txt, r, stdin); // freopen(out.txt, w, stdout);freopen(fn.in, r, stdin);freopen(fn.out, w, stdout); // Tread(); // while(T--) { // // }nread(); mread(); sqsqrt(m); // printf(# %lld\n, sq); for(i1; in; i) a[i]read(); for(i1; in; i) {uread(); vread(); G[u].pb(v); G[v].pb(u); }sumn; mx[rt0]1e9; dfs(1, 0); Rtrt; // printf(%lld\n, rt); dfz(rt);printf(%lld, (ans%momo)%mo); return 0; }