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

大概就是推式子的时候注意有两个边界都是 n ,考虑变成 2*... 之类的。

分块维护 f[ ] 的前缀和。很好的思路是修改一个位置后前缀和数组需要区间加,整块地打上加法标记就行了。

自己本来想维护整块之间的前缀和,还有块内的前缀和;却WA得不行。之后再探究为什么WA吧。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define ll long long
using namespace std;
const int N=4e6+,M=,mod=1e9+;
int n,g[N],phi[N],f[N],pri[N];bool vis[N];
int base,bh[N],s[M],si[M],fl[N];
void upd(int &x){while(x>=mod)x-=mod;while(x<)x+=mod;}
int gcd(int a,int b){return b?gcd(b,a%b):a;}
int pw(int x,int k)
{int ret=;while(k){if(k&)ret=(ll)ret*x%mod;x=(ll)x*x%mod;k>>=;}return ret;}
void init()
{
phi[]=g[]=; int cnt=;
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;break;}
else phi[i*pri[j]]=(ll)phi[i]*phi[pri[j]]%mod;
}
g[i]=(g[i-]+(ll)i*i%mod*phi[i])%mod;//presum!!
}
base=sqrt(n);
for(int i=;i<=n;i++)f[i]=(ll)i*i%mod,fl[i]=(fl[i-]+f[i])%mod;
for(int i=,j=,k=;i<=n;i++,k++)
{
bh[i]=j;if(k==base)k=,j++;
}
/*
for(int i=1,j=1,k=base;i<=n;i++)
{
f[i]=(ll)i*i%mod;bh[i]=j;
si[j]+=f[i]; upd(si[j]); fl[i]=si[j];
if(i==k)s[j]=s[j-1]+si[j],j++,k+=base;
}
*/
}
int calc(int x){int ret=fl[x]+s[bh[x]];upd(ret);return ret;}
int main()
{
int T;scanf("%d%d",&T,&n);init();
int x,y,tn; ll w;
while(T--)
{
scanf("%d%d%lld%d",&x,&y,&w,&tn);//w not %mod!!!
int u=gcd(x,y),d=bh[u];
int tf=w/(x/u)/(y/u)%mod;//not inv /*//also ok
int chg=tf+calc(u-1)-calc(u);upd(chg);
for(int i=u;bh[i]==bh[u];i++)
fl[i]+=chg,upd(fl[i]);
for(int i=bh[u]+1;i<=bh[n];i++)//n not tn!!!
s[i]+=chg,upd(s[i]);
*/
int pl=tf-f[u];upd(pl);f[u]=tf;
for(int i=u;bh[i]==bh[u];i++)
fl[i]+=pl,upd(fl[i]);
for(int i=bh[u]+;i<=bh[n];i++)
s[i]+=pl,upd(s[i]);
/*
si[d]=si[d]-f[u]+tf;upd(si[d]);
for(int i=d;i<=bh[n];i++)
s[i]=s[i-1]+si[i],upd(s[i]);
f[u]=tf; fl[u]=(bh[u-1]==bh[u]?fl[u-1]:0)+f[u];
for(int i=u+1,j=d*base;i<=j;i++)
fl[i]=fl[i-1]+f[i],upd(fl[i]);
*/
int ans=;
for(int i=,j;i<=tn;i=j+)
{
int d=tn/i,sm=; j=tn/d;
/*
if(bh[j]-bh[i]<=1)
for(int l=i;l<=j;l++)
sm+=f[l],upd(sm);
else
{
sm=s[bh[j]-1]-s[bh[i]-1]+fl[j];
upd(sm);
if(bh[i-1]==bh[i])sm-=fl[i-1],upd(sm);
}
ans=(ans+(ll)sm*g[d])%mod;
*/
ans=(ans+(ll)(calc(j)-calc(i-))*g[d])%mod;upd(ans);
}
printf("%d\n",ans);
}
return ;
}

