SPOJ GCDEX (数论)

转载请注明出处，谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents by---cxlove

题意：求sigma (gcd (i , j)) 1 <= i < j <= n

和POJ 2480类似，如果枚举j，求的话，还是会TLE的。。。

考虑sigma(gcd (i , n)) = sigma (d * phi[n / d]) d | n。

做法同样是先预处理出phi，然后枚举gcd = d。

O(n)预处理，O(1)查询

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <set>
#include <map>
#include <string>
#include <vector>
#include <queue>
#include <stack>
#define lowbit(x) (x & (-x))
#define Key_value ch[ch[root][1]][0]
#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;
typedef long long LL;
const int N = 1000005;
int phi[N];
LL ans[N];
void init () {for (int i = 2 ; i < N ; i ++) {if (phi[i]) continue;for (int j = i ; j < N ; j += i) {if (phi[j] == 0) phi[j] = j;phi[j] = phi[j] / i * (i - 1);}}for (int i = 2 ; i < N ; i ++)ans[i] = phi[i];for (int i = 2 ; i * i < N ; i ++) {ans[i * i] += 1LL * phi[i] * i;for (int j = i + 1 ; j * i < N ; j ++) {ans[j * i] += 1LL * phi[i] * j + 1LL * phi[j] * i;}}for (int i = 1 ; i < N ; i ++)ans[i] += ans[i - 1];
}
int main () {#ifndef ONLINE_JUDGEfreopen ("input.txt" , "r" , stdin);// freopen ("output.txt" , "w" , stdout);#endifint n ;init ();while (cin >> n && n) {cout << ans[n] << endl;}return 0;
}

转载于:https://www.cnblogs.com/pangblog/p/3331354.html

SPOJ GCDEX (数论)相关推荐

BZOJ 2226 [Spoj 5971] LCMSum 最大公约数之和 | 数论
BZOJ 2226 [Spoj 5971] LCMSum 这道题和上一道题十分类似. \[\begin{align*} \sum_{i = 1}^{n}\operatorname{LCM}(i, n) ...
数论（一）——素数，GCD，LCM
这是一个数论系列:) 一.素数 ×费马小定理 Theorem: 设 p 是一个素数,a 是一个整数且不是 p 的倍数,那么很遗憾,费马小定理的逆定理是不成立的.对 a = 2,满足的非素数 n 是存 ...
bzoj 2588 Spoj 10628. Count on a tree （可持久化线段树）
Spoj 10628. Count on a tree Time Limit: 12 Sec Memory Limit: 128 MB Submit: 7669 Solved: 1894 [Sub ...
BZOJ 2780: [Spoj]8093 Sevenk Love Oimaster( 后缀数组 + 二分 + RMQ + 树状数组 )
全部串起来做SA, 在按字典序排序的后缀中, 包含每个询问串必定是1段连续的区间, 对每个询问串s二分+RMQ求出包含s的区间. 然后就是求区间的不同的数的个数(经典问题), sort queries ...
SPOJ GSS3-Can you answer these queries III-分治+线段树区间合并
Can you answer these queries III SPOJ - GSS3 这道题和洛谷的小白逛公园一样的题目. 传送门: 洛谷 P4513 小白逛公园-区间最大子段和-分治+线段树区间 ...
【数论总结】-----励志写好一篇数论总结↖(^ω^)↗//正在施工...未完工
近期学了学数论,来写一波总结吧. (1)排列组合,比较基础的东西了吧.//只写个概念吧,(逃: 概念:就是从n个不同元素中,任取m(m≤n)个元素并成一组,叫做从n个不同元素中取出m个元素的一个组合: ...
SPOJ ATOMS - Atoms in the Lab
题目链接:http://www.spoj.com/problems/ATOMS/ 题目大意:有N个原子,他们每秒分裂成K个新原子,新原子也能继续分裂.问如果要控制他的数量为M以内,应在什么时候使其停止 ...
解题报告（十八）数论题目泛做（Codeforces 难度：2000 ~ 3000 + ）
整理的算法模板合集: ACM模板点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划繁凡出品的全新系列:解题报告系列 -- 超高质量算法题单,配套我写的超高质量的题解和代码,题目难度不一 ...
《算法竞赛中的初等数论》（五）正文 0x50筛法（ACM / OI / MO）（十五万字符数论书）
整理的算法模板合集: ACM模板点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划写在最前面:本文部分内容来自网上各大博客或是各类图书,由我个人整理,增加些许见解,仅做学习交流使用,无 ...

SPOJ GCDEX (数论)

SPOJ GCDEX (数论)相关推荐

最新文章

热门文章