BZOJ3529: [Sdoi2014]数表

对程序看半天才发现公式打错了我怎么又把×打成+了….
Mod 2^31 可以利用int的自然溢出，最后&(2^31-1)—–PoPoQQQ
这题我公式好像有多种推法但是能保证不大于a的条件的好像只有一种

---

---

离线处理询问，要计算不大于a的，将询问按a排序，推出一个和每个位置权值有关的函数按照从小到大插值
先不考虑a的限制，但最后推得的式子要有方法能确保对于每个询问计算的是不大于a的和
因为对于一对i，j，题目要求的其实是gcd(i,j)的约数和，所以设$F(n)$表示n的约数和
$F(n)=\sum_{d|n}d=1*n$（$*$代表狄利克雷卷积，下文不再解释）
那么$F(n)$是积性函数，可以筛出来，我们再回到题目
设d为i，j的gcd，那么题目要求的就是：
$\sum_{d=1}^{min(n,m)}F(d)\sum_{i=1}^{\lfloor n/d\rfloor}\sum_{j=1}^{\lfloor m/d\rfloor}1[gcd(i,j)=1]->$
$\sum_{d=1}^{min(n,m)}F(d)\sum_{i=1}^{\lfloor n/d\rfloor}\sum_{j=1}^{\lfloor m/d\rfloor}\sum_{k|gcd(i,j)}μ(k)->$
$\sum_{d=1}^{min(n,m)}F(d)\sum_{k}μ(k)\sum_{i=1}^{\lfloor n/dk\rfloor}\sum_{j=1}^{\lfloor m/dk\rfloor}1->$
$\sum_{d=1}^{min(n,m)}F(d)\sum_{k}μ(k)\lfloor \frac{n}{dk}\rfloor\lfloor \frac{m}{dk}\rfloor->$
枚举D=dk
$\sum_{D=1}^{min(n,m)}\sum_{d|D}F(d)×μ(\frac D d)\lfloor \frac{n}{D}\rfloor\lfloor \frac{m}{D}\rfloor$ 前面的部分又是一个卷积的形式
所以设$g=F*μ$
得：$\sum_{D=1}^{min(n,m)}g(D)\lfloor \frac{n}{D}\rfloor\lfloor \frac{m}{D}\rfloor$ 然后就可以$\sqrt n$求解了

再考虑a，题目要求的是$F$不大于a，所以将所有的$F$排序，从小到大插入，插入一个$F$就暴力更新他对所有能贡献到的$g$的贡献，在树状数组上加，询问时树状数组询问区间和，这样复杂度是$O(nlog^2n + Q\sqrt n log\ n)$

---

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

const int maxn = 100001;
struct node
{
    int i,f;
}a[maxn];
bool cmp(node x,node y){return x.f<y.f;}
int p[maxn],pri,miu[maxn];
int f[maxn],g[maxn];
bool v[maxn];

int lowbit(int x){return x&-x;}
void add(int x,int c){while(x<maxn){g[x]+=c;x+=lowbit(x);}}
int query(int x){int ret=0;while(x){ret+=g[x];x-=lowbit(x);}return ret;}

void get_prime()
{
    f[1]=miu[1]=1; a[1].f=1;a[1].i=1;
    for(int i=2;i<maxn;i++)
    {
        if(!v[i])
        {
            p[++pri]=i;
            miu[i]=-1;
            f[i]=1+i;
        }
        for(int j=1;j<=pri&&p[j]*i<maxn;j++)
        {
            int k=p[j]*i;
            v[k]=true;
            if(i%p[j]==0)
            {
                miu[k]=0;
                int temp=i,tmp=p[j];
                while(temp%p[j]==0) temp/=p[j],tmp*=p[j];
                if(temp==1) f[k]=f[i]+k;
                else f[k]=f[temp]*f[tmp];
                break;
            }
            miu[k]=-miu[i];
            f[k]=f[i]*f[p[j]];
        }
        a[i].i=i;a[i].f=f[i];
    }
}
struct Query
{
    int n,m,a,i,ans;
}q[maxn]; int Q;
bool cmp2(Query x,Query y){return x.a<y.a;}
bool cmp3(Query x,Query y){return x.i<y.i;}
int n,m;

int main()
{
    get_prime(); sort(a+1,a+maxn,cmp);

    scanf("%d",&Q);
    for(int i=1;i<=Q;i++) 
        scanf("%d%d%d",&q[i].n,&q[i].m,&q[i].a),q[i].i=i;
    sort(q+1,q+Q+1,cmp2);

    int pos=0;
    for(int i=1;i<=Q;i++)
    {
        while(a[pos+1].f<=q[i].a&&pos<maxn-1)
        {
            pos++;
            for(int j=a[pos].i;j<maxn;j+=a[pos].i)
                add(j,a[pos].f*miu[j/a[pos].i]);
        }
        n=q[i].n; m=q[i].m;
        int l=1,r=min(n,m);
        while(l<=r)
        {
            int nex=min(n/(n/l),m/(m/l));
            int s1=(query(nex)-query(l-1));
            q[i].ans+=(n/l)*(m/l)*s1;
            l=nex+1;
        }
    }

    sort(q+1,q+Q+1,cmp3);
    for(int i=1;i<=Q;i++)printf("%d\n",q[i].ans&INT_MAX);

    return 0;
}