Solution

计算每个质因子在哪些点出现，然后在树上只保存这些点，通过计算大于000的g(i,j)g(i,j)g(i,j)个数来计算贡献
如果质因子次数>1>1>1，那么把p2p^2p2，p3p^3p3，…，pkp^kpk都按同样方法做一遍
复杂度O((n+ai)logai)O((n+a_i)loga_i)O((n+ai)logai)

Code

代码中*q就相当于q[0]q[0]q[0]，*vis，*vis1同理

#include<bits/stdc++.h>
using namespace std;
#define fi first
#define se second
typedef pair<int,int> pi;
const int N=100002,M=1e9+7,P=200002;
struct node{int to,ne;
}e[P<<1];
int a[P],b[P],cnt,tot,ans,vis[P],vis1[P],q[P],h[P],n,i,x,y,num,j,p,t,k;
vector<int>d[N],v[N][17];
bool fl[N];
inline char gc(){static char buf[100000],*p1=buf,*p2=buf;return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int rd(){int x=0,fl=1;char ch=gc();for(;ch<48||ch>57;ch=gc())if(ch=='-')fl=-1;for(;48<=ch&&ch<=57;ch=gc())x=(x<<3)+(x<<1)+(ch^48);return x*fl;
}
pi dfs(int u,int fa){vis1[u]=*vis1;pi res(b[u]>0,b[u]==0);for (int i=h[u],v;i;i=e[i].ne)if ((v=e[i].to)!=fa && vis[v]==*vis){pi sub=dfs(v,u);if (b[u]) sub.fi+=sub.se,sub.se=0;tot+=res.fi*sub.fi+res.fi*sub.se+res.se*sub.fi;res.fi+=sub.fi,res.se+=sub.se;}return res;
}
int pw(int x,int y){int z=1;for (;y;y>>=1,x=1ll*x*x%M)if (y&1) z=1ll*z*x%M;return z;
}
void add(int x,int y){e[++tot]=(node){y,h[x]},h[x]=tot;}
int main(){n=rd();for (i=1;i<n;i++) x=rd(),y=rd(),add(x,y),add(y,x);for (i=1;i<=n;i++) a[i]=rd();for (i=1;i<=n;i++) b[i]=rd();for (i=2;i<N;i++)if (!fl[i]){d[i].push_back(i);//i也是i的约数，不能忘记for (j=i<<1;j<N;j+=i) fl[j]=1,d[j].push_back(i);}ans=1;for (i=1;i<=n;i++){if (b[i]) ans=1ll*ans*a[i]%M;for (j=0;j<d[a[i]].size();j++){p=d[a[i]][j];for (t=a[i],cnt=0;t%p==0;t/=p,cnt++) v[p][cnt].push_back(i);}}for (i=2;i<N;i++) if (!fl[i]){p=i;for (t=N-1,cnt=0;t;t/=p,cnt++);tot=0;for (j=0;j<cnt-1;j++) if (v[i][j].size()){++*vis,++*vis1,*q=0;for (k=0;k<v[i][j].size();k++) vis[q[++*q]=v[i][j][k]]=*vis;for (k=1;k<=*q;k++)if (vis1[q[k]]!=*vis1) dfs(q[k],0);}ans=1ll*ans*pw(i,tot)%M;}printf("%d",ans);
}

一中模拟赛3.15——树上gcd相关推荐

2018.7.17 绍兴一中模拟赛解题报告
成绩: 100+100+50=250(考得不错,只可惜跑得太慢,是250中最慢的一个) T1:交换(swap) 满分解法1: 这题是真的水. 先介绍一个比较麻烦的做法,也是我考试时的做法. 大致思路就 ...
[wikioi 1418]铃仙•优昙华院稻叶（东方幻想乡系列模拟赛）（树上递推）
题目:http://www.wikioi.com/problem/1418/ 分析: 一看就肯定是树上的递推设f[i][j][k]表示第i秒在k点(从j点走过来的)的概率则f[i][j][k]=f ...
数字对（长乐一中模拟赛day2T2）
2.数字对 [题目描述] 小H是个善于思考的学生,现在她又在思考一个有关序列的问题. 她的面前浮现出一个长度为n的序列{ai},她想找出一段区间[L, R](1 <= L <= R < ...
改造二叉树 (长乐一中模拟赛day2T1)
1.改造二叉树 [题目描述] 小Y在学树论时看到了有关二叉树的介绍:在计算机科学中,二叉树是每个结点最多有两个子结点的有序树.通常子结点被称作"左孩子"和"右孩子&quo ...
绍兴一中模拟赛3.19——时光流转
Description Solution 离线以后点分对于每个点,都用这个点的祖先把这个点的子树更新一遍, 考虑到操作时间早的才能更新晚的和题目中说的"路径上边权都大于等于 v a l v ...
2020.03.11模拟赛15（第一题）
1.水果盛宴(fruit) 题目描述贝茜又再一次地闯入了 Farmer John 的房子!她在厨房发现了一堆柠檬和一堆橘子(每堆都有无限多个),并且,她希望尽可能地多吃. 贝茜的有一个饱腹值上限 T ...
【模拟赛】10.15—10.17
背景: 继国庆自闭四联考之后,确实,是很自闭,还做了些比较疯狂的事情,然后顿悟了很多,然后刚知道怎么去努力,刚知道努力的方法和方向就赶快滚回来写博客了,毕竟是在将近考试的时候出现了问题,还是赶快狗回来 ...
省选模拟赛记录(越往下越新哦~~~)
LOG 模拟赛第一次见尼玛这么给数据范围的-- 开考有点困,迷迷糊糊看完了三道题,真的是像老吕说的那样,一道都不会-- 思考T1,感觉有点感觉,但是太困了,就先码了暴力,发现打表可以50分,于是就大 ...
第十三届蓝桥杯模拟赛（第三期）试题与题解 C++
文章目录第十三届蓝桥杯模拟赛(第三期)试题与题解 1.试题A 题解:数制转换 2.试题B 题解:枚举 3.试题C 题解:枚举 4.试题D 题解:最小生成树 5.试题E 方法一:暴力求和方法二:一维 ...

一中模拟赛3.15——树上gcd

Solution

Code

一中模拟赛3.15——树上gcd相关推荐

最新文章

热门文章