标签:str lin 实现 tor mission new describe return val
Time Limit: 8000/4000 MS (Java/Others) Memory Limit: 524288/524288 K (Java/Others)
Total Submission(s): 722 Accepted Submission(s): 255
#include <iostream> #include <vector> #define rep(i,a,b) for(int i=a;i<b;i++) #define per(i,a,b) for(int i=a-1;i>=b;i--) const int MX = 1e5; const int MXX = 400*MX; using namespace std; int n; vector<int> G[MX+5],vv[MX+5]; // init函数 实现将i因数分解 (i < 1e5) void init() { rep(i,1,MX+1) vv[i].push_back(1); rep(i,2,MX+1) { vv[i].push_back(i); for(int j=i+i;j<=MX;j+=i) vv[j].push_back(i); } // rep(i,1,MX+1) { // rep(j,0,vv[i].size()) { // printf("%d ",vv[i][j]); // }puts(""); // } } int root[MX+5],ls[MXX],rs[MXX],sum[MXX],rear,ans[MX]; inline void push_up(int rt) { if(ls[rt] && rs[rt]) sum[rt] = max(sum[ls[rt]],sum[rs[rt]]); else if(ls[rt]) sum[rt] = sum[ls[rt]]; else if(rs[rt]) sum[rt] = sum[rs[rt]]; } void update(int &rt,int l,int r,int p) { if(rt==0) rt = ++rear; if(l == r) { sum[rt] = p; return ; } int m = (l+r)>>1; if(p <= m) update(ls[rt],l,m,p); else update(rs[rt],m+1,r,p); push_up(rt); } int merge(int rt, int prt, int &ans) { if(rt==0 || prt==0) return rt^prt; //这里维护最大的gcd if(sum[rt] == sum[prt]) ans = max(ans,sum[rt]); //这里只有有因子,就归并到rt上面 if(ls[rt] | ls[prt]) ls[rt] = merge(ls[rt], ls[prt], ans); if(rs[rt] | rs[prt]) rs[rt] = merge(rs[rt], rs[prt], ans); push_up(rt); return rt; } void dfs(int u) { ans[u] = -1; rep(i, 0, G[u].size()) { int v = G[u][i]; dfs(v); root[u] = merge(root[u],root[v],ans[u]); } } int main () { freopen("in.txt" ,"r",stdin); freopen("out.txt","w",stdout); init(); scanf("%d", &n); //建边 rep(i,2,n+1) { int fa; scanf("%d",&fa); G[fa].push_back(i); } //对每个v[i]建线段树 rear=0; rep(i,1,n+1) { int x; scanf("%d", &x); root[i]=0; rep(j, 0, vv[x].size()) { update(root[i], 1, MX, vv[x][j]); } } //暴力更新 gcd dfs(1); //输出答案 rep(i,1,n+1) printf("%d\n", ans[i]); return 0; }
标签:str lin 实现 tor mission new describe return val
原文地址:https://www.cnblogs.com/Draymonder/p/9524241.html