时光真疯狂, 我一路执迷于匆忙.

忘了的时候来抄模板（我记性差~）scanf("%d%d",&n,&m);if (n>m) swap(n,m);mu[1]=1;fo(i,2,n) { if (!f[i]) p[++tot]=i,mu[i]=-1; fo(j,1,tot) { int k=p[j]*i;if (k>n) break;f[k]=1; if (!(i%p[j])) break;m

Description给出t组询问，每组询问给出n个数，a1~an，和模数p，求a1^a2^….an mod p的值。 tSolution这样我们只需要快速计算axmodpa^x\mod p的值就可以了。 如果gcd(a,p)=1\gcd(a,p)=1的话，那么aφ(p)≡1modpa^{\varphi(p)}≡1\mod p 所以ax≡axmodφ(p)modpa^x≡a^{x

Description求∑i=1n∑d|igcd(d,id)\sum_{i=1}^{n}\sum_{d|i}gcd(d,{i\over d}) n<=10^11Solution首先，我们发现gcd中的两个东西是所有乘积不超过n的数对，即 Ans=∑i∑j,i∗j<=ngcd(i,j)Ans=\sum_{i}\sum_{j,i*j<=n}gcd(i,j) 然后Ans=∑i=1n∑j=1⌊ni⌋g

Description在平面直角坐标系中，点(x,y)的代价定义为它和原点的连线中经过多少的其他整点个数*2+1。求横坐标在1~n且纵坐标在1~m的所有点的代价和。 n,m<=10^5Solution首先，一个点的代价就是gcd(x,y)*2-1（蒟蒻不会证，各路大犇们带带我呗） 然后，我们就变成了要求∑i=1n∑j=1mgcd(i,j)\sum_{i=1}^{n}\sum_{j=1}^{m}g

Description求∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i,j) n,m<=10^7Solution（注意，以下内容默认n<=mn<=m） 看着lcm不爽，把它变一变： ∑i=1n∑j=1mijgcd(i,j)\sum_{i=1}^{n}\sum_{j=1}^{m}{ij\over \gcd(i,j)} 莫比乌斯反演常用，枚

Description给出n和m，求n!的范围中与m!互质的数的个数。答案对r取模。 n,m<=1e7，r>n>=m，且r为质数。Solution首先，如果gcd(a,b)=1，那么gcd(a+b,b)=1. 所以，1~km与m互质的数的个数为k∗φ(m)k*\varphi(m) 答案就是n!m!φ(m!){n!\over m!}\varphi(m!) 然后，φ(m!)=m!∏p<=n且p为

Description求∑i=1n(−1)∑mj=1d(i∗j)\sum_{i=1}^{n}(-1)^{\sum_{j=1}^{m}d(i*j)} 其中d(i)表示i的因子个数。 n<=10^7,m<=10^14Solution既然是-1的次幂，那么我们就来分析一下奇偶性吧。。。 这里有一个很（不）显然的性质，d(n)是奇数当且仅当n是一个完全平方数。 然而我比赛3个小时都没有想出来233

Description求∑i=1n∑j=1me(gcd(i,j)是质数)\sum_{i=1}^{n}\sum_{j=1}^{m}e(gcd(i,j)是质数) T<=1000,n,m<=5*10^6Solution经典反演套路题，貌似没有其他做法（其他做法的大爷不要鄙视蒟蒻=w=） 复习一下莫比乌斯反演，所以就来打了这道题 首先约定n<mn<m 设FdF_d表示∑ni=1∑mj=1e(gcd(