【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$ 第一 ...
随机推荐
- A2dp初始化流程源码分析
蓝牙启动的时候,会涉及到各个profile 的启动.这篇文章分析一下,蓝牙中a2dp profile的初始化流程. 我们从AdapterState.java中对于USER_TURN_ON 消息的处理说 ...
- 分布式全文搜索引擎ElasticSearch
一 什么是 ElasticSearch Elasticsearch 是一个分布式可扩展的实时搜索和分析引擎,一个建立在全文搜索引擎 Apache Lucene(TM) 基础上的搜索引擎.当然 Elas ...
- item 12: 把重写函数声明为“override”的
本文翻译自modern effective C++,由于水平有限,故无法保证翻译完全正确,欢迎指出错误.谢谢! 博客已经迁移到这里啦 C++中的面向对象编程总是围绕着类,继承,以及虚函数.这个世界中, ...
- 2018年高教社杯全国大学生数学建模竞赛A题解题思路
题目 先贴一下A的题目吧 A题 高温作业专用服装设计 在高温环境下工作时,人们需要穿着专用服装以避免灼伤.专用服装通常由三层织物材料构成,记为I.II.III层,其中I层与外界环境接触,III层与 ...
- [2019BUAA软件工程]第1次阅读作业
[2019BUAA软件工程]第1次阅读作业 Tips Link 作业连接 [2019BUAA软件工程]第1次阅读作业 读<构建之法>的疑惑 个人开发流程(Personal Software ...
- phpcms全站搜索
这篇博客已经移至http://www.cnblogs.com/nuanai/p/8028562.html中~~~~~~
- 接口(interface)与多态
1. 接口(interface)是抽象方法与常量值的集合: 2. 从本质上来讲,接口是一种特殊的抽象类,这种抽象类中只包含常量与方法的定义,而没有变量和方法的实现: 3. 接口中声明的属性默认为:pu ...
- 正则表达式(java)
概念: 正则表达式,又称规则表达式.(英语:Regular Expression,在代码中常简写为regex.regexp或RE),计算机科学的一个概念. 正则表通常被用来检索.替换那些符合某个模式( ...
- java.lang.Exception: No tests found matching Method tes(com.bw.test.Testrefiect) from org.junit.vintage.engine.descriptor.RunnerRequest@3bfdc050 at org.junit.internal.requests.FilterRequest.getRunner
junit 方法 没有加上注解 @Test java.lang.Exception: No tests found matching Method tes(com.bw.test.Testre ...
- 使用phpstudy创建本地虚拟主机
在使用php开发网站的时候,每次测试自己的网站时,通常都是用localhost/dirname/filename.php来访问自己所写的程序 但是有时候我们需要模拟真实的场景,如通过域名访问时,如果你 ...