SPOJ LCMSUM&GCDEX

爱神博客上的两道数论题。SPOJ做的苦不堪言。。。灰常难受的说，各种卡时。。。。。然后因为一个小错误，WA了很长时间。。。。PS:代码能力已经真心不忍直视了，导致长沙网络赛没有过前70。。。。。。。不多说

先是LCMSUM这道题

f(n) = sigma{lcm(i,n)} =sigma{ (i*n)/gcd(i,n)}

分母都不一样，所以要去分母，将n与gcd约分后，可以看到，对于gcd(i,n) = d的i而言

每一个i都是与n/d互质且小于n/d的数的d倍

根据数论知识，令G(x) 为 比x小且与它互质的数的和 G(x) = phi(x)*x/2

所以，所求的f(n) = sigma{t * t * phi(t) /2} * g,其中t =n/g

g是有可能成为n的最大公约数的数，那就是n的所有约数了。

但是，由于T太多，复杂度是O(T*sqrt(n)),超时。。。。555555

也就是说对于每一个查询必须是O(1)

就要预处理f(1)~f(n)的所有值。。。。。（好难的样子）

这个时候不妨用图像解决问题，f(i) = F(t,g)，其中t*g<=N,1<=g<N我们要求所有t*g<=N的F(t,g)

好了，反函数图像出来了，根据反函数对称性，预处理复杂度O(n+n/2+n/3+....+sqrt(n))，这样的话就OK了。。

附代码：

#include <iostream>
#include <cstdio>
  
#define N 1000005
#define ll long long
#define ss(a) scanf("%d",&a)

using namespace std;

ll euler[N],f[N],ans[N];

void prime()
{
    int i,j;
    for (i=2;i<N;i++)
        if (!euler[i])
        {
            for (j=i;j<N;j+=i)
            {
                if (!euler[j]) euler[j]=j;
                euler[j]=euler[j]*(i-1)/i;
            } 
        }       
}

void init()
{
    int i,j;
    for (i=1;i<N;i++) f[i]=f[i]+(ll)i*i*euler[i]/2+i;
    for (i=2;i*i<N;i++)
    {
        f[i*i]+=(ll)i*i*i*euler[i]/2;
        for (j=i+1;j*i<N;j++)
        {
            f[i*j]=f[i*j]+(ll)j*j*i*euler[j]/2+(ll)i*i*j*euler[i]/2;
        }
    }
}

int main()
{
    prime();
    init();
    int T,n;
    ss(T);
    while (T--)
    {
        ss(n);
        printf("%lld\n",f[n]);
    }
    return 0;
}

GCDEX:

求sigma{gcd(i,j)},1<=i,j<=n

先固定一个j，求sigma{gcd(i,n)}

思路和上一道题差不多，都是枚举gcd(i,n)，发现符合要求的i都是与n/g互质的个数

sigma{gcd(i,n)} = sigma{euler(x)},其中x|n

Ans = sigma{gcd(i,j)} = sigma{euler(d)*(n/d)}  在1~n中,要统计n/d次euler(d)

同样，由于时间复杂度问题，要计算ANS[i]

Ans[i] = sigma{euler(d) * t } =G(d,t)，而我们必须把平面上所有(d,t)所对应的函数值统计到ans[i]中

由于d*t<=n,和上道题差不多，还是反函数图像处理法。附代码：

#include <iostream>
#include <cstdio>
#include <cmath>
 
#define N 1000005
#define ll long long
#define ss(a) scanf("%d",&a)

using namespace std;
 
ll euler[N],f[N],ans[N];

void prime()
{
    int i,j;
    for (i=2;i<N;i++)
        if (!euler[i])
        {
            for (j=i;j<N;j+=i) 
            {
                if (!euler[j]) euler[j]=j;
                euler[j]=euler[j]*(i-1)/i;
            }
        }
}

void init()
{
    int i,j;
    ans[0]=0;
    for (i=1;i<N;i++) f[i]=euler[i];
    for (i=2;i*i<N;i++)
    {
        f[i*i]+=i*euler[i];
        for (j=i+1;j*i<N;j++)
        {
            f[i*j]+=i*euler[j]+j*euler[i];
        }
    }
    for (i=1;i<N;i++) ans[i]=ans[i-1]+f[i];
}       

int main()
{
    int n;
    prime();
    init();
    while (ss(n)!=EOF&&n) printf("%lld\n",ans[n]);
    return 0;
}