摘要： d模f的乘法逆元int Extend_Eulid(int d,int f){ int x1,x2,x3,y1,y2,y3 ; x1=1,x2=0,x3=f,y1=0,y2=1,y3=d ; while(y3 && y3!=1) { int q=x3/y3 ; int t1,t2,t3 ; t1=x1-q*y1,t2=x2-q*y2,t3=x3-q*y3 ; x1=y1,x2=y2,x3=y3 ; y1=t1,y2=t2,y3=t3 ; } if(!y3)return -1 ; ret... 阅读全文

摘要： http://acm.hdu.edu.cn/showproblem.php?pid=1452原来真心没见过这种题，不会做，非常帅gcd(a,b)==1 && s(a,b)==s(a)*s(b)满足这种条件的s叫做积性函数，本题求的因子和就是一个积性函数接着有一个结论if(prime[p])s(p^n)=1+p^1+p^2+p^n=(p^(n+1)-1)/(p-1)s(2004^n)=s(2^(2n))*s(3^n)*s(167^n)其中，167和22关于29同余所以，s(2004^n)=s(2^(2n))*s(3^n)*s(2^n)a=s(2^(2n))=(2^(2n+1)-1 阅读全文