#include <iostream>
using namespace std;
typedef long long ll;
const int mx=1000000; const ll mod=998244353;
ll pow(ll x, int e)
{
    if(!e) return 1; ll t=pow(x,e>>1); t=(t*t)%mod;
    if(e&1) t=(t*x)%mod; return t;
}
int u[mx+1]; bool p[mx+1];
int main()
{
    //count no. gcd 1 n-tuples of -k..k
    //t|gcd: (2(k/t)+1)^n; PIE: mobius inversion
    int n; ll k, ans=1; cin>>n>>k;
    for(int i=2; i<=mx; i++) p[i]=1;
    for(int i=2; i<=1000; i++) if(p[i]) for(int j=i+i; j<=mx; j+=i) p[j]=0;
    for(int i=2; i<=mx; i++) if(p[i])
    {
        for(int j=i; j<=mx; j+=i) u[j]++;
        if(i<=1000) for(int j=i*i; j<=mx; j+=i*i) u[j]=-mx;
    }
    for(int i=1; i<=mx; i++) u[i]=(u[i]<0)?0:((u[i]&1)?-1:1);
    for(int i=1; i<=k; i++) if(u[i]) ans+=(pow(2*(k/i)+1,n)-1)*u[i];
    cout<<((ans%mod)+mod)%mod<<'\n';
}