2016年06月_ramay7

原创 HDU 2521 反素数（区间反素数）

题目链接： HDU 2521 反素数题意：求出区间[a,b]约数最多的数。求出区间[a,b]约数最多的数。分析：先求[1,b][1,b]中约数最多的数ansans，如果ans∉[a,b]ans\notin [a,b]那就暴力跑。本身这题数据也不大。#include <iostream>#include <cstdio>#include <cstring>#include <s

2016-06-11 23:42:40 584

原创 Codeforeces 27 E Number With The Given Amount Of Divisors(反素数)

题目链接： Codeforeces 27 E Number With The Given Amount Of Divisors 题意：给出nn找到最小的数使得其约数个数恰为nn。#include <iostream>#include <cstdio>#include <cstring>#include <string>#include <algorithm>#include <c

2016-06-11 23:37:59 332

原创 HDU 4542 小明系列故事——未知剩余系(反素数+最小的n使得n的约数为n-k（k已知）)

题目链接： HDU 4542 小明系列故事——未知剩余系#include <iostream>#include <cstdio>#include <cstring>#include <string>#include <algorithm>#include <climits>#include <cmath>#include <ctime>#include <cassert>#def

2016-06-11 23:34:58 470

原创 HDU 3923 Invoker（Polya原理+费马小定理）

