这题真是过了n年才a。最早是在2010年北大培训比赛上看到的这题,当时我不会,竹教主也不会,但他记下来了,研究一段时间后就会了,还把这题加到我校oj上。过了这么多年,我上网搜,关于这个问题的解题报告还是没有,于是我花了几天时间做了出来,发布此解题报告。

题目是要求从1到n的所有数与n的最小公倍数的和,再除了n。第一眼看到数据范围,就知道是不能硬做(即从1到n依次算一遍)的。那么,为了找规律,得把公式化一化。

f(n) = [lcm(1,n)+lcm(2,n)+……+lcm(n,n)]/n
   = 1/gcd(1,n)+2/gcd(2,n)+……+n/gcd(n,n)

于是f(n)化成了一堆数的和。

显然,每一个gcd(k,n)都是能整除n的,按gcd(k,n)的值进行分类合并,就可以把f(n)整理成以下形式

f(n) = (……)/1 + (……)/2 +……+ (……)/d

这里d就是所有能整除n的整数。而(……)/1的分子为所有与n互质的数的和;(……)/2的分子为所有与n的最大公约数为2,也就是与n/2互质的数的和;依次类推,得到公式

有了这个公式,虽然是前进了一大步,但依然没法解题。因为这还是必须枚举n的所有约数,而这个复杂度依然为O(n)。

接下来就是比较核心的一步,构造。

f(n) = (g(n) + 1) / 2,求得g(n)就立即可得f(n)了。

这里构造g(n)的原因,是因为g(n)有积性性质。

对于任意gcd(m, n) = 1,有g(m * n) = g(m) * g(n)。证明略掉,读者自证。

有了积性性质就好办了,将n质因数分解成pi^ci相乘的形式,n = ∏(pi^ci),则g(n) = g(∏(pi^ci)) = ∏(g(pi^ci))

g(pi^ci)很好计算。这里的d就是1, pi, pi^2, pi^3, ..., pi^ci,而phi(pi^ci)=(pi-1)*pi^(ci-1),合起来就是个等比数列,可以推出来一个公式。也就是说可以O(1)时间算出g(pi^ci)。

如此一来,整个问题就可解了。

不过这题还是很变态的,因为数据就有50000组,仅分解质因数写得不好,都可能挂掉。最后我还是过了,代码如下:

/*
* bjfu1109
* Author : ben
*/
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <ctime>
#include <iostream>
#include <algorithm>
#include <queue>
#include <set>
#include <map>
#include <stack>
#include <string>
#include <vector>
#include <deque>
#include <list>
#include <functional>
#include <numeric>
#include <cctype>
using namespace std;
#ifdef ON_LOCAL_DEBUG
#else
#endif
typedef long long LL;
typedef int typec; LL getPow(int a, int b) {
LL res, temp;
res = , temp = (LL) a;
while (b) {
if (b & ) {
res = res * temp;
}
b >>= ;
temp = temp * temp;
}
return res;
} int get_int() {
int res = , ch;
while (!((ch = getchar()) >= '' && ch <= '')) {
if (ch == EOF)
return -;
}
res = ch - '';
while ((ch = getchar()) >= '' && ch <= '')
res = res * + (ch - '');
return res;
} const int N = ;
bool isPrime[N + ];//多用两个元素以免判断边界
vector<int> pt;
void init_prime_table() {
memset(isPrime, true, sizeof(isPrime));
int p = , q, del;
double temp;
while (p <= N) {
while (!isPrime[p]) { p++; }
if (p > N) {//已经结束
break; }
temp = (double) p;
temp *= p;
if (temp > N)
break;
while (temp <= N) {
del = (int) temp; isPrime[del] = false;
temp *= p; }
q = p + ;
while (q < N) {
while (!isPrime[q]) { q++; }
if (q >= N) { break;}
temp = (double) p;
temp *= q;
if (temp > N) break;
while (temp <= N) {
del = (int) temp;
isPrime[del] = false;
temp *= p;
}
q++;
}
p++;
}
for (int i = ; i <= N; i++) {
if (isPrime[i]) {
pt.push_back(i);
}
}
} /**
* p为素数表,至少应该包括sqrt(N)以内的所有素数,f保存结果
* f[i].first表示N的一个素因数,f[i].second为这个素因子的个数
* typec可以是int、long、long long等
*/
typedef vector<pair<typec, int> > FactorList;
void get_prime_factor(const typec &N_, FactorList &f, const vector<int> &p) {
int i, t, n, pl = p.size();
typec N = N_;
f.clear();
for(i = ; i < pl; i++) {
t = p[i];
if (t * t > N_) {
break;
}
if(N % t == ) {
n = ;
while(N % t == ) {
n++; N /= t;
}
f.push_back(make_pair(t, n));
}
if(N == ) { break; }
}
if(N > ) {
f.push_back(make_pair(N, ));
}
} inline LL g(LL p, int c) {
LL ret = getPow(p, * c) - ;
ret = ret / (p + ) * p + ;
return ret;
} FactorList fl;
LL g(int n) {
get_prime_factor(n, fl, pt);
int len = fl.size();
LL ret = ;
for (int i = ; i < len; i++) {
ret *= g(fl[i].first, fl[i].second);
}
return ret;
} int main() {
#ifdef ON_LOCAL_DEBUG
freopen("data.in", "r", stdin);
#endif
LL ans, N;
init_prime_table();
while ((N = get_int()) > ) {
ans = g(N) + ;
ans /= ;
printf("%I64d\n", ans);
}
return ;
}