bzoj 4815 [Cqoi2017]小Q的表格——反演+分块的更多相关文章

  1. bzoj 4815: [Cqoi2017]小Q的表格 [数论]

    4815: [Cqoi2017]小Q的表格 题意: 单点修改,查询前缀正方形和.修改后要求满足条件f(a,b)=f(b,a), b×f(a,a+b)=(a+b)*f(a,b) 一开始sb了认为一次只会 ...

  2. bzoj 4815: [Cqoi2017]小Q的表格【欧拉函数+分块】

    参考:http://blog.csdn.net/qq_33229466/article/details/70174227 看这个等式的形式就像高精gcd嘛-所以随便算一下就发现每次修改(a,b)影响到 ...

  3. BZOJ 4815 CQOI2017 小Q的表格 欧拉函数+分块

    题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=4815 题意概述:要认真概述的话这个题就出来了... 分析: 首先分析题目,认真研究一下修 ...

  4. BZOJ 4815 [Cqoi2017]小Q的表格 ——欧拉函数

    把式子化简一波. 发现一个比较厉害的性质:每个点只能影响到行列下标$gcd$与它相同的点. 然后就可以计算$\sum_{g<=k}f(g,g)*\sum_{i<=k}\sum_{j< ...

  5. [CQOI2017]小Q的表格(数论+分块)

    题目描述 小Q是个程序员. 作为一个年轻的程序员,小Q总是被老C欺负,老C经常把一些麻烦的任务交给小Q来处理.每当小Q不知道如何解决时,就只好向你求助. 为了完成任务,小Q需要列一个表格,表格有无穷多 ...

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

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4815 思路就和这里一样:https://blog.csdn.net/leolyun/arti ...

  7. 4815: [Cqoi2017]小Q的表格 莫比乌斯反演 分块

    (Updated 2018.04.28 : 发现公式效果不好,重新处理图片)国际惯例的题面:看到这两个公式,很多人都会想到与gcd有关.没错,最终的结论就是f(a,b)=f(gcd(a,b))*(a/ ...

  8. [CQOI2017]小Q的表格——反演好题

    zhoutb2333的题解 难得一见的新颖反演题. 一眼看可能不是反演题. 修改影响别的,很恶心. 所以考虑化简f的联系式,发现和gcd有关 于是考虑用gcd来表示所有的gcd(a,b)=g的所有f( ...

  9. [BZOJ4815][CQOI2017]小Q的表格 数论+分块

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=4815 题目中所给条件中的$(a,a+b)$和$(a,b)$的关系很瞩目. 然后大家都知道$ ...

随机推荐

  1. 我的json

    { "firstName":[ "xMan" ], "members":[ { "name":"X教授&quo ...

  2. dom&bom的起源,方法,内容,应用

    Document Object Model的历史可以追溯至1990年代后期微软与Netscape的"浏览器大战"(browser wars),双方为了在JavaScript与JSc ...

  3. 牛客-https://www.nowcoder.com/acm/contest/96/H

    链接:https://www.nowcoder.com/acm/contest/96/H来源:牛客网 题目描述 今天qwb要参加一个数学考试,这套试卷一共有n道题,每道题qwb能获得的分数为ai,qw ...

  4. socketserver源码简介

    一.socketserver流程简介 +------------+ | BaseServer | +------------+ | v +-----------+ +----------------- ...

  5. day6-面向对象基础篇

    一.面向对象引子及概念 结合编程的一些理论知识和实践,可以总结出目前存在以下编程模式: 1. 面向过程 按照业务逻辑和实现过程步骤来逐步垒代码,代码编写的逻辑即对应于实际实现的步骤过程,核心是过程两个 ...

  6. 由浅入深了解EventBus:(三)

    原理 EventBus的核心工作机制如下图 在EventBus3.0架构图: EventBus类 在EventBus3.0框架的内部,核心类就是EventBus,订阅者的注册/订阅,解除注册,以及事件 ...

  7. 二进制包形式安装mysql5.7.17

    mysql5.7与mysql5.5的安装方式有些区别,如果按照5.5的方式来安装5.7就会有很多报错信息,下面介绍mysql5.7.17的二进制安装方式. 首先安装依赖包: yum -y instal ...

  8. PHP实现登录功能DEMO

    PHP实现登录的原理是什么呢?就是利用Session实现的,用户访问网站,系统会自动在服务器生成一个Session文件,这个Session可以用来存储用户的登录信息.好了,这是基本储备,我们下面来实现 ...

  9. java.util.Collection List与其子类 Set与其子类

    package com.Collection; import java.util.ArrayList; import java.util.Collection; import java.util.It ...

  10. vue.js 源代码学习笔记 ----- 工具方法 option

    /* @flow */ import Vue from '../instance/index' import config from '../config' import { warn } from ...