tju_virus的专栏

题目描述：求出第i个非平方数M，以及小于等于M的所有数的开方之和                                                                   纯数学题，把公式推出来就OK了#include #include #include using namespace std;int main(

题目大意：给一个数N（10^8），求 C(n,0),C(n,1),C(n,2)...C(n,n).中有多少个奇数。解体思路：本题为Lucas定理推导题，我们分析一下 C(n,m)%2,那么由lucas定理，我们可以写成二进制的形式观察，比如 n=1001101，m是从000000到1001101的枚举，我们知道在该定理中C(0,1)=0,因此如果n=1001101的0对应位置的

题目描述：有N个元素的一个集合经过K次置换能变回原来的集合，求k的个数。解题思路：k为置换群中各独立置换群长度的最小公倍数，k的个数即N以内的总和小于10的质数的乘积（包括1）比如N=10时，k可为：1，2，3，2*2，5，2*3，7，2*2*2，3*3，2*5，2*2*3，2*7，3*5，2*2*5，3*7，2*3*5，共16个解法：DP      状态： dp[

题目相当于求n个数的和不超过m的方案数。如果和恰好等于m，那么就等价于方程x1+x2+...+xn = m的解的个数，利用插板法可以得到方案数为：(m+1)*(m+2)...(m+n-1)  = C(m+n-1,n-1) = C(m+n-1,m)现在就需要求不大于m的，相当于对i = 0,1...,m对C(n+i-1,i)求和，根据公式C(n,k) = C(n-1,k)+C(n-1,k

#include #include #include #include typedef long long ll;const int M = 100005; //所取模的上界const int N = 405;int p[N],e[N];int len;int cover[N];int pri[N];int top;int f[M];int finv[M];

给定几个数，问其中素数的个数。由于数量很大，打表肯定过不了，所以必须用Miller_Rabbin随机素数测试的方法，下面是我的C++ Miller_Rabbin模板#include #include #include using namespace std;typedef long long ll;ll exp_mod(ll a,ll b,ll mod){ ll

经典的数学题题意：判断一个数是不是素数，如果是输出Prime，否则输出这个数的质因数中最小的一个题解：用Miller_rabbin算法判断这个数是不是素数，如果不是就用Pollard_rho算法进行因式分解以下为代码，但是在poj上不知道为什么如果我用G++交会报runtime error而用C++交就能过，求大牛指教#include #include #include

题意：给4个数n，r，k，m，求从n中选r个数（至少间隔k ）并将这r个数分成至多m组一共有多少种情况。输入为5 2 3 2时Sample input means you can choose 1 and 4,1 and 5,2 and 5 in the same day.And you can make the machines in the same group or in th

该题是Pollard_rho分解质因数的模板题1、所求长度是N的因子（且大于1小于N）,集合中的元素得两两互质 2、为了尽可能多的选出，每个L的质因子应当只包含N的一个质因子，L是一个质因子的整数次，所以K的值就是N中不同质因子的个数3、要想和最大，那么使得每个L最大，只要使得质因子的指数最大即可如果N本身是质数的话，直接输出1 1用pollard_rho分