2017年04月_孤鸿子_

原创 2016江苏省赛H题(莫比乌斯反演+分块求和)

题目描述就是让你计算这个东西 ∑i=1n∑j=1mi2∗j2gcd(i,j)n,m≤1e6\sum_{i = 1}^n\sum_{j=1}^mi^2*j^2gcd(i,j)\\n,m\le 1e6分析和经典的这个公式很像∑ni=1∑mj=1gcd(i,j)\sum_{i = 1}^n\sum_{j=1}^mgcd(i,j)只需做一些简单的转化. ∑i=1n∑j=1mi2∗j2gcd(i,j

2017-04-29 13:03:12 1056

原创 10692 - Huge Mods(指数循环节)

题目链接10692 - Huge Mods分析其实这个问题很简单,有以下定理可用若(a,n)=1,ax≡ax mod ϕ(m)(mod m)(a,n) = 1,a^x \equiv a ^{x\ mod\ \phi(m)}(mod\ m)若(a,n)≠1,ax≡ax mod ϕ(m)+ϕ(m)(mod ϕ(m))(a,n) \neq 1,a^x\equiv a^{x\ mod\ \phi(m)+

2017-04-26 22:27:41 547

原创 numpy学习笔记

slicinga[:-1]#去掉最后一个a[:7:2]#[0,7]以2为增值累加索引In [48]: a[::-1]Out[48]: array([9, 8, 7, 6, 5, 4, 3, 2, 1, 0])In [49]: a[::-2]Out[49]: array([9, 7, 5, 3, 1])In [53]: a[:4:-2] #到4截至Out[53]: array([9, 7

2017-04-25 09:38:16 382

原创 c++ string常用函数总结

构造函数 string(const char *s); //用c字符串s初始化string(int n,char c); //用n个字符c初始化substr string substr(int pos = 0,int n = npos) const;//返回pos开始的n个字符组成的字符串只有当pos = npos即最后一个空时才返回空串.枚举的时候一定要小心

2017-04-24 22:18:06 474

原创 HDu 6000 Wash(贪心)

题目链接Wash分析每台洗衣机的工作时间是独立的可以先预处理出每件衣服洗完的时间.然后用贪心的想法把最后洗完的衣服安排给最早烘干的.#include <cstdio>#include <iostream>#include <vector>#include <queue>#include <algorithm>#include<cmath>#include <cstring>#inclu

2017-04-23 23:47:13 657

原创 codeforces #410 div2 题解(数学)

官方题解http://codeforces.com/blog/entry/51652 这一次的题目基本都有很强的思维性。编写代码很简单.C给你一个序列a1,a2,…,ana_1,a_2,\dots ,a_n,有一个操作将(ai,ai+1)变为(ai−ai+1,ai+ai+1)(a_i,a_{i+1}) 变为(a_i-a_{i+1},a_i+a_{i+1}),问最少的操作次数使得 gcd(a1,a

2017-04-22 10:52:08 451

原创 UVA 11728 - Alternate Task(线性筛)

题目链接11728 - Alternate Task分析其实就是一个水题……由于因子和函数是积性函数，我们可以筛出所有1−10001-1000以内的因子和函数的值，然后二分查找一下就可以了.#include <cstdio>#include <iostream>#include <vector>#include <queue>#include <algorithm>#include<cma

2017-04-21 21:01:00 437

原创离散对数和原根

原根与阶阶定义设 (a,m)=1(a,m)=1, 满足 ax≡1(mod m)a^x \equiv 1(mod\ m) 的最小的 xx，称为a对m的阶，记为 ordm(a)ord_{m}(a) 当 ordm(a)=ϕ(m)ord_{m}(a)=\phi(m) 时称为a为m的原根.简单性质ax≡1⇔ordm(a)∣xa^x\equiv 1\Leftrightarrow ord_m(a) \mid

2017-04-19 23:41:59 4858 4

原创 zoj 3940E - Modulo Query (数学)

题目链接E - Modulo Query 分析不要问我为什么，看了题解再yy一下你就懂了. 题解AC code#include <cstdio>#include <iostream>#include <vector>#include <queue>#include <algorithm>#include<cmath>#include <cstring>#include <map>#

2017-04-19 20:34:47 822

原创 codeforces 801 E. Vulnerable Kerbals(DAG最长路+解线性同于方程)

题目链接codeforces 801 E. Vulnerable Kerbals分析假设我们求出的前缀乘积的顺序是这样的p1,p2,…,pkp_1,p_2,\dots ,p_k.产生的序列为a1,a2,…,aka_1,a_2,\dots ,a_k,那麽就一定有,pi−1∗ai≡pi(MOD m)p_{i-1}*a_i\equiv p_i(MOD\ m)则一定会有 gcd(pi−1,m)∣pigcd(

2017-04-17 13:52:41 952

原创莫比乌斯(mobius)笔记

mobius反演的本质是容斥原理，这在《组合数学》里面有提到基本公式定理mobius反演公式 F(n)f(n)=∑d∣nf(d)=∑d∣nμ(n/d)F(d)=∑d∣nμ(d)F(n/d) \begin{align} F(n) &= \sum_{d\mid n}f(d)\\ f(n)&=\sum_{d\mid n}\mu(n/d)F(d)\\ &=\sum_{d\

2017-04-15 14:44:42 1048

原创《The Art of R Programming》R初学者笔记(1-6章)

向量seq与length length(x)返回向量x的长度，seq(x)返回x的长度的一个向量.可能为空。它们的区别在于> 1:length(x)[1] 1 0> seq(x)integer(0)循环补齐当用矩阵加法或者向量加法是如果有一个向量很短，回将其扩展成与另一个向量一样长.其中补齐是按照重复当前向量进行.向量索引> x <- c(1,2,3,4,5)> x[c(1,2)

2017-04-12 00:20:50 1060

原创 UVA 10791 - Minimum Sum LCM(因子分解)

题目链接10791 - Minimum Sum LCM题目分析水题，刚开始看错题目了，以为和经典题目GCD&LCM一样，原来只需要求至少这么多个因数和相加最小.显然就是把所有的priip_i^{r_i}加起来就行了.,注意，因子数目只有一个的时候，还要加上1.AC code#include <cstdio>#include <iostream>#include <vector>#include

2017-04-06 23:14:58 497

原创 Uva 12716 - GCD XOR(数域筛法+qijiyinqiao)

题目链接GCD XOR题目分析这个题是《刘汝佳算法经典》上面的题目，一个数域筛法的经典题目.我做一个搬运工~ 一个基本的想法是枚举每个公约数c,及其倍数ac,及其倍数a,由异或的特性我们有，a XOR b=c,b=a XOR c,a\ XOR\ b=c,b = a\ XOR\ c,，然后在判断一下是否满足gcd(a,b)==cgcd(a,b)==c,把答案加上.可以达到O(n(lgn)2)O(n

2017-04-06 19:44:46 1579

原创 poj 1150 The Last Non-zero Digit(阶乘取模运算的简单性质及递归技巧)

题目链接The Last Non-zero Digit题目分析参见http://www.cppblog.com/abilitytao/archive/2009/10/31/99907.html 题解写的非常好.我就简要总结一下自己的收获.对于n!分解成一个数x的幂次我们可以递归的去计算.因为只有x的倍数相乘才会产生x的幂次.及x,x2,x3,…,x,x^2,x^3,\dots , 因此 f(n,

2017-04-04 21:10:43 1095

原创同余方程总结

题目链接Recurrent Function题目分析可以参考具体数学的公式设mm的dd进制数为mm[0,1,2,…,n−1]k的d进制数为kk[0,1,2,…,m−1]mm[0,1,2,\dots,n-1]k 的d进制数为kk[0,1,2,\dots,m-1],由于置换数组a为{1,2,…,d−1},没有0所以如果说k和m的d进制长度不相等那么一定是不能变换过去的a为\{1,2,\dots,

2017-04-02 18:04:37 7131

taotao 的大学墓志