数学基础
文章平均质量分 62
ZMOIYNLP
嗨嗨嗨嗨
展开
-
为什么线性筛欧拉函数i%prime[j]==0的时候phi[i*prime[j]]=phi[i]*prime[j]
看贾志鹏线性筛的时候想起来的。 我有一个繁琐的证明- -。 证明ϕ(pm)=p×ϕ(m),p为素数,m∈Z\phi(pm)=p\times\phi(m),p为素数,m\in \Bbb Z. 设m=pα⋅m′,α,m′∈N,(pα,m′)=1.m=p^{\alpha}\cdot m',\;\alpha,m'\in\Bbb N,(p^{\alpha},m')=1. 那么ϕ(m)=ϕ(m′)⋅ϕ(原创 2015-03-23 17:34:40 · 1964 阅读 · 0 评论 -
【UVa12298】 Super Joker II 【FFT】【生成函数】
居然因为精度问题WA掉了。。。 这题要用long double,并且交c++11…… 每种花色的牌开一个系数向量,有这种牌那么系数就是1,否则为0.然后FFT乘起来就行了。#include<cstdio> #include<cstring> #include<iostream> #include<cmath> using namespace std; typedef long double re原创 2015-04-14 19:51:59 · 2149 阅读 · 0 评论 -
一点数论题目
今年省选真是悲催第二天第二题没有人A,导致一堆二百五(我是说分数二百五,没有其他的意思~~)。。 bzoj1101:求∑ai=1∑bj=1[gcd(i,j)==d]\sum_{i=1}^a\sum_{j=1}^b[gcd(i,j)==d]。 bzoj2005:等价于求∑ai=1∑bj=1(2gcd(i,j)−1)\sum_{i=1}^a\sum_{j=1}^b \left( 2gcd(i,j)-原创 2015-04-21 19:53:48 · 1660 阅读 · 0 评论 -
【bzoj1101】Zap【神奇的∑】
传送门: http://server.mclscloud.com:5230/JudgeOnline/problem.php?id=1101 求∑ai=1∑bj=1[gcd(i,j)==d]\sum_{i=1}^a\sum_{j=1}^b[gcd(i,j)==d]. 设a≤ba \le b. 令a′=⌊ad⌋,b′=⌊bd⌋a'=\lfloor \frac ad \rfloor,b'=\lfl原创 2015-04-21 20:16:32 · 2011 阅读 · 0 评论 -
【bzoj2005】能量采集【GCD】
为么很多这种题把∑\sum顺序换一下就得到答案了。。。 一个植物(坐标(x,y))到原点的路线上经过的植物数是gcd(x,y)(包括那个植物本身) ∑((gcd−1)∗2+1)=∑(gcd∗2−1)\sum\left((gcd-1)*2+1\right)=\sum(gcd*2-1) 因此把gcd∗2−1gcd*2-1的和求出来即可。 下面说一说如何快速求∑ai=1∑bj=1gcd(i,j)\原创 2015-04-21 20:44:49 · 1717 阅读 · 0 评论 -
【bzoj1004】Cards【Polya计数定理】【递推】
传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=1004 这是一道Polya好题~ 根据那个什么引理,本质不同的方案数等于每个置换下不同的方案数的平均值。 但是Polya定理是: l=1|G|∑f∈Gkm(f)l=\frac 1{|G|} \sum_{f \in G}{k^{m(f)}} 而这里是有三种颜色,求每个置换里每个循环涂同样原创 2015-03-22 20:48:02 · 948 阅读 · 0 评论