#ifdef MAIN bool __multi__ = 0; namespace XK { const int MD = 998244353; const ll mod = 998244353; struct mm { ll x; mm(ll x_ = 0) : x(x_ % mod) { if(x < 0) x += mod; } friend mm operator+(mm a, mm b) { return a.x + b.x; } friend mm operator-(mm a, mm b) { return a.x - b.x; } friend mm operator*(mm a, mm b) { return a.x * b.x; } friend mm operator/(mm a, mm b) { return a * b.inv(); } friend mm& operator+=(mm& a, mm b) { return a = a.x + b.x; } friend mm& operator-=(mm& a, mm b) { return a = a.x - b.x; } friend mm& operator*=(mm& a, mm b) { return a = a.x * b.x; } friend mm& operator/=(mm& a, mm b) { return a = a * b.inv(); } mm inv() const { return pow(mod - 2); } mm pow(ll b) const { mm a = *this, c = 1; while(b) { if(b & 1) c *= a; a *= a; b >>= 1; } return c; } }; void solve() { int N, M; in(N, M); mm ans = 0; rep(i, N) rep(j, M) { ans += (j + 1) * (i + 1) / __gcd(i + 1, j + 1); } out(ans.x); } }; signed main() {ios::sync_with_stdio(0);cin.tie(0);fixed(cout).precision(12);int t = 1;if(__multi__) cin >> t;while(t --) XK::solve();} #else #include "cassert" #include "cmath" #include "cstdint" #include "cstdio" #include "cstdlib" #include "cstring" #include "algorithm" #include "bitset" #include "chrono" #include "complex" #include "deque" #include "functional" #include "iostream" #include "limits" #include "map" #include "numeric" #include "queue" #include "random" #include "set" #include "sstream" #include "string" #include "unordered_map" #include "unordered_set" #include "utility" #include "vector" #include "array" using namespace std; #define int long long using ll = long long; using ull = unsigned long long; const ll INF = 1ll << 60; const ll LINF = 0x1fffffffffffffff; const ll MINF = 0x7fffffffffff; template bool chmax(A& l, const B& r){ return r > l ? l = r, 1 : 0; } template bool chmin(A& l, const B& r){ return r < l ? l = r, 1 : 0; } #define sz(x) ssize(x) #define rep(i, a) for(ll i = 0; i < (a); i ++) #define Rep(i, a, b) for(ll i = (a); i < (b); i ++) #define rrep(i, a, b) for(ll i = (b); i --> (a); ) #define all(x) begin(x), end(x) #define fst first #define snd second #define pb push_back template using V = vector; template using AR = array; namespace IO { template void in(Ts&... t); [[maybe_unused]] void print(){} template void print(const T& t, const Ts&... ts); template void out(const Ts&... ts){ print(ts...); cout << '\n'; } namespace IO{ #define VOID(a) decltype(void(a)) struct S{ S(){ cin.tie(nullptr)->sync_with_stdio(0); fixed(cout).precision(12); } }S; template struct P : P{}; template<> struct P<0>{}; template void i(T& t){ i(t, P<3>{}); } void i(vector::reference t, P<3>){ int a; i(a); t = a; } template auto i(T& t, P<2>) -> VOID(cin >> t){ cin >> t; } template auto i(T& t, P<1>) -> VOID(begin(t)){ for(auto&& x : t) i(x); } template void ituple(T& t, index_sequence){ in(get(t)...); } template auto i(T& t, P<0>) -> VOID(tuple_size{}){ ituple(t, make_index_sequence::value>{}); } template void o(const T& t){ o(t, P<4>{}); } template void o(const char (&t)[N], P<4>){ cout << t; } template void o(const T (&t)[N], P<3>){ o(t[0]); for(size_t i = 1; i < N; i++){ o(' '); o(t[i]); } } template auto o(const T& t, P<2>) -> VOID(cout << t){ cout << t; } template auto o(const T& t, P<1>) -> VOID(begin(t)){ bool first = 1; for(auto&& x : t) { if(first) first = 0; else o(' '); o(x); } } template void otuple(const T& t, index_sequence){ print(get(t)...); } template auto o(T& t, P<0>) -> VOID(tuple_size{}){ otuple(t, make_index_sequence::value>{}); } #undef VOID } template void in(Ts&... t){ (IO::i(t), ...); } template void print(const T& t, const Ts&... ts){ IO::o(t); (IO::o((cout << ' ', ts)), ...); } #undef unpack constexpr ll debug_const(ll judge, ll debug) { #ifdef DEBUG return debug; #else return judge; #endif } #ifdef DEBUG ll __lg(ull x){ return 63 - __builtin_clzll(x); } #define debug(...) { print(#__VA_ARGS__); print(":"); out(__VA_ARGS__); } #else #define debug(...) void(0) #endif #define dbg debug } using namespace IO; #define MAIN #include __FILE__ #endif