shiyuan的专栏

题目：hdu 4767 Bell思路：两个公式 B(n+p)%p=(B(n)+B(n+1))%p  B(p^m+n)%p=(m*B(n)+B(n+1))%p首先明确的是，对于取模数，质因数分解 ，直接套用模板了。#include #include #include #include #include #include #include using namespace

题目：CF K. Killer Challenge思路：状态dp前8项素因子的积就超过了100w，所以状态至多只有2^8个#include #include #include #include #include #include using namespace std;#define maxn 1000010bool vis[maxn];vecto

题目：hdu 2814 Interesting Fibonacci题意：求F（a^b）^（F(a^b) ^ (n-1)）%c  （F（n）表示斐波那契数列，1思路：指数循环节和斐波那契数列循环节F(n)%c 的情况，根据斐波那契数列的循环节求出来，做法跟12年四川现场赛那题一样。式子中存在两个F（a^b） ，对于第一个，求对于mod c的循环节，第二个根据指数循环节，求mod ph

题目：hdu 2447 K-dimension number题意：一个数有k个不同的约数定义为一个k-dimension数，输入n,k，求出第n大的k-dimension数，题目说的有点含糊，保证输入的k满足k是一个3-dimension数，也就是说k只能是1或者p，或者p^2，且最大为97*97.思路：对于答案，质因数分解来考虑，ans=x1^p1 * x2^p2 * ... * xn^

题目：spoj 3871. GCD Extreme思路：欧拉函数 做法 同 poj2480时限好像是3s 我的AC代码运行的好像是接近5s  记得用%lld 用%I64d会Wa#include #include #include #include #include using namespace std;const int maxn=1000001;long

题目：Light oj 1132 - Summing up Powers题意：求这个式子        思路：构造矩阵。。。。。#include #include #include #include #include using namespace std;const long long mod = (

题目：virtual contest 里面的 Repackaging      或者 uva 上的10089 - Repackaging题意：输入n组包的大小，有三个参数a,b,c，现在每个包拿若干个，然后拆分，组成新的包要使新的a,b,c全部相等。相当于把问题转换成了下面这个：           思路：       好吧，这题我不会，第一反应

题目：《挑战程序设计》P143有 n 个物品的重量和价值分别是 wi 和 vi 。从中选出 k 个物品使得单位重量的价值最大11思路：二分设最后选定物品的集合为S用 sigma(i in S) V[i]   所有的价值和用 sigma(i in S) W[i]  所有的质量和ans = sigma(i in S) V[i] /  sigma(i in S

题目： Ipad,IPhone思路：这题。。。好坑反正，WA了刚好20次。。。题目给的式子，分成两部分，对于第一部分，直接快速幂即可，对于后面的部分，后面括号里是整数的证明 出题人 AekdyCoin 用二次剩余进行了解说，请点击 here 对于后面的部分，是不是整数我给直接忽略了。。。不然怎么进行取余呢。。事实证明，前半部分是有存在感的，根据欧拉准则，即二次剩余的

题目：1951: [Sdoi2010]古代猪文tag：Lucas定理，快速幂，扩展欧几里得，中国剩余定理第一次做这么恶心的数论题，真心写的想吐了。题目还是粘一下吧，还是蛮好的一个题呀：1951: [Sdoi2010]古代猪文Time Limit: 1 Sec  Memory Limit: 64 MBSubmit: 681  Solved: 236[Sub

题目：hdu 3369 Robot题意：告诉你从周几开始，到第n天总的学会单词的个数，其中第i天学会i^k个单词，周末学会单词为0思路：给的k的范围比较小，n的范围比较大开始想的是，枚举k，但是这样对于n并没有影响所以，我们要处理的应该是，i^k+(i+7)^k+(i+14)^k+....(i+7*j)^k这样的一个通式由于k的范围比较小，展开成k+1项二项式表示之后可以构造矩