#include #include using namespace std; using ll = long long; using mint = atcoder::modint; int main() { ll x,y,m; cin >> x >> y >> m; mint::set_mod(m); mint ans = 0; ll mx = max(x, y) / 2 + 1; for (ll m=0; m<=mx; m++) { for (ll n=m+1; n<=mx; n++) { if (m % 2 == n % 2 || gcd(n, m) != 1) continue; ll ai = n*n - m*m; ll bi = 2*n*m; if (ai > mx || bi > mx) break; ll mul = 1; while (mul*max(ai, bi) <= mx) { ll a = ai*mul; ll b = bi*mul; if (a*2 < x && b*2 < y) { if (b == 0) ans += (x - 2 * a) * (y - 2 * b) * 2; else ans += (x - 2 * a) * (y - 2 * b) * 4; } if (a*2 < y && b*2 < x) { if (b == 0) ans += (y - 2 * a) * (x - 2 * b) * 2; else ans += (y - 2 * a) * (x - 2 * b) * 4; } mul++; } } } cout << ans.val() << "\n"; }