P6860 象棋与马

$Result$

---

$Hyperlink$

【LGR-078】洛谷 10 月月赛 II T4
题目链接

---

$Description$

一张无限大的网格，若从$(x,y)$出发，每一步可以类似于马一样移动，长度分别为$a,b$（对于经常提起的中国象棋中的马的话，则有$a=1,b=2$或$a=2,b=1$），若能到达这个网格中的任意一节节点，则$p(a,b)==1$

求${\sum }_{i=1}^n{\sum }_{j=1}^{n}p(i,j)$

数据范围：
$20\%,n\le 3000$
$50\%,n\le 10^7$
$100\%,n\le 10^{11}$

---

$Solution$

$20pts$

能从$(x,y)$移动到任意一个格子的充要条件是它能位移到$(x+1,y)$或$(x,y+1)$
也就是要找出一种位移方案，使得这么多步后，位移到一个格子$(x+ma+nb,y+ea+fb)$，使得$ma+nb=1$且$ea+fb=0$，对于前者，我们很容易得到$gcd(a,b)=1$

对于后者，略加思索，可以发现：$a,b$不能同时为奇数

如此这般，我们预处理3000以内的这样的$a,b$，为了节省效率，我们只算$a>b$的，最后再乘二即可

时间复杂度：$O(n^{2}logn)$ 实际上常数较小，可以拿到$20pts$

#include<cctype>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
#define ULL unsigned long long
using namespace std;int Gcd[3010][3010],k,T;
ULL sum,s[3010],n;
inline int gcd(int x,int y)
{
	if(Gcd[x][y]) return Gcd[x][y];//似乎没有多大用处的记忆化
	return Gcd[x][y]=y?gcd(y,x%y):x;
}
inline LL read()
{
	char c;LL d=1,f=0;
	while(c=getchar(),!isdigit(c)) if(c=='-') d=-1;f=(f<<3)+(f<<1)+c-48;
	while(c=getchar(),isdigit(c)) f=(f<<3)+(f<<1)+c-48;
	return d*f;
}
signed main()
{
	for(register int i=1;i<=3000;i++)
	{
		for(register int j=1;j<=i;j++)
		{
			if((i&1)&&(j&1)) continue;
			sum+=gcd(i,j)==1;
		}
		s[i]=sum;
	}
	T=read();
	while(T--)
	{
		n=read();
		printf("%llu\n",s[n]*2);
	}
}

$50pts$

假设我们已经确定了$a$，考虑怎样的$b$是合法的
刚刚我们已经得出了$a,b$互质，但又有$a,b$不同时为奇数这个限制

那么假设$a$是一个偶数，则它的贡献是不大于它且与它互质的自然数，即$\phi (a)$
假设$a$是一个奇数，则它的贡献是不大于它且与它互质的偶数，这貌似很难办，考场的时候我是推规律推出它为$\frac{\phi (2a)}{2}$【有能力证明的dalao可以在下方提一下】

如此，我们只需要线性筛出$2\times 10^7$内的$\phi$，求个前缀和即可

时间复杂度：$O(n)$，可以拿到$50pts$

#include<cctype>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
#define ULL unsigned long long
using namespace std;int T,phi[20000010],prime[2000010],vis[20000010],m;
ULL n,s[10000010];
inline LL read()
{
	char c;LL d=1,f=0;
	while(c=getchar(),!isdigit(c)) if(c=='-') d=-1;f=(f<<3)+(f<<1)+c-48;
	while(c=getchar(),isdigit(c)) f=(f<<3)+(f<<1)+c-48;
	return d*f;
}
inline void prework()
{
	phi[1]=1;
	for(register int i=2;i<=20000000;i++)
	{
		if(vis[i]==0) {vis[i]=i;prime[++m]=i;phi[i]=i-1;}
		for(register int j=1;j<=m&&prime[j]*i<=20000000;j++)
		{
			if(prime[j]>vis[i]) break;
			vis[i*prime[j]]=prime[j];
			phi[i*prime[j]]=phi[i]*(i%prime[j]?prime[j]-1:prime[j]);
		}
	}
	return;
}
signed main()
{
	prework();
	for(register int i=1;i<=1e7;i++)
	{
		if(i&1) s[i]=phi[i*2]/2;
		else s[i]=phi[i];
		s[i]+=s[i-1];
	}
	T=read();
	while(T--)
	{
		n=read();
		printf("%llu\n",s[n]*2);
	}
}

$100pts$

基于$50pts$的思路，多利用些欧拉函数的性质即可

若$i$是奇数，$\phi (2i)=\phi (i)$
若$i$是偶数，$\phi (2i)=2\phi (i)$

当$i$为偶数时，$w_i=\phi (i)$
当$i$为奇数时，$w_i=\frac{\phi (2i)}{2}=\frac{\phi (i)}{2}$

设$S(n)$表示答案

则$S(n)=2{\sum }_{i=1}^n{w}_i$，那么$S(n)$就等于所有偶数的$\phi$和的两倍+所有奇数的$\phi$的和

我们可以直接用杜教筛算出奇数+偶数的$\phi$和（即$\phi$的前缀和），然后考虑补上没算的偶数的贡献，这一段的贡献即为$S(\lfloor \frac{n}{2}\rfloor )$，递归处理即可，层数显然是$log$层的

时间复杂度：$O(n^{\frac{2}{3}}logn)$，可以通过本题

---

$Code$

#include<map>
#include<cctype>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ULL unsigned long long
using namespace std;
ULL T,prime[2000010],vis[20000010],m,n,phi[20000010],s[20000010];
map<ULL,ULL>Sphi; 
inline ULL read()
{
	char c;ULL d=1,f=0;
	while(c=getchar(),!isdigit(c)) if(c=='-') d=-1;f=(f<<3)+(f<<1)+c-48;
	while(c=getchar(),isdigit(c)) f=(f<<3)+(f<<1)+c-48;
	return d*f;
}
inline void prework()
{
	phi[1]=1;
	for(register int i=2;i<=20000000;i++)
	{
		if(vis[i]==0) {vis[i]=i;prime[++m]=i;phi[i]=i-1;}
		for(register int j=1;j<=m&&prime[j]*i<=20000000;j++)
		{
			if(prime[j]>vis[i]) break;
			vis[i*prime[j]]=prime[j];
			phi[i*prime[j]]=phi[i]*(i%prime[j]?prime[j]-1:prime[j]);
		}
	}
	for(register int i=1;i<=20000000;i++) 
	{
		if(i&1) s[i]=phi[i*2]/2;
		else s[i]=phi[i];
		s[i]+=s[i-1];
		phi[i]+=phi[i-1];
	}
	return;
}
inline ULL Getphi(ULL n)
{
	if(n<=20000000) return phi[n];
	if(Sphi[n]) return Sphi[n];
	ULL res=0;
	if(n%2==0) res=n/2*(n+1);
	else res=(n+1)/2*n;//这一段要这样子弄一弄，不然可能会爆
	for(ULL l=2,r;l<=n;l=r+1)
	{
		r=n/(n/l);
		res-=(r-l+1)*Getphi(n/l);
	}
	return Sphi[n]=res;
}
inline ULL S(ULL n)
{
	if(n<=1) return 0;
	return Getphi(n)+S(n/2);
}
signed main()
{
	prework();
	T=read();
	while(T--)
	{
		n=read();
		printf("%llu\n",S(n));
	}
}