BZOJ2831(小强的金字塔系列问题--区域整点数求法)
题意就是给出A,B,C,R,L,然后求
这里其实用到扩展欧几里德。(基本上参照clj的解题报告才理解的)
分析:我们先来分析一般情况:
这里我们假设A<C和B<C,否则我们可以把它化成A<C,B<C的情况
我们令:,所以上式就等价于:
,如果
,那么
的值是1,否则为0
然后我们交换顺序:,由于
,所以
推导过程:因为,所以
,那么
,进一步有:
,所以
令:
那么式子就转化为:
这样一看我们就交换了A,C了,也就是扩展欧几里德算法,详见金斌2009论文。
我们继续化简,发现又把问题转化成求: 和
这样的话与上面同样的思路继续化简:对于比较简单,就不写过程了。
下面来说说的化简过程:
然后分析就完毕!
由于本题的数很大,要用高精度,用Python最好,因为一是Python高精度可以直接算,二是Python直接返回元组方便。
sum = lambda n : n*(n+1)/2
sqrtsum = lambda n : n*(n+1)*(2*n+1)/6 def f(a,b,c,r):
if r<0:
return [0]*3 if a>=c:
t=f(a%c,b,c,r)
a/=c
return [t[0]+a*sqrtsum(r),t[1]+a*sum(r),t[2]+a*a*sqrtsum(r)+2*a*t[0]] elif b>=c:
t=f(a,b%c,c,r)
b/=c
return [t[0]+b*sum(r),t[1]+b*(r+1),t[2]+b*b*(r+1)+2*b*t[1]] else:
if a==0:
return [0]*3
y=(a*r+b)/c
t=f(c,c-b-1,a,y-1)
ans=[y*sum(r)-(t[1]+t[2])/2,y*r-t[1],y*y*r-2*t[0]-t[1]]
return ans a,c,b,l,r=map(int,raw_input().split())
print f(a,b,c,r)[0]-f(a,b,c,l-1)[0]
典型题目一:SPOJ4717. Grid Points in a Triangle
题意:求满足 y <= ax / b and x <= n的整点个数,x,y,a,b,n都是非负的整数。
分析:实际上就是求:,就是上题的思路。
#include <iostream>
#include <stdio.h>
#include <string.h> using namespace std;
typedef long long LL; LL gcd(LL a,LL b)
{
return b? gcd(b,a%b):a;
} LL dfs(LL n,LL a,LL b)
{
LL t=n*(n+1)/2*(a/b);
a%=b;
if(a==0) return n+1+t;
LL d=a*n/b;
t+=(n+1)*(d+1)+d/a+1;
return t-dfs(d,b,a);
} int main()
{
LL t,n,a,b,p;
scanf("%lld",&t);
while(t--)
{
scanf("%lld%lld%lld",&n,&a,&b);
p=gcd(a,b);
printf("%lld\n",dfs(n,a/p,b/p));
}
return 0;
}
典型题目二:2013年ACM全国邀请赛南京赛区一题:http://icpc.njust.edu.cn/Local/1742
分析:本题比较难,其实方法跟上题差不多,要分析的过程更加复杂些。
BZOJ2831(小强的金字塔系列问题--区域整点数求法)的更多相关文章
- POJ 1265 /// 皮克定理+多边形边上整点数+多边形面积
题目大意: 默认从零点开始 给定n次x y上的移动距离 组成一个n边形(可能为凹多边形) 输出其 内部整点数 边上整点数 面积 皮克定理 多边形面积s = 其内部整点in + 其边上整点li / 2 ...
- .NET Core 3 WPF MVVM框架 Prism系列之区域管理器
本文将介绍如何在.NET Core3环境下使用MVVM框架Prism的使用区域管理器对于View的管理 一.区域管理器 我们在之前的Prism系列构建了一个标准式Prism项目,这篇文章将会讲解之前项 ...
- MVVM架构~knockoutjs系列之文本框数符长度动态统计功能
返回目录 这个功能为什么要写呢,因为在之前做了一个前端的页面效果,使用JS写的,感觉很累,真的,对于一个文本框长度动态统计,你要写blur,press,down什么的事件,太麻烦了,这时,我想到了kn ...
- 剑指offer系列59---寻找丑数
[题目]把只包含因子2.3和5的数称作丑数(Ugly Number). * 例如6.8都是丑数,但14不是,因为它包含因子7. 习惯上我们把1当做是第一个丑数.求按从小到大的顺序的第N个丑数. 解法一 ...
- COJ 0985 WZJ的数据结构(负十五)(限定区域不同数)
传送门:http://oj.cnuschool.org.cn/oj/home/addSolution.htm?problemID=955 试题描述: CHX有一个问题想问问大家.给你一个长度为N的数列 ...
- CCF系列之有趣的数(201312-4)
题目链接: http://115.28.138.223:81/view.page?opid=4 试题名称: 有趣的数 时间限制: 1.0s 内存限制: 256.0MB 问题描述: 问题描述 我们把一个 ...
- Codeforces 1036E Covered Points (线段覆盖的整点数)【计算几何】
<题目链接> <转载于 >>> > 题目大意: 在二维平面上给出n条不共线的线段(线段端点是整数),问这些线段总共覆盖到了多少个整数点. 解题分析: 用GC ...
- C#进阶系列——WebApi 接口返回值不困惑:返回值类型详解
前言:已经有一个月没写点什么了,感觉心里空落落的.今天再来篇干货,想要学习Webapi的园友们速速动起来,跟着博主一起来学习吧.之前分享过一篇 C#进阶系列——WebApi接口传参不再困惑:传参详解 ...
- C#进阶系列——WebApi 接口参数不再困惑:传参详解
前言:还记得刚使用WebApi那会儿,被它的传参机制折腾了好久,查阅了半天资料.如今,使用WebApi也有段时间了,今天就记录下API接口传参的一些方式方法,算是一个笔记,也希望能帮初学者少走弯路.本 ...
随机推荐
- align="absmiddle" 的意义
align=absmiddle表示图像的中间与同一行中最大元素的中间对齐 AbsBottom 图像的下边缘与同一行中最大元素的下边缘对齐. AbsMiddle 图像的中间与同一行中最大元素的中间对 ...
- java基础37 集合框架工具类Collections和数组操作工具类Arrays
一.集合框架工具类:Collections 1.1.Collections类的特点 该工具类中所有的方法都是静态的 1.2.Collections类的常用方法 binarySearch(List< ...
- set,list,map分析
想了下集合,列表,映射三者关系 set,list,map ArrayList采用数组方式存储数据,继承List; LinkedList采用链表方式存储数据,继承List; 所以数组方式都有下表,以及每 ...
- python之uinttest,用例执行顺序
unittest单元测试框架, 以test开头的测试用例,默认执行顺序是按照ASC码来执行 如果有类,先排序执行类,在执行类中,再排序用例顺序执行 如果想要按照指定的顺序执行测试用例. 那么就需要用到 ...
- 洛谷P1186玛丽卡
传送门啦 先跑一遍最短路,将最短路的路径记录下来,然后枚举每一条最短路的边,将其断掉,记录此时的1-n的时间,取其中最大的一个时间即为所求. (通过 $ cut[][] $ 和 $ f[] $ 进行操 ...
- algorithm之unique
#include<iostream> #include<algorithm> usingnamespace std; /**< unique函数的算法思想 */ vect ...
- springMVC源码分析--HttpMessageConverter写write操作(三)
上一篇博客springMVC源码分析--HttpMessageConverter参数read操作中我们已经简单介绍了参数值转换的read操作,接下来我们介绍一下返回值的处理操作.同样返回值的操作操作也 ...
- sql server中扩展存储过程
--列出服务器上安装的所有OLEDB提供的程序 execute master..xp_enum_oledb_providers --得到硬盘文件信息 --参数说明:目录名,目录深度,是否显示文件 (少 ...
- 容器计划任务大坑:在alpine容器里,想用非root帐号执行crontab任务
我只能说抱歉,我前前后后测试了七天, 将自己预想的配置错误,一个一个去验证. 非root帐号在alpine容器里执行crontab任务,还是失败, 输出依旧是一片空白~ stackoverflow里, ...
- Python类和人类
Python中的类 俗话说,物以类聚,人以群分,类是什么,类就是一组相同属性的集合.下面来结合人,探讨Python中类和人类的关系. 首先,我们定义一个人的类,如下: class People(obj ...