LightOJ 1375 - LCM Extreme 莫比乌斯反演或欧拉扩展

题意:给出n [1,3*1e6] 求
并模2^64。
思路:先手写出算式

观察发现可以化成

那么关键在于如何求得i为1~n的lcm(i,n)之和。可以知道lcm(a,b)为ab/gcd(a,b)
变换得(a/gcd) (b/gcd)gcd 由于GCD的性质，可以知道a/gcd 与 b/gcd是互质的两个质数。由此可以想到应用欧拉函数，并且由性质能够证明 nphi(n)/2为小于n所有与n互质数之和（证明:已知一个质数p那么显然n-p与它互质，那么phi(n)中有phi(n)/2对数，每对数和为n）
故
设n/gcd(I,n)为d则
由此题目化成枚举d即可。还需注意格式的控制转换，本题需要模2^64 只需设unsigned long long 溢出即模，内存限制是刚好卡住的。

#include <stdio.h>
#include <iostream>
#include <string.h>
#include <algorithm>
#include <utility>
#include <vector>
#include <map>
#include <set>
#include <string>
#include <stack>
#include <queue>
#define LL unsigned long long
#define MMF(x) memset((x),0,sizeof(x))
#define MMI(x) memset((x), INF, sizeof(x))
using namespace std;

const int INF = 0x3f3f3f3f;
const int N = 1e6+10;

int eul[3*N];
LL fa[3*N];
LL ans[3*N];
void eular()
{
    MMF(eul);
    MMF(fa);

    eul[1] = 1;
    for(int i = 2; i < 3*N; i++)
    {
        if(!eul[i])
        {
            for(int j = i; j < 3*N; j+=i)
            {
                if(!eul[j])
                    eul[j] = j;
                eul[j] = eul[j]/i * (i-1);
            }
        }
    }
    ans[0] = ans[1] = 0;
    for(LL i = 2; i < 3*N; i++)
    {
        for(LL j = i; j < 3*N; j += i)
        {
            LL t =  j * eul[i] / 2;
            fa[j] += i* t;
        }
        ans[i] = ans[i-1] + fa[i];
    }
}


int main()
{
    eular();
    int T;
    int cnt = 0;
    scanf("%d", &T);
    while(T--)
    {
        LL n;
        scanf("%llu", &n);
        printf("Case %d: %llu\n", ++cnt, ans[n]);
        //printf("%d\n",eul[3000000]);
    }
    return 0;
}
/*
5
2
10
13
100000
3000000
**/

转载于:https://www.cnblogs.com/Yumesenya/p/6008026.html