首页 > 其他分享 >类欧几里德

类欧几里德

时间:2024-09-17 16:45:35浏览次数:1  
标签:lfloor frac ai 欧几里德 rfloor sum

1.用途

在\(O(logn)\)的时间复杂度内求出\(\sum_{i=0}^{n}\lfloor\frac{ai+b}{c}\rfloor\)

2.推导

1.当a>c且b>c时

\(\sum_{i=0}^{n}\lfloor\frac{ai+b}{c}\rfloor\)
\(=\sum_{i=0}^{n}\lfloor\frac{(\lfloor\frac{a}{c}\rfloor*c+a\%c)i+ \lfloor\frac{b}{c}\rfloor*c+b\%c}{c}\rfloor\)
\(=\sum_{i=0}^{n}\lfloor \lfloor\frac{a}{c}\rfloor*i+ \lfloor\frac{b}{c}\rfloor+\frac{a\%c*i+b\%c}{c}\rfloor\)
\(=\sum_{i=0}^{n} \lfloor\frac{a}{c}\rfloor*i+ \lfloor\frac{b}{c}\rfloor+\lfloor\frac{a\%c*i+b\%c}{c}\rfloor\)
\(= \lfloor\frac{a}{c}\rfloor*\frac{n*(n+1)}{2}+ \lfloor\frac{b}{c}\rfloor*(n+1)+\sum_{i=0}^{n}\lfloor\frac{a\%c*i+b\%c}{c}\rfloor\)

2.当a<c且b<c时

\(\sum_{i=0}^{n}\lfloor\frac{ai+b}{c}\rfloor\)
\(=\sum_{i=0}^{n}\sum_{j=0}^{\lfloor\frac{ai+b}{c}\rfloor-1}1\)
\(= \sum_{j=0}^{\lfloor\frac{an+b}{c}\rfloor-1}\sum_{i=0}^{n}[j<= \lfloor\frac{ai+b}{c}\rfloor-1]\)
\(= \sum_{j=0}^{\lfloor\frac{an+b}{c}\rfloor-1}\sum_{i=0}^{n}[j< \frac{ai+b}{c}]\)
\(= \sum_{j=0}^{\lfloor\frac{an+b}{c}\rfloor-1}\sum_{i=0}^{n}[\lfloor\frac{j*c-b}{a}\rfloor< i]\)
$= \sum_{j=0}^{\lfloor\frac{an+b}{c}\rfloor-1}n-\lfloor\frac{j*c-b}{a}\rfloor $

$= n* \lfloor\frac{an+b}{c}\rfloor-\sum_{j=0}^{\lfloor\frac{an+b}{c}\rfloor-1}\lfloor\frac{c*j-b}{a}\rfloor $

标签:lfloor,frac,ai,欧几里德,rfloor,sum
From: https://www.cnblogs.com/Owen1234/p/18416955

相关文章

  • 推式子——拓展欧几里德算法exgcd
    试求方程\(ax+by=\gcd(a,b)\)的一组整数解,其中\(a\)和\(b\)是给定的数提前准备好一个式子:辗转相除法\[\gcd(a,b)=\gcd(b,a\bmodb)\]同时我们可以注意到,\(\bmod\)的等价版本:\[a\bmodb=a-b\lfloor\frac{a}{b}\rfloor\]开始推式子首先考虑\(b=0\)的情况,因为......
  • dp-双调欧几里德旅行商问题
    双调欧几里德旅行商问题目录双调欧几里德旅行商问题问题描述分析递推关系程序算法导论3rd-15.3问题描述平面上n个点,确定一条连接各点的最短闭合旅程。这个解的一般形式为NP的(在多项式时间内可以求出)J.L.Bentley建议通过只考虑双调旅程(bitonictour)来简化问题,这种旅程即......
  • 关于扩展欧几里德算法的一点小思考
    求ax=c%b的最小正整数解首先我们应该知道所有的同余方程,均是从不定方程来的例如ax=c%b,就是将不定方程ax+by=c的左右两边,同时模b得来的于是我们先试着对ax+by=c这个方程进......
  • 洛谷P1290 欧几里德的游戏
    链接:https://www.luogu.com.cn/problem/P1290不妨假设\(b\leqa\)。显然,当\(a\)是\(b\)的倍数时,为必胜态。接下来考虑\(a\)不为\(b\)的倍数时:1.\(a\)小于\(2*b\)时,当前......