Pentium.Labs

hdu 4609 FFT

摘要：题意：给出一堆数，问从这些数中取3个能组成三角形的概率？sol：其实就是问从这些数里取3个组成三角形有多少种取法脑洞大开的解法：用FFT设一开始的数是1 3 3 4作一个向量x，其中x[i]=边长为i的边的个数那么就有x=[0 1 0 2 1 0 0 0 0]令y=x，对x和y作DFT，得到dx和d... 阅读全文

posted @ 2015-04-28 18:56 Pentium.Labs 阅读(295) 评论(0) 推荐(0) 编辑

hdu1402 FFT入门

摘要：参考这里：http://www.cnblogs.com/pdev/p/4354705.html http://www.cnblogs.com/pdev/p/4354629.html题意：求大数乘法A*BA和B位数很长。裸高精度时间复杂度是O(nm)，会完蛋不妨回忆下裸高精度的过程：其实乘法的... 阅读全文

posted @ 2015-04-27 21:47 Pentium.Labs 阅读(678) 评论(0) 推荐(0) 编辑

多项式乘法快速算法

摘要：多项式乘法优化算法：设有如下两个多项式：把它们的系数分别做成向量X=[x0,x1,x2,x3,......]的形式，得 F=[2,3,1] G=[5,2,0]那么根据卷积公式可以求得向量F和G的卷积S=[10,19,11,2]而由多项式乘法可算出各项系数和上面的卷积结果正好一一对应。所以说，多项... 阅读全文

posted @ 2015-03-20 21:58 Pentium.Labs 阅读(1196) 评论(0) 推荐(0) 编辑

FFT

