BZOJ3853: GCD Array

对原序列中所有(x,m)=d的ax+v，询问前缀和
对这个修改推一波柿子（以下省略下取整）
$v\sum_{x=1}^n(x,m)=d->$
$v\sum_{x=1}^{n/d}(x,m/d)=1->$
$v\sum_{x=1}^{n/d}\sum_{D|(x,M)}μ(D)->$
$v\sum_{D|M}μ(D)\sum_{x=1}^{n/Dd}$
因为是对序列中所有符合条件的位置修改，我们维护序列中所有$Dd|x$的x的共同增量$f(Dd)$
询问就变成了$\sum_{i=1}^j\sum_{d|i}f(d)->$
$\sum_{d=1}^jf(d)\sum_{i=1}^{j/d}$
树状数组资瓷单点修改，求前缀和，查询的时候用分块的思想跳，就可以做到修改查询$O(n\sqrt nlogn)$

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
#define lowbit(x) x&(-x)
using namespace std;

inline void read(int &x)
{
    char c; while(!((c=getchar())>='0'&&c<='9'));
    x=c-'0';
    while((c=getchar())>='0'&&c<='9') (x*=10)+=c-'0';
}
const int maxn = 210000;

int mu[maxn],p[maxn],pri;
bool v[maxn];
void pre()
{
    mu[1]=1;
    for(int i=2;i<maxn;i++)
    {
        if(!v[i]) p[++pri]=i,mu[i]=-1;
        for(int j=1,k=p[j]*i;k<maxn;j++,k=p[j]*i)
        {
            v[k]=true;
            if(i%p[j]==0) break;
            mu[k]=-mu[i];
        }
    }
}

int n,m;
ll s[maxn];
void add(int x,int c){for(;x<=n;x+=lowbit(x))s[x]+=(ll)c;return;}
ll query(int x){ll re=0;for(;x;x-=lowbit(x))re+=s[x];return re;}

ll ans;

int main()
{
    pre();
    int tcase=0;
    while(scanf("%d%d",&n,&m)!=EOF&&n&&m)
    {
        printf("Case #%d:\n",++tcase);
        for(int i=1;i<=n;i++) s[i]=0;
        for(int i=1;i<=m;i++)
        {
            int ti; read(ti);
            if(ti==1)
            {
                int x,d,v; read(x); read(d); read(v);
                if(x%d) continue;
                x/=d;
                int u=sqrt(x);
                for(int j=1;j<=u;j++) if(x%j==0)
                {
                    int k=j;
                    add(k*d,mu[k]*v);
                    if(j*j==x)break;
                    k=x/j;
                    add(k*d,mu[k]*v);
                }
            }
            else
            {
                ans=0;
                int x; read(x);
                int l,r;
                for(l=1;l<=x;l=r+1)
                {
                    r=x/(x/l);
                    ans+=(query(r)-query(l-1))*(x/l);
                }
                printf("%lld\n",ans);
            }
        }
    }

    return 0;
}