mtxxxx的blog

#include #include #include using namespace std; const int N=1000005; long long eu[N],ans[N]; void eular() { int i,j; for(i=2;i<N;i++) if(!eu[i]) { for(j=i;j<N;j

//题目网址 http://acm.hdu.edu.cn/showproblem.php?pid=5673//以下参考过网上代码 http://m.blog.youkuaiyun.com/article/details?id=51227707#include #include #include #include using namespace std; const int mod=1000000007

//在求形如C(n,k)%p这样的式子时，如果p值是一个可变的素数， //而且不是特别大，就可以使用Lucas定理 #include #include #include #include using namespace std; long long pow_mod(long long a,long long i,long long p) { if(i==0)return 1%p

数学题真的好难啊...证明参考博客   http://blog.youkuaiyun.com/guhaiteng/article/details/52400679 #include #include #include using namespace std; typedef long long ll; const int MOD=1e9+7; ll pow_mod(ll a,ll b,ll mod)