题意：

给你一颗树，其中根是111，每个点有一个点权，求每个点到根的所有路径的gcdgcdgcd之和。
n≤1e5n\le1e5n≤1e5

思路：

一看到以为是个点分治，让后发现不是任意点对所以显然不能用点分治来做。
后来想了半天也只想到一个暴力，没想到还真是个暴力。
考虑从根到点iii的路径上同的gcdgcdgcd个数，可以证明最多有log2(1012)log_2(10^{12})log2(1012)个，所以每个点都继承一下父节点的值即可，由于需要用mapmapmap维护，所以复杂度O(nlog22(1012))O(nlog^2_2(10^{12}))O(nlog22(1012))。

// Problem: E. Kamil and Making a Stream
// Contest: Codeforces - Codeforces Round #588 (Div. 2)
// URL: https://codeforces.com/contest/1230/problem/E
// Memory Limit: 768 MB
// Time Limit: 4000 ms
//
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#include<random>
#include<cassert>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid ((tr[u].l+tr[u].r)>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=1000010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;int n;
vector<int>v[N];
vector<pair<LL,LL>>all[N];
LL a[N],ans;
map<LL,int>mp;void del() {if(ans>=mod) ans-=mod;
}void dfs(int u,int f) {mp.clear();ans+=a[u]; del();for(auto y:all[f]) {ans+=y.Y*__gcd(y.X,a[u])%mod; ans%=mod;mp[__gcd(y.X,a[u])]+=y.Y;}mp[a[u]]++;for(auto x:mp) all[u].pb({x.X,x.Y});for(auto x:v[u]) {if(x==f) continue;dfs(x,u);}
}int main()
{//  ios::sync_with_stdio(false);
//  cin.tie(0);cin>>n;for(int i=1;i<=n;i++) scanf("%lld",&a[i]);for(int i=1;i<=n-1;i++) {int a,b; scanf("%d%d",&a,&b);v[a].pb(b); v[b].pb(a);}dfs(1,0);printf("%lld\n",ans%mod);return 0;
}
/**/

Codeforces Round #588 (Div. 2) E. Kamil and Making a Stream 数学 + 暴力相关推荐

Codeforces Round #588 (Div. 2) F. Konrad and Company Evaluation 图论 + 建反图好题
传送门文章目录题意: 思路: 题意: 给你一张nnn个点mmm条边的图,其中每个点iii初始编号为iii,边是有向的,方向为从编号大的指向编号小的.定义一个贡献为存在某三个点a,b,ca,b,ca ...
Codeforces Round #588 (Div. 2) D. Marcin and Training Camp 思维
传送门文章目录题意: 思路: 题意: 有nnn个人,每个人都有一个能力值bib_ibi以及他会的技能aia_iai,当他会第xxx个技能的时候,aia_iai的第xxx位是111.定义当xx ...
Codeforces Round #369 (Div. 2) B. Chris and Magic Square【数学，模拟】
B. Chris and Magic Square time limit per test 2 seconds memory limit per test 256 megabytes input st ...
Codeforces Round #447 (Div. 2) B. Ralph And His Magic Field 数学
题目链接题意:给你三个数n,m,k;让你构造出一个nm的矩阵,矩阵元素只有两个值(1,-1),且满足每行每列的乘积为k,问你多少个矩阵. 解法:首先,如果n,m奇偶不同,且k=-1时,必然无解: 设 ...
Codeforces Round #315 (Div. 2C) 568A Primes or Palindromes? 素数打表+暴力
题目:Click here 题意:π(n)表示不大于n的素数个数,rub(n)表示不大于n的回文数个数,求最大n,满足π(n) ≤ A·rub(n).A=p/q; 分析:由于这个题A是给定范围的,所以 ...
Codeforces Round #657 (Div. 2) B. Dubious Cyrpto（思维，数学）
题目链接题意: m=n⋅a+b−c(n为任意正整数),给出m的值a,b,c的范围l,r(l<=a,b,c<=r),求出a,b,c. 思路: 由推倒知0<=|b-c|<=r-l ...
Codeforces Round #670 (Div. 2) D. Three Sequences（差分序列、数学思维分析）
D. Three Sequences D. Three Sequences 题意:给定一个长度为nnn的序列a[]a[]a[],要求把他分成两个序列b[],c[]b[],c[]b[],c[] a[i] ...
Codeforces Round #568 (Div. 2)C2. Exam in BerSU (hard version)
Codeforces Round #568 (Div. 2)C2. Exam in BerSU (hard version) 贪心+暴力大致题意:N个人考试,每个人花费的时间是a[i],他们总共花费 ...
Codeforces Round #506 (Div. 3)
Codeforces Round #506 (Div. 3) 实习期间事不多,对div3 面向题解和数据编程了一波 A. Many Equal Substrings 题目链接 A题就是找后缀和前缀重合 ...

Codeforces Round #588 (Div. 2) E. Kamil and Making a Stream 数学 + 暴力

文章目录

题意：

思路：

Codeforces Round #588 (Div. 2) E. Kamil and Making a Stream 数学 + 暴力相关推荐

最新文章

热门文章