bzoj3309

题意：
定义f(n)为n所含质因子的最大幂指数。求
$\sum_{i=1}^n\sum_{j=1}^mf(gcd(i,j))$
T<=10000
1<=n,m<=10^7

#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<iostream>
#define N 11000000
#define LL long long
using namespace std;
int p[N],pl,a[N],d[N],n,g[N];
LL ans;
bool b[N];
void get_p()
{
    for(int i=2;i<N;i++)
    {
        if(b[i]==0) {p[++pl]=i;a[i]=1;d[i]=i;g[i]=1;}
        for(int j=1;j<=pl;j++)
        {
            if(i*p[j]>=N) break;
            b[i*p[j]]=1;
            if(i*p[j]==4)
            {int oo=1;}
            if(i%p[j])
            {
                a[i*p[j]]=1;d[i*p[j]]=p[j];
                if(a[i]!=1) g[i*p[j]]=0;
                else g[i*p[j]]=-g[i];
            }
            else
            {
                a[i*p[j]]=a[i]+1;d[i*p[j]]=d[i]*p[j];
                int t=i/d[i];
                if(t==1) g[i*p[j]]=1;
                else if(a[t]==a[i*p[j]]) g[i*p[j]]=-g[t];
                else g[i*p[j]]=0;
                break;
            }
        }
    }
    for(int i=1;i<N;i++) g[i]=g[i]+g[i-1];
}
int main()
{
    get_p();
    int z;scanf("%d",&z);
    while(z--)
    {
        int n,m,nex;scanf("%d%d",&n,&m);
        ans=0;
        if(n>m) swap(n,m);
        for(int i=1;i<=n;i=nex+1)
        {
            nex=min(n/(n/i),m/(m/i));
            ans+=(LL)(n/i)*(m/i)*(g[nex]-g[i-1]);
        }
        printf("%lld\n",ans);
    }
    return 0;
}

题解：
先大力化式子
枚举i，j的gcd
$\sum_{d=1}^nf(d)\sum_{i=1}^{\frac{n}{d}}\sum_{j=1}^{\frac{m}{d}}\epsilon((i,j))$
$\sum_{d=1}^nf(d)\sum_{i=1}^{\frac{n}{d}}\sum_{j=1}^{\frac{m}{d}}\sum_{k|i,k|j}\mu(k)$
$\sum_{k=1}^n\mu(k)\sum_{d=1}^{\frac{n}{k}}f(d)\sum_{i=1}^{\frac{n}{dk}}\sum_{j=1}^{\frac{m}{dk}}$
枚举D=dk
$\sum_{D=1}^n\frac{n}{D}\frac{m}{D}\sum_{d|D}f(d)\mu(\frac{D}{d})$
设$g=f*\mu$
$\sum_{D=1}^n\frac{n}{D}\frac{m}{D}g(D)$
算出g后就可以分块解决。
觉得g不是一个积性函数，看了一眼题解(捂脸)“由$\mu$和$f$性质可得”
那就大力推一下吧。。
设$n=p_1^{k_1}p_2^{k_2}...p_t^{k_t}$
k1>=k2>=…>=kt
我们只考虑非0的μ和对应的f
当多个ki值相等的情况，我们把它当作用最小那个i的情况
使f取值为k1，强制p1不选，贡献为
$k_1\sum_{i=0}^{t-1}(-1)^iC_{t-1}^i$
$k_1(1-1)^{t-1}$
似乎为0
注意如果对于kj，存在ki>kj，那kj永远不会作为f，因为ki-1>=kj
所以，kj能做f的情况为k1=k2=…=kj且k1到kj-1都被选了
类似，贡献为
$(-1)^{j-1}kj\sum_{i=0}^{t-j}(-1)^iC_{t-1}^i$
$(-1)^{j-1}kj(1-1)^{t-j}$
如果存在k1=k2=…=kj,kj>kj+1那强制p1到pj都选后，后面所有情况的f依然为k1-1，后面所有情况贡献为
$(-1)^j(k_1-1)\sum_{i=0}^{t-k}(-1)^iC_{t-1}^i$
$(-1)^{j}kj(1-1)^{t-j}$
上面所有式子在j不等于t时都为0
唯一不为0的情况就是k1=k2=…=kt时的
$(-1)^{t-1}k_t(1-1)^{t-t}$
以及全部p都被选后的
$(-1)^{t}(k_t-1)$
相加即为$(-1)^{t+1}$
于是可以用线性筛求出g。
具体来说就是对于每个数n求出他的最小质因子p最多有k次方，记为a(n)
比较$a(n)和a(\frac{n}{p^k})$即可从$g(\frac{n}{p^k})推出g(n)$
这是一个递推的关系，用心感受一下就知道很对了。。