对原序列中所有(x,m)=d的ax+v,询问前缀和
对这个修改推一波柿子(以下省略下取整)
v∑nx=1(x,m)=d−>
v∑n/dx=1(x,m/d)=1−>
v∑n/dx=1∑D|(x,M)μ(D)−>
v∑D|Mμ(D)∑n/Ddx=1
因为是对序列中所有符合条件的位置修改,我们维护序列中所有
Dd|x
的x的共同增量
f(Dd)
询问就变成了
∑ji=1∑d|if(d)−>
∑jd=1f(d)∑j/di=1
树状数组资瓷单点修改,求前缀和,查询的时候用分块的思想跳,就可以做到修改查询
O(nn√logn)
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;
}