摘要：
题意：给出a、p、d、m 求a^x=d(mod p) 0m则输出0如果x 2 #include 3 #include 4 typedef long long ll; 5 typedef unsigned long long ull; 6 const ll mo=200003,N=20000; 7 ll a,p,d,m,sum=1,sq,hash[mo],hat[mo],pri[N],save[N],bo[N+1],tot,add; 8 ll gcd(ll a,ll b){ return b ? gcd(b,a%b) : a; } 9 ll extgcd(ll &x,ll &.. 阅读全文
posted @ 2013-10-18 23:55
g_word
阅读(282)
评论(0)
推荐(0)
摘要：
题意：求前n项的欧拉函数之和题解：预处理出所有欧拉函数 赤裸裸的模版题- - 没什么好说的代码： 1 #include 2 typedef long long ll; 3 const ll N=1000001; 4 ll n,phi[N],pri[N],bo[N]; 5 void makepri(){ 6 for (ll i=2;i<=N;i++){ 7 if (!bo[i]){ 8 pri[++pri[0]]=i; 9 phi[i]=i-1;10 }11 for (ll j=1;j<=p... 阅读全文
posted @ 2013-10-18 23:37
g_word
阅读(228)
评论(0)
推荐(0)
摘要：
题意：给出n个模方程x=a(mod r) 求x的最小解题解：这就是个线性模方程组的模版题- - 但是有一些要注意的地方extgcd算出来的解x可能负数 要让x=(x%mo+mo)%mo而且mo不是等于lcm(r1,r2) 而是r2/gcd(r1,r2)代码： 1 #include 2 typedef long long ll; 3 ll n,a,r; 4 ll extgcd(ll &x,ll &y,ll a,ll b){ 5 if (!b){ 6 x=1,y=0; 7 return a; 8 }else{ 9 ll res=ex... 阅读全文
posted @ 2013-10-18 23:35
g_word
阅读(180)
评论(0)
推荐(0)
摘要：
经过两天的努力 终于把AC大神的课件都看完了 感动啊啊啊TAT 顿时感觉智力上升了一个层次 之后看到数论的题目终于可以不用放弃治疗了~~~ 我会说我刚刚用了十分钟就把欧拉函数看完了吗~~~（其实之前就会-。-）欧拉函数也就是phi(n) 表示小于等于n的且与n互质的数的个数欧拉函数的公式是 phi(n)=n*(1-1/p1)*(1-1/p2)*...*(1-1/pn)要证明这个公式 我们需要先证明一些简单的性质1.若p为素数 phi(p)=p-1小于某个素数的数都和这个素数互质- - 显然2.phi(p^k)=(p-1)*p^(k-1)小于等于p^k的数中不与p互质的数肯定都是p的倍数这些数有 阅读全文
posted @ 2013-10-18 00:10
g_word
阅读(473)
评论(6)
推荐(1)
浙公网安备 33010602011771号