Mail.Ru Cup 2018 Round 3 B. Divide Candies

题目链接
分析一下题意可以得到题目要求的是满足下面这个公式的不同的 $i ， j$ 的方案数;
即 $(i^{2} + j^{2}) m o d <mtext> </mtext> m = 0 <mtext> </mtext> (n <mtext> </mtext> \leq <mtext> </mtext> i, j \leq <mtext> </mtext> n)$ .
当然 $n^{2}$ 是可以解决的，但是数据范围不允许我们这么做，于是我们可以换个思路，如果一个数 $i$ 满足 $<mtext> </mtext> i m o d m = 0$ ，那么 $(i + k * m)$ 也必然是等于 $0$ 的，因为任意一个数 $x$ 的结果必然是小于 $m$ 的，于是我们可以 $m^{2}$ 的复杂度来解决这个问题。在 $m$ 的范围内找出所有满足题意的 $i, j$ 。然后对每一对 $i, j$ 进行累加贡献。然后这个问题就解决了。

#include<bits/stdc++.h>

#define LL long long
#define fi first
#define se second
#define mp make_pair
#define pb push_back

using namespace std;

LL gcd(LL a,LL b){return b?gcd(b,a%b):a;}
LL lcm(LL a,LL b){return a/gcd(a,b)*b;}
LL powmod(LL a,LL b,LL MOD){LL ans=1;while(b){if(b%2)ans=ans*a%MOD;a=a*a%MOD;b/=2;}return ans;}
const int N = 2e5 +11;
LL n,m,ans;
int main(){
	ios::sync_with_stdio(false);
	cin>>n>>m;
	for(int i=0;i<m;i++){
		for(int j=0;j<m;j++){
			if((i*i+j*j)%m==0&&i<=n&&j<=n){
				LL a=(n-i)/m+1,b=(n-j)/m+1;
				if(!i)a--;
				if(!j)b--;
				ans+=a*b;
			}
		}
	}
	cout<<ans<<endl;
	return 0;
}