bzoj 4815 小Q的表格 —— 反演+分块

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=4815

思路就和这里一样：https://blog.csdn.net/leolyun/article/details/70146612

不知为何乘逆元就错了，必须直接除...不过题目保证了是整数，所以直接除也没问题；

然后重新学习了一下分块的简洁写法，就能A了hhh

代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
int const xn=4e6+,mod=1e9+,base=,xb=;
int n,cnt,pri[xn],phi[xn],g[xn],blk[xn],s[xn],pls[xb];
bool vis[xn];
ll rd()
{
  ll ret=,f=; char ch=getchar();
  while(ch<''||ch>''){if(ch=='-')f=; ch=getchar();}
  while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();
  return f?ret:-ret;
}
ll pw(ll a,int b)
{
  ll ret=;
  for(;b;b>>=,a=(a*a)%mod)if(b&)ret=(ret*a)%mod;
  return ret;
}
int upt(int x){while(x>=mod)x-=mod; while(x<)x+=mod; return x;}
void init()
{
  phi[]=;
  for(int i=;i<=n;i++)
    {
      if(!vis[i])pri[++cnt]=i,phi[i]=i-;
      for(int j=;j<=cnt&&(ll)i*pri[j]<=n;j++)
    {
      vis[i*pri[j]]=;
      if(i%pri[j])phi[i*pri[j]]=(ll)phi[i]*(pri[j]-)%mod;
      else {phi[i*pri[j]]=(ll)phi[i]*pri[j]%mod; break;}
    }
    }
  for(int i=;i<=n;i++)g[i]=(ll)i*i%mod*phi[i]%mod;
  for(int i=;i<=n;i++)g[i]=upt(g[i-]+g[i]);
}
void init2()
{
  for(int i=;i<=n;i++)blk[i]=(i-)/base+;
  for(int i=;i<=n;i++)s[i]=(s[i-]+(ll)i*i)%mod;
}
void change(int d,ll x)//add x
{
  for(int i=d;blk[i]==blk[d]&&i<=n;i++)s[i]=upt(s[i]+x);
  for(int i=blk[d]+;i<=blk[n];i++)pls[i]=upt(pls[i]+x);
}
int cal(int ps)
{
  return upt(pls[blk[ps]]+s[ps]);
}
int gcd(int a,int b){return b?gcd(b,a%b):a;}
int main()
{
  int m=rd(); n=rd();
  init(); init2(); ll x;
  for(int i=,a,b,k;i<=m;i++)
    {
      a=rd(); b=rd(); x=rd(); k=rd();
      int d=gcd(a,b);
      //int t1=(ll)d*pw(a,mod-2)%mod,t2=(ll)d*pw(b,mod-2)%mod;
      ll t=((ll)a/d)*(b/d); x=(x/t)%mod;//
      change(d,upt(x-cal(d)+cal(d-)));
      //change(d,x*t1%mod*t2%mod);
      int ans=;
      for(int j=,nxt;j<=k;j=nxt+)
    {
      nxt=k/(k/j);
      ans=(ans+(ll)g[k/j]*upt(cal(nxt)-cal(j-)))%mod;
    }
      printf("%d\n",ans);
    }
  return ;
}