题目链接： HDU 3923 Invoker 题意： tt种颜色来涂nn颗珠子的项链，问本质不同的珠子有多少个?考虑翻转和旋转。结果mod1000000007mod 1000000007 分析：模版题。需要注意最后需要除以2n2n,又因为结果%1000000007，可以用费马小定理，相当于乘以quickpow(2∗n,1000000005).quick_pow(2*n, 10000000

2016-06-10 23:26:00 577

原创 UVA 11255 Necklace（每种颜色珠子个数限制、Polya原理、组合数）

题目链接： UVA 11255 Necklace 题意：各有a,b,c(a,b,c≥0,a+b+c≤40)a,b,c(a,b,c \geq 0, a+b+c \leq 40)颗三种颜色，问这些珠子能串成的项链有多少种?考虑翻转和旋转。分析：令∑3i=1color[i]=n\sum_{i=1}^{3}color[i]=n,即珠子总数。考虑旋转置换。我们考虑旋转ii颗珠子的间距，则

2016-06-10 22:38:57 1667 1

原创 SGU 294 He's Circle(Polya原理+高精度)

题目链接： SGU 294 He’s Circle 题意：用两种颜色涂n颗珠子的项链，问本质不同的方案数有多少？考虑旋转。分析：很容易根据PolyaPolya原理得到：Ans=1n∑d|nn2dϕ(nd)Ans=\frac{1}{n}\sum_{d|n}^{n}2^d\phi(\frac{n}{d}) 但是因为n≤200000n \leq 200000,并且答案不需要模数。那么难点

2016-06-10 17:23:04 1038

原创 POJ 2154 Color(Polya原理+欧拉函数)

题目链接： POJ 2154 Color 题意： nn个颜色涂nn个珠子的项链，考虑旋转，求本质不同的项链数？，分析：根据PolyaPolya原理可得Ans=∑i=1nngcd(n,i)=∑d=1d=nnd∑i=1i=n[gcd(n,i)==d]=∑d=1nnd∑i=1n[gcd(nd,id)==1]=∑d=1nnd∗ϕ(nd) Ans

2016-06-09 21:37:34 545

原创 HDU 1812 Count the Tetris(Polya原理+高精度)

题目链接： HDU 1812 Count the Tetris 题意：有一个n*n的网格，用m种颜色来涂，求本质不同的方案数？分析； Polya原理+高精度。从原题中描述可以看出要考虑旋转，但是翻转的情况好像没有说清吧。。。还有需要记录已经算出来的答案，当再次询问时，直接输出，否则会TLE，这也太恶心了吧。。。#include <cstdio>#include <cstrin

2016-06-09 20:50:07 707

原创 POJ 1286 Necklace of Beads(Polya原理)

题目链接： POJ 1286 Necklace of Beads 题意：有3种颜色来涂n颗珠子的项链，考虑翻转和旋转，问不同的项链个数?有3种颜色来涂n颗珠子的项链，考虑翻转和旋转，问不同的项链个数? 分析：假设有tt种颜色，nn颗珠子。旋转如果逆时针旋转ii颗珠子的间距，则珠子0,i,2i,...0,i,2i,...构成一个循环，这个循环有i∗ngcd(i,n\frac{i*n}

2016-06-08 23:40:37 475

原创 ZOJ 3435 Ideal Puzzle Bobble(gcd(i,j,k)=1/莫比乌斯反演)

题目链接： ZOJ 3435 Ideal Puzzle Bobble 题意： ∑i=0i=a∑j=0j=b∑k=0k=c[gcd(i,j,k)==1],a,b,c∈[1,1000000]\sum_{i=0}^{i=a}\sum_{j=0}^{j=b}\sum_{k=0}^{k=c}[gcd(i,j,k)==1],a,b,c\in [1,1000000] 分析； 1.1.当i=j=k=0i=

2016-06-06 21:57:05 1970

原创 BZOJ 2154 Crash的数字表格（sigma(lcm(i,j))，莫比乌斯反演）

题目链接；BZOJ 2154 Crash的数字表格题意： ∑i=1i=n∑j=1j=mlcm(i,j) % 100000009 (n,m≤107)\sum_{i=1}^{i=n} \sum_{j=1}^{j=m}lcm(i,j)\ \%\ 100000009\ \ (n,m\leq 10^{7}) 分析： ans=∑d=1d=n∑i=1i=n∑j=1j=mi∗jd (gcd(i,j)=

2016-06-06 21:52:10 799

原创 POJ 2773 Happy 2006(求第k个和m互素的数/欧拉函数)

题目链接： POJ 2773 Happy 2006 题意：求第kk个和mm互素的数。m<=106,k<=108m<=10^6,k<=10^8 分析：如果gcd(a,m)=1gcd(a,m)=1,那么gcd(a+k∗m,m)=1,k∈Zgcd(a+k*m,m)=1,k\in Z证明:只需要证明gcd(a+m,m)=1gcd(a+m, m)=1.令d=gcd(a+m,m)d=gcd(a+m

2016-06-06 21:39:46 574

原创 Codeforces 235 E Number Challenge(莫比乌斯反演)

题目链接： Codeforces 235 E Number Challenge 题意: 记d(i)表示i的约数个数，计算：∑i=1a∑j=1b∑k=1cd(ijk)记d(i)表示i的约数个数，计算：\sum_{i=1}^{a}\sum_{j=1}^{b}\sum_{k=1}^{c}d(ijk) 分析： Ans = ∑i=1a∑j=1b∑k=1cd(ijk) = ∑i=1a⌊ai⌋∑j=1b

2016-06-06 16:00:29 1110

原创 BZOJ 3994 约数个数和(莫比乌斯反演)

题目链接： BZOJ 3994 约数个数和题意：求∑ni=1∑mj=1d(i∗j),定义d(i)为i的约数个数.n,m∈[1,50000]求\sum_{i=1}^{n}\sum_{j=1}^{m}d(i*j),定义d(i)为i的约数个数.n,m\in [1, 50000] 分析： ans=∑gcd(i,j)=1⌊ni⌋⌊mj⌋=∑i=1n⌊ni⌋∑j=1m⌊mj⌋ans =\sum_{

2016-06-06 02:25:26 839

原创 HDU 4746 Mophues(有趣的前缀和/莫比乌斯反演)

题目链接： HDU 4746 Mophues 题意：令gcd(a,b)=d，对d质因子分解得到质因子个数为k如果k≤C，则称k为C的lucky数字，给出n,m,p，求a∈[1,n],b∈[1,m]使得gcd(a,b)是p的lucky数字的数字对(a,b)有多少对？令gcd(a, b) = d，对d质因子分解得到质因子个数为k如果k\leq C，则称k为C的lucky数字，给出n,m,p，求a

2016-06-05 21:56:31 804

原创 BZOJ 2818 Gcd（gcd(x,y)为素数/欧拉函数/莫比乌斯反演）

题目链接： BZOJ 2818 Gcd 题意： x∈[1,N],y∈[1,N]，gcd(x,y)=素数的有序对(x,y)的对数。x\in [1,N],y\in [1, N]，gcd(x,y)=素数的有序对(x,y)的对数。分析：对于一个素数p，如果gcd(x,y)=p，那么相当于x∈[1,Np],y∈[1,Ny]的(x,y)的对数，又因为是有序对，需要乘以2，那么就是∑Npi=12∗ϕ

2016-06-04 14:55:04 1528

原创 HDU 1695 GCD（gcd(x,y)=k无序对数/容斥原理）

题目链接： HDU 1695 GCD 题意: 求x∈[a,b],y∈[c,d]，且gcd(x,y)=k的(x,y)的无序对对数。a=c=1x\in [a,b],y\in[c,d]，且gcd(x,y)=k的(x,y)的无序对对数。a=c=1 分析：和BZOJ 2301的区别就是这里是无序对。不妨设b≤d,,利用容斥原理那么多计算的部分就是x∈[1,b],y∈[1,b]的部分，那么减掉这

2016-06-04 14:13:57 1059

原创 BZOJ 2301 Problem B（x属于[a,b],y属于[c,d]满足gcd(x,y)=k的(x,y)的有序对数）

题目链接： BZOJ 2301 Problem B 题意: 区间x∈[a,b],y∈[c,d]满足gcd(x,y)=k的(x,y)的对数x\in [a,b],y \in[c,d]满足gcd(x,y)=k的(x,y)的对数分析：莫比乌斯反演+优化。#include <cstdio>#include <cmath>#include <cstring>#include <string>

2016-06-04 14:01:53 1285

原创 SPOJ - VLATTICE Visible Lattice Points（gcd(x,y,z)=1的对数/莫比乌斯反演）

题目链接： SPOJ - VLATTICE Visible Lattice Points 题意：一个n*n*n的方格，从最左下角(0, 0, 0)最多可以看到多少个点？（不被遮挡）包括方格内部。分析：假设能看到的点的坐标为(x,y,z)则必须满足：gcd(x,y,z)=1。(0≤x,y,z≤n)。当x=y=z=0时是不成立的。当x,y,z中有两个为0时,只有三种情况(0,0,1),(

2016-06-04 13:57:32 820

原创 BZOJ 2005 能量采集（莫比乌斯反演）

题目链接： BZOJ 2005 能量采集题意：一块n*m的土地，能量采集器位与(0, 0)，如果一棵植物与能量采集器连接成的线段上有k棵植物，那么能量损失为2*k+1,如果没有植物能量损失为1.求总的能量损失。分析：定义f(d)为gcd(x,y)=d(x∈[1,n],y∈[1,n]的(x,y)对数，则：ans=∑i=ni=1(2∗i−1)∗f(i)定义f(d) 为gcd(x, y)

2016-06-04 13:50:40 512

ramay7