阿凡达(类欧几里得算法)

六眼飞鱼酱① 提交于 2020-01-23 20:52:03

一、题目

在这里插入图片描述

二、解法

由于nn很大,但是操作数很小,且一开始没有初值,很容易想到动态开点,时间复杂度O(nlogn)O(n\log n)

剩下的问题是如何计算一个修改段内的和,注意到(il+1)x%y=(il+1)x+(il+1)xy(i-l+1)\cdot x\% y=(i-l+1)\cdot x+\frac{(i-l+1)\cdot x}{y}(本题解中的除号均为整除),前者等差数列,后者用经典的类欧几里得算法求和,总时间复杂度O(logn)O(\log n),下面详细讲一下这个算法。

类欧几里得算法一般用于解决此类问题:f(a,b,c,n)=i=0nia+bcf(a,b,c,n)=\sum_{i=0}^n \frac{ia+b}{c},给定a,b,c,na,b,c,nff,下面给出算法过程及推导。

ac or bca\geq c\space or\space b\geq c时,f(a,b,c,n)=f(a%c,b%c,c,n)+acn(n+1)2+(n+1)bcf(a,b,c,n)=f(a\%c,b\%c,c,n)+\frac{a}{c}\frac{n(n+1)}{2}+(n+1)\frac{b}{c}

a,bca,b\leq c,设m=an+bcm=\frac{an+b}{c},我们开始推式子:
f(a,b,c,n)=i=0nj=0m1[j<ia+bc]f(a,b,c,n)=\sum_{i=0}^n\sum_{j=0}^{m-1}[j<\frac{ia+b}{c}]=j=0m1i=0n[i>cjb+c1a]=\sum_{j=0}^{m-1}\sum_{i=0}^n[i>\frac{cj-b+c-1}{a}]上面那一步怎么来的呢?我们对括号内的内容做推导:
j<ai+bcj<\frac{ai+b}{c}jai+bc1j\leq \frac{ai+b}{c}-1cjai+bccj\leq ai+b-ccj<ai+bc+1cj<ai+b-c+1i>cjb+c1ai>\frac{cj-b+c-1}{a}然后就推出来了,我们继续推导:
=j=0m1ncjb+c1a=\sum_{j=0}^{m-1}n-\frac{cj-b+c-1}{a}=nmf(c,cb1,a,m1)=nm-f(c,c-b-1,a,m-1)发现上述算法过程类似于辗转相除法,故时间复杂度为log\log

咕咕咕
标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!