bjfu1109 最小公倍数和的更多相关文章

  1. 求N个数的最大公约数和最小公倍数(转)

    除了分解质因数,还有另一种适用于求几个较小数的最大公约数.最小公倍数的方法 下面是数学证明及算法实现 令[a1,a2,..,an] 表示a1,a2,..,an的最小公倍数,(a1,a2,..,an)表 ...

  2. C语言 · 最小公倍数

    问题描述 编写一函数lcm,求两个正整数的最小公倍数. 样例输入 一个满足题目要求的输入范例.例:3 5 样例输出 与上面的样例输入对应的输出.例: 数据规模和约定 输入数据中每一个数的范围. 例:两 ...

  3. Java程序设计之最大公约数和最小公倍数

    题目:输入两个正整数number1和number2,求其最大公约数和最小公倍数. 算法:较大数和较小数取余,较小数除余数,一直到余数为0时,为最大公约数(辗转相除法):最大公倍数numbe1*numb ...

  4. 最大公约数和最小公倍数--java实现

    代码: //最大公约数 public int gcd(int p,int q){ if(q == 0) return p; return gcd(q, p % q); } //最小公倍数 public ...

  5. python 最小公倍数

    最小公倍数 求解两个整数(不能是负数)的最小公倍数 方法一:穷举法 def LCM(m, n): if m*n == 0: return 0 if m > n: lcm = m else: lc ...

  6. 输入两个正整数m和n,求其最大公约数和最小公倍数

    public static void main(String[] args){  Scanner sc = new Scanner (System.in);  int a,b;  System.out ...

  7. Java编写最大公约数和最小公倍数

    package javaapplication24; class NegativeIntegerException extends Exception{ String message; public ...

  8. poj 3101Astronomy(圆周追击+分数最小公倍数)

    /* 本题属于圆周追击问题: 假设已知两个圆周运动的物体的周期分别是a ,b, 设每隔时间t就会在同一条直线上 在同一条直线上的条件是 角度之差为 PI ! 那么就有方程 (2PI/a - 2PI/b ...

  9. 【codevs1012】最大公约数和最小公倍数

    题目描述 Description 输入二个正整数x0,y0(2<=x0<100000,2<=y0<=1000000),求出满足下列条件的P,Q的个数 条件:  1.P,Q是正整 ...

随机推荐

  1. js模块化开发

    主要有两个:一个是sea.js,另一个是require.js

  2. java.util.concurrent包API学习笔记

    newFixedThreadPool 创建一个固定大小的线程池. shutdown():用于关闭启动线程,如果不调用该语句,jvm不会关闭. awaitTermination():用于等待子线程结束, ...

  3. BCB常用文件与字符串函数

    VCL库函数简介 一.BORLAND C++ BUILDER VCL的内存管理函数 1. AllocMem 在队中分配指定字节的内存块,并将分配的每一个字节初始化为 0.函数原型如下: void * ...

  4. PHP裁剪图片

    PHP裁剪图片 $src_path = '1.jpg'; //创建源图的实例 $src = imagecreatefromstring(file_get_contents($src_path)); / ...

  5. LA 6047 Perfect Matching 字符串哈希

    一开始我用的Trie+计数,但是不是计多了就是计少了,后来暴力暴过去的…… 看了别人的代码知道是字符串哈希,但是仍有几个地方不理解: 1.26^500溢出问题 2.没考虑哈希碰撞? 跪求指点! #in ...

  6. [转载]浅析Java中的final关键字

    浅析Java中的final关键字 谈到final关键字,想必很多人都不陌生,在使用匿名内部类的时候可能会经常用到final关键字.另外,Java中的String类就是一个final类,那么今天我们就来 ...

  7. apk反编译(5)用apktool重新生成一个未签名的apk

    用apktool反编译apk后,得到一个目录,里面有smali文件,可以对其修改,然后用apktool重新生成一个未签名的apk. 如,把smali文件中的广告部分去掉或改成自己的. 命令如下:与破解 ...

  8. Android下 ionic view 无法登录

    ionic view一个超棒工具,它是测试 ionic 框架搭建项目的app软件. 在它的官网有iphone 和 android 版本的下载地址.但是,这里只有在 google play 里面才有,而 ...

  9. eclipse mingw cpp开发环境

    Eclipse开发c++ 对比:微软的VC++6.0:太老了,对win7兼容不好, 现在微软的Visual Studio:安装包太大,好几个G,装了一堆你不需要的东西,要钱,教 育版申请麻烦 DOS下 ...

  10. POJ 2828 Buy Tickets (线段树 单点更新 变形)

    题目链接 题意:有N个人排队,给出各个人想插队的位置和标识,要求输出最后的序列. 分析:因为之前的序列会因为插队而变化,如果直接算时间复杂度很高,所以可以用 线段树逆序插入,把序列都插到最后一层,le ...