一、题目
二、解法
由于n很大,但是操作数很小,且一开始没有初值,很容易想到动态开点,时间复杂度O(nlogn)。
剩下的问题是如何计算一个修改段内的和,注意到(i−l+1)⋅x%y=(i−l+1)⋅x+y(i−l+1)⋅x(本题解中的除号均为整除),前者等差数列,后者用经典的类欧几里得算法求和,总时间复杂度O(logn),下面详细讲一下这个算法。
类欧几里得算法一般用于解决此类问题:f(a,b,c,n)=∑i=0ncia+b,给定a,b,c,n求f,下面给出算法过程及推导。
当a≥c or b≥c时,f(a,b,c,n)=f(a%c,b%c,c,n)+ca2n(n+1)+(n+1)cb
当a,b≤c,设m=can+b,我们开始推式子:
f(a,b,c,n)=i=0∑nj=0∑m−1[j<cia+b]=j=0∑m−1i=0∑n[i>acj−b+c−1]上面那一步怎么来的呢?我们对括号内的内容做推导:
j<cai+bj≤cai+b−1cj≤ai+b−ccj<ai+b−c+1i>acj−b+c−1然后就推出来了,我们继续推导:
=j=0∑m−1n−acj−b+c−1=nm−f(c,c−b−1,a,m−1)发现上述算法过程类似于辗转相除法,故时间复杂度为log。
咕咕咕