【BZOJ3817/UOJ42】Sum(类欧)
【BZOJ3817/UOJ42】Sum(类欧)
题面
题解
令\(x=\sqrt r\),那么要求的式子是$$\sum_{d=1}n(-1){[dx]}$$
不难发现,对于每个\(d\)而言的取值只和\([dx]\)的奇偶性相关。
如果\(x\)是个整数,也就是\(r\)是完全平方数的时候,显然是可以直接算答案的。
计算答案的时候显然之和有几个奇数或者几个偶数相关(只要求一个另外一个就是补集)
比如说我们来求有几个是偶数,那么要满足的条件就是\([dx]=2*[\frac{dx}{2}]\)。
先重新写下式子,我们写成这个样子$$\sum_{d=1}^n(1-2*([dx]\ mod\ 2)$$
这个显然成立,当\([dx]\)是偶数的时候贡献是\(1\),奇数的时候贡献是\(0\)。
而\([dx]\ mod\ 2\)可以写成减法的形式。所以原式可以写成
\]
化简之后得到了
\]
现在把模型转化一下,变得更加一般一点,那么要求解的东西是$$\sum _{d=1}^n[\frac{ax+b}{c}d]$$
听说这个玩意叫做类欧?类欧几里得算法。
令\(k=\frac{ax+b}{c}\),那么式子可以化简为\(\sum[kd]\)
这里进行分类讨论
1.当\(k>=1\)的时候
\]
显然后面那一半是可以直接计算的,而前面这一半令\(k'=k-[\frac{ax+b}{c}]\)
我们可以递归处理。
2.\(k<1\)
转化到一个坐标系上面去,我们要求的东西本质上就是有一条直线\(y=kx\),要求解\(x\in[1,n]\)时,与\(x=n\)、\(x\)正半轴围成的三角形内部整点的个数。
我们把这个三角形逆时针旋转\(90\)度,再沿着\(y\)轴翻转过来,让长度为\(n\)的那条边靠着\(y\)轴,这样子翻转过来之后,\(k\)变成了倒数,\(n\)变成了\([nk]\),然后拿矩形减去多出来的部分,不就是矩形减去\(k>=1\)的情况了吗?
再讨论一下取倒数的结果,斜率\(k=\frac{ax+b}{c}\),取倒数之后\(k=\frac{c}{ax+b}\),然后分母有理化一下\(k=\frac{c(ax-b)}{a^2r-b^2}\)。
总的时间复杂度是\(log\)级别的。
#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
#define ll long long
inline int read()
{
int x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return t?-x:x;
}
double x;ll n,r;
ll Calc(ll a,ll b,ll c,ll n)
{
if(!n)return 0;
ll d=__gcd(a,__gcd(b,c));a/=d;b/=d;c/=d;
ll k=(a*x+b)/c;
if(!k)
{
ll N=(a*x+b)/c*n;
return n*N-Calc(a*c,-b*c,a*a*r-b*b,N);
}
else
return k*n*(n+1)/2+Calc(a,b-c*k,c,n);
}
int main()
{
int T=read();
while(T--)
{
n=read();r=read();x=sqrt(r);
ll k=x;
if(k*k==r)
{
if(k&1)printf("%lld\n",n-2*((n+1)/2));
else printf("%lld\n",n);
}
else printf("%lld\n",n+Calc(1,0,2,n)*4-Calc(1,0,1,n)*2);
}
return 0;
}
【BZOJ3817/UOJ42】Sum(类欧)的更多相关文章
- BZOJ3817 清华集训2014 Sum 类欧几里得
传送门 令\(\sqrt r = x\) 考虑将\(-1^{\lfloor d \sqrt r \rfloor}\)魔改一下 它等于\(1-2 \times (\lfloor dx \rfloor \ ...
- UOJ#42. 【清华集训2014】Sum 类欧几里德算法
原文链接https://www.cnblogs.com/zhouzhendong/p/UOJ42.html 题解 首先我们把式子改写一下: $$(-1)^{\lfloor a\rfloor} \\=1 ...
- JZOJ3492数数&&GDOI2018超级异或绵羊——位&&类欧几里得
JZOJ3492 数数(count) 我们知道,一个等差数列可以用三个数A,B,N表示成如下形式: B+A,B+2A,B+3A⋯B+NA ztxz16想知道对于一个给定的等差数列,把其中每一项用二进 ...
- 2019牛客多校九 I. KM and M (类欧几里得)
大意: 给定$N,M$, 求$\sum\limits_{K=1}^N \text{(KM)&M}$ 考虑第$i$位的贡献, 显然为$\lfloor\frac{KM}{2^i}\rfloor$为 ...
- Kattis - itsamodmodmodmodworld It's a Mod, Mod, Mod, Mod World (类欧几里得)
题意:计算$\sum\limits_{i=1}^n[(p{\cdot }i)\bmod{q}]$ 类欧模板题,首先作转化$\sum\limits_{i=1}^n[(p{\cdot}i)\bmod{q} ...
- [ZZOJ#31]类欧几里得
[ZZOJ#31]类欧几里得 试题描述 这是一道模板题. 给出 \(a, b, c, n\),请你求出 \(\sum_{x=0}^n{\lfloor \frac{a \cdot x + b}{c} \ ...
- 算法马拉松35 E 数论只会Gcd - 类欧几里得 - Stern-Brocot Tree - 莫比乌斯反演
题目传送门 传送门 这个官方题解除了讲了个结论,感觉啥都没说,不知道是因为我太菜了,还是因为它真的啥都没说. 如果 $x \geqslant y$,显然 gcd(x, y) 只会被调用一次. 否则考虑 ...
- bzoj2187 fraction&&hdu3637 Find a Fraction——类欧几里得
bzoj2187 多组询问,每次给出 $a, b, c, d$,求满足 $\frac{a}{b} < \frac{p}{q} < \frac{c}{d}$ 的所有二元组 $(p, q)$ ...
- CF1182F Maximum Sine【类欧,扩欧】
题目链接:洛谷 题目描述:求整数$x\in [a,b]$使得$|2px \ mod \ 2q-q|$最小,如果有多个$x$输出最小的. 数据范围:$1\leq a,b,p,q\leq 10^9$ 第一 ...
随机推荐
- 从harbor部署到在k8s中使用
一.概述 harbor是什么呢?英文单词的意思是:港湾.港湾用来存放集装箱(货物的),而docker的由来正是借鉴了集装箱的原理,所以harbor是用于存放docker的镜像,作为镜像仓库使用.官方的 ...
- ubuntu12.04安装mininet
网上安装mininet教程有很多,都是通过git命令安装,但有一个坑,安装到./install.sh时会报错,记录下来 1.通过git 下载mininet git clone git://github ...
- ecna2017-Sheba’s Amoebas
很简单的深搜的一道题,由于这道题要找环的个数,并且认为相连当一个点的8个方向种中有一个方向和这个点相连. 这个题做法无非就是暴力每个点,然后满足条件的深搜即可. 感觉我自己的代码写的很无趣,大佬的代码 ...
- Codeforces Round #505 (rated, Div. 1 + Div. 2, based on VK Cup 2018 Final)-C. Plasticine zebra
问了学长,感觉还是很迷啊,不过懂了个大概,这个翻转操作,实质不就是在序列后面加上前面部分比如 bw | wwbwwbw 操作过后 wbwbwwbww 而 bw | wwbwwbwbw 这样我们就知道 ...
- 洛谷P1004 方格取数-四维DP
题目描述 设有 N \times NN×N 的方格图 (N \le 9)(N≤9) ,我们将其中的某些方格中填入正整数,而其他的方格中则放入数字 00 .如下图所示(见样例): A 0 0 0 0 0 ...
- 1013 B. And
链接 [http://codeforces.com/contest/1013/problem/B] 题意 给你一个n和x,再给n个数,有一种操作用x&a[i]取代,a[i],问使其中至少两个数 ...
- 《面向对象程序设计》第三次作业 Calculator
c++第三次作业 Calculator git上的作业展示点这里. ps:有一点不是很明确,作业要求:将数字和符号提取出来,得到一组string,然后才将这些string存入队列中.按我的理解是需要将 ...
- APP推广(预期方案)
首先,在推广过程中有一些定的弊端:我们这个O2O平台暂时只能适用于学校局域网. 因为我们的APP才刚刚“出炉”不久,在网络上还是属于一篇空白的状态,我们想过可以在百度百科上进行相应的推广,如果有用户搜 ...
- SDN网络虚拟化、资源映射等相关论文粗读
1. Control Plane Latency with SDN Network Hypervisors: The Cost of Virtualization 年份:2016 来源:IEEE NE ...
- Vue的filter属性
Vue.filter('sss',function(data){ //可以使用逻辑处理数据 data += '123'; return data*10; }) new Vue({ el:'#app', ...