#include <bits/stdc++.h>
#define LL long long
using namespace std;
const LL MOD = 19940417;
const LL inv6 = 3323403;
const LL inv2 = 9970209;
LL n, m;

LL Calc1(LL x, LL y) {
	LL ret = 0;
	for (LL l = 1, r; l <= y; l = r + 1) {
		r = min(y, x / (x / l));
		ret = (ret + (x / l) * (l + r) % MOD * (r - l + 1) % MOD * inv2 % MOD) % MOD;
	}
	return ret;
}

LL f(LL x) { return x * (x + 1) % MOD * (x * 2 + 1) % MOD * inv6 % MOD; }

LL Calc2(LL x, LL y) {
	LL ret = 0;
	for (LL l = 1, r; l <= x && l <= y; l = r + 1) {
		r = min(y / (y / l), x / (x / l));
		ret = (ret + (x / l) * (y / l) % MOD * (f(r) - f(l - 1)) % MOD) % MOD;
	}
	return ret;
}

LL Calc(LL x) {
	return x * x % MOD - Calc1(x, x);
}

int main() {
	ios :: sync_with_stdio(0); cin.tie(0); cout.tie(0);
	cin >> n >> m; 
	if (n > m) swap(n, m);
	LL Ans = Calc(n) * Calc(m) % MOD - n * n % MOD * m % MOD + m * Calc1(n, n)
		+ n * Calc1(m, n) - Calc2(n, m);
	Ans = (Ans + MOD) % MOD;
	cout << Ans;
	return 0;
}