摘要： Fast Fourier Transform首先介绍几个概念：1.卷积：设f(x)和g(x)是两个可积函数，令那么h(x)就是f(x)和g(x)的卷积但更常用的是离散的情况：卷积的意义：如下图首先将两个函数都用T来表示。对其中一个函数做水平翻转：g(T)→g(... 阅读全文

posted @ 2015-03-20 21:24 Pentium.Labs 阅读(286) 评论(0) 推荐(0) 编辑

GDUT校赛

摘要：题目链接：http://4.gdutcode.sinaapp.com/contest.php?cid=1021F题意：给出n和m，要求满足gcd(x,y)=n && lcm(x,y)=m的pair(x,y)的个数sol：先YY一下：设 gcd(x,y)=n，lcm(x,y)=m那么有x=a*n，y... 阅读全文

posted @ 2015-03-19 21:25 Pentium.Labs 阅读(325) 评论(0) 推荐(0) 编辑

light oj 1236 分解质因数

摘要：题目链接：http://acm.hust.edu.cn/vjudge/contest/view.action?cid=70017#problem/H题意：求满足1的数目一开始我是这么想的：既然lcm(i,j)=n，那么n=x*i=y*j，且x和y一定互质。若i和j固定了，那么x和y也固定了。那么问题... 阅读全文

posted @ 2015-03-16 01:38 Pentium.Labs 阅读(735) 评论(0) 推荐(0) 编辑

hdu5187 奇怪题

摘要：本来很水的，答案就是(2^n)-2，但是写坑了QAQ因为原题要求答案要mod P，一开始我是这么干的： LL ans=pow_mod(2,N,P); ans=(ans-2)%P; if (N==1) ans=1%P; printf("%I6... 阅读全文

posted @ 2015-03-15 11:56 Pentium.Labs 阅读(261) 评论(0) 推荐(0) 编辑

hdu5183 hash大法

摘要：维护前缀和sum[i]=a[0]-a[1]+a[2]-a[3]+…+(-1)^i*a[i]枚举结尾i，然后在hash表中查询是否存在sum[i]-K的值。如果当前i为奇数，则将sum[i]插入到hash表中。上面考虑的是从i为偶数为开头的情况。然后再考虑以奇数开头的情况，按照上述方法再做一次即可。不... 阅读全文

posted @ 2015-03-08 18:27 Pentium.Labs 阅读(211) 评论(0) 推荐(0) 编辑

poj3696 快速幂的优化+欧拉函数+gcd的优化+互质

摘要：这题满满的黑科技orz题意：给出L，要求求出最小的全部由8组成的数(eg: 8,88,888,8888,88888,.......)，且这个数是L的倍数sol：全部由8组成的数可以这样表示：((10^x)-1)*(8/9)那么有m=((10^x)-1)*(8/9)=k*L，answer即满足条件的最... 阅读全文

posted @ 2015-03-05 19:28 Pentium.Labs 阅读(299) 评论(0) 推荐(0) 编辑

hdu5175 gcd 求约数

摘要：题意：求满足条件GCD(N,M) = N XOR M的M的个数sol：和uva那题挺像的。若gcd(a,b)=a xor b=c，则b=a-c暴力枚举N的所有约数K，令M=NxorK,再判断gcd(N,M)是不是等于K。注意枚举约数时传统方法是O(N)的，会完蛋有个O(sqrt(N))的方法：注意一... 阅读全文

posted @ 2015-03-03 21:32 Pentium.Labs 阅读(159) 评论(0) 推荐(0) 编辑

uva12716 gcd

摘要：题意：给出N，1 2 #include 3 #define LL long long 4 #define MX 30000005 5 6 int N,T; 7 LL S[MX],f[MX]; 8 9 int main()10 {11 memset(f,0,sizeof(f));12 ... 阅读全文

posted @ 2015-03-02 18:05 Pentium.Labs 阅读(223) 评论(0) 推荐(0) 编辑

uva11426 gcd、欧拉函数

摘要：题意：给出N，求所有满足i gcd(x/i,n/i)=1 -> x/i与n/i互质 -> 满足该条件的x/i有phi(n/i)个那么再用欧拉函数就可以求出每一个f(n)啦~如果找n的每一个约数i会有点慢，可以枚举i，令n=2*i，3*i，........（n是i的所有倍数且小于MAXN） for (... 阅读全文

posted @ 2015-03-02 17:27 Pentium.Labs 阅读(144) 评论(0) 推荐(0) 编辑

hdu3982 直线切多边形【WA中...】

摘要：题意：有一块蛋糕，上面有一颗cherry。用刀子切n次，求切完之后有cherry的那部分的面积My solution：先做一个大矩形，使cake内切于这个大矩形。如图：然后不断切这个大矩形，每次切割的时候保留与cherry同侧的那部分。最后剩下的就是一个多边形。求该多边形与圆的面积交即可。在切割的时... 阅读全文

posted @ 2015-02-11 20:08 Pentium.Labs 阅读(656) 评论(0) 推荐(0) 编辑

bc#29 做题笔记

摘要：昨天的bc被坑惨了= =本来能涨rating的大好机会又浪费了。。。大号已弃号A：第一反应是高精度，结果模板找不到了= =，然后现学现卖拍了个java的BigInteger+快速幂，调了好半天不说还TLE。貌似这题就在卡java实际上尼玛等号两边取log不就完了么。。。卒B：A题调了半天，开始做B的... 阅读全文

posted @ 2015-02-08 16:55 Pentium.Labs 阅读(184) 评论(0) 推荐(0) 编辑

poj3675 求多边形与圆的面积交

摘要：题意：给出多边形的顶点坐标、圆的圆心坐标和半径，求面积交sol：又是模板题啦= =注意poj的C++好像认不出hypot函数，要稍微改写一下。hypot(double x,double y)：即返回sqrt(x*x+y*y)的值 1 #include 2 #include 3 #include... 阅读全文

posted @ 2015-02-07 17:42 Pentium.Labs 阅读(365) 评论(0) 推荐(0) 编辑

poj1279 半平面交

摘要：题意：没看懂= =sol：在纸上随便画两下就可以看出，答案即按逆时针方向建立line，求它们的半平面交的面积。模板题。注意输出答案时输出ans+eps，否则可能会出现结果为-0.00的情况。 1 #include 2 #include 3 #include 4 #include 5 #in... 阅读全文

posted @ 2015-02-06 18:05 Pentium.Labs 阅读(192) 评论(0) 推荐(0) 编辑

poj3335 半平面交

摘要：题意：给出一多边形。判断多边形是否存在一点，使得多边形边界上的所有点都能看见该点。sol：在纸上随手画画就可以找出规律：按逆时针顺序连接所有点。然后找出这些line的半平面交。题中给出的点已经按顺时针排好序了，所以只要倒过来一下就可以了。很简单的模板题。 1 #include 2 #includ... 阅读全文

posted @ 2015-02-06 17:29 Pentium.Labs 阅读(170) 评论(0) 推荐(0) 编辑

POJ3525 半平面交

摘要：题意：求某凸多边形内部离边界最远的点到边界的距离首先介绍半平面、半平面交的概念：半平面：对于一条有向直线，它的方向的左手侧就是它所划定的半平面范围。如图所示：半平面交：多个半平面的交集。有点类似二元函数的线性规划。如图求半平面交：用的kuangbin模板= =sol：二分答案二分距离值，按这个值... 阅读全文

posted @ 2015-02-05 20:47 Pentium.Labs 阅读(643) 评论(0) 推荐(0) 编辑

poj2187 旋转卡(qia)壳(ke)

摘要：题意：求凸包的直径关于对踵点对、旋转卡壳算法的介绍可以参考这里：http://www.cnblogs.com/Booble/archive/2011/04/03/2004865.htmlhttp://www.cppblog.com/staryjy/archive/2009/11/19/101412.... 阅读全文

posted @ 2015-02-01 17:24 Pentium.Labs 阅读(237) 评论(0) 推荐(0) 编辑

poj2284 欧拉公式

摘要：题意：给出一图形，求该图形把平面分成了几部分欧拉公式： http://blog.csdn.net/wangxiaojun911/article/details/4586550对于二维平面上的情况。设图形上有V个点，E条边，把平面分成了F个独立的部分，那么满足V+F-E=2如下图：那么求... 阅读全文

posted @ 2015-01-29 16:26 Pentium.Labs 阅读(375) 评论(2) 推荐(0) 编辑

导航

公告