AtCoder - abc152(D，E)

mldl_

已于 2022-05-17 00:51:23 修改

阅读量380

点赞数 1

分类专栏：数据结构与算法

于 2022-05-17 00:44:48 首次发布

本文链接：https://blog.youkuaiyun.com/mldl_/article/details/124811424

版权

c++ 算法哈希算法

数据结构与算法专栏收录该内容

102 篇文章

订阅专栏

这两篇博客探讨了数学和编程在实际问题解决中的关键作用。第一篇D-Handstand2介绍了如何利用数字首尾统计解决特定问题，通过建立数组统计规律并进行求和计算。第二篇E-Flatten关注于计算数列的调和平均数，采用欧拉筛法找出质因数分布，并利用最小公倍数来解决问题。这两个实例展示了数学思维和高效编程在数据处理中的强大能力。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

D - Handstand 2

思路：
s(x) 表示 x 的首数字，e(x) 表示 x 的尾数字，v[a][b] 表示以 a开头b结尾的数组的个数。
$\sum_{1}^{n}v[e(i)][s(i)]$

#include<bits/stdc++.h>
using namespace std;

int v[10][10];

int s(int x){
	while(x > 9) x /= 10;
	return x;
}

int e(int x){
	return x % 10;
}

int main(){
	
	int n;
	cin >> n;
	for(int i = 1; i <= n; i++) v[s(i)][e(i)]++;
	
	long long res = 0;
	for(int i = 1; i <= n; i++){
		res += v[e(i)][s(i)];
	}
	cout << res << endl;
	
	return 0;
}

E - Flatten

由题意得： $\sum_{1}^{n}\frac{lcm(a_1, ..., a_n)}{a_i}$
$lcm(a_1, ..., a_n)$ 的求法：定义 $cnt_{i,j}$ 表示 $a_i$ 中第 $j$ 个质因数的个数， $\prod_{i=1}^{tot}{pow(pri[i], max(cnt_{1,i}, ... cnt_{n,j}))}$ 。
( $t o t$ 表示欧拉筛之后得到的质数的个数，第 $j$ 个质数指欧拉筛的 $p r i m e [j]$ )

#include<bits/stdc++.h>
#define ll long long
using namespace std;

const ll mod = 1e9 + 7;
const int N = 1e6 + 5;

ll vis[N], pri[N], cnt = 0;
ll a[N], v[N];

ll Qpow(ll a, ll b){
	ll res = 1;
	while(b){
		if(b&1) res *= a;
		res %= mod;
		a *= a;
		a %= mod;
		b >>= 1;
	}
	return res;
}

int main(){
	
	for(int i = 2; i < N; i++){
		if(vis[i] == 0) pri[++cnt] = i;
		for(int j = 1; j <= cnt && i * pri[j] < N; j++){
			vis[i * pri[j]] = 1;
			if(i % pri[j] == 0) break;
		} 
	}
	
	int n;
	cin >> n;
	for(int i = 1; i <= n; i++) cin >> a[i];
	
	for(int i = 1; i <= n; i++){
		ll x = a[i];
		for(int j = 1; pri[j] * pri[j] <= x; j++){
			ll num = 0;
			while(x % pri[j] == 0) num++, x /= pri[j];
			v[pri[j]] = max(num, v[pri[j]]);
		}
		if(x != 0) v[x] = max(v[x], 1ll);
	}

	ll lcm = 1, res = 0;
	for(int i = 1; i <= cnt; i++) lcm = lcm * Qpow(pri[i], v[pri[i]]) % mod;
	for(int i = 1; i <= n; i++) res = (res + lcm * Qpow(a[i], mod-2) % mod) % mod;
	
	cout << res << endl;
		
	return 0;
}