#include<bits/stdc++.h>
#define Gc() getchar()
#define Me(x,y) memset(x,y,sizeof(x))
#define Mc(x,y) memcpy(x,y,sizeof(x))
#define d(x,y) ((m)*(x-1)+(y))
#define R(n) (rnd()%(n)+1)
#define Pc(x) putchar(x)
#define LB lower_bound
#define UB upper_bound
#define fi first
#define se second
#define eb emplace_back
#define all(x) x.begin(),x.end()
using namespace std;using ll=long long;using db=double;using lb=long db;using ui=unsigned;using ull=unsigned long long;using pii=pair<int,int>;
const int N=1e6+5,M=N*120+5,K=1000+5,mod=998244353,Mod=mod-1;const db eps=1e-9;const int INF=1e9+7;mt19937 rnd(263082);
#define Tp template<typename T>
#define Ts template<typename T,typename... Ar>
namespace Debug{
	Tp void _debug(char* f,T t){cerr<<f<<'='<<t<<endl;}
	Ts void _debug(char* f,T x,Ar... y){while(*f!=',') cerr<<*f++;cerr<<'='<<x<<",";_debug(f+1,y...);}
	#ifdef LOCAL
	#define gdb(...) _debug((char*)#__VA_ARGS__,__VA_ARGS__)
	#else 
	#define gdb(...) void()
	#endif
}using namespace Debug;
ll mpow(ll x,int y=mod-2){ll ans=1;while(y) y&1&&(ans=ans*x%mod),y>>=1,x=x*x%mod;return ans;}
ll n,k;
int getid(ll x){
	return x<=k?x:k+n/x;
}
namespace Mu{
	ll f[N];
	ll id[N];int ih;
	void init(){
		for(ll l=1,r;l<=n;l=r+1){
			r=n/(n/l);
			id[++ih]=n/l;
		}
		reverse(id+1,id+ih+1);
		for(int i=1;i<=ih;i++){
			int v=getid(id[i]);
			f[v]=1;
			for(ll l=2,r;l<=id[i];l=r+1){
				r=id[i]/(id[i]/l);
				f[v]-=f[getid(id[i]/l)]*(r-l+1);
			}
		}
	}
}
int vis[N];ll f[N],g[N],dp[N];
int flag[N],pr[N],ph,mu[N];
void init(){
	mu[1]=1;
	for(int i=2;i<=n;i++){
		if(!flag[i]) pr[++ph]=i,mu[i]=-1;
		for(int j=1;j<=ph&&i*pr[j]<=n;j++){
			flag[i*pr[j]]=1;
			if(i%pr[j]==0) break;
			mu[i*pr[j]]=-mu[i];
		}
	}
}
void calc(ll n){
	int v=getid(n);
	if(vis[v]) return;
	vis[v]=1;
	for(ll l=1,r;l<=n;l=r+1){
		r=n/(n/l);
		g[v]+=(Mu::f[getid(r)]-Mu::f[getid(l-1)])*(mpow(2,n/l)-1)%mod;
	}
	g[v]=(g[v]%mod+mod)%mod;
	for(int i=2;i<=n;i++){
		calc(n/i);int u=getid(n/i);
		f[v]+=(mpow(2,i/2-1)-1)*(f[u]+f[u]-dp[u])*2%mod;
		f[v]+=(mpow(2,i/2-1)-1)*(g[u]*2-1)%mod;
		f[v]+=f[u]-dp[u];
		f[v]+=g[u]-mpow(2,n/i-1);
		dp[v]+=f[u]*3-dp[u];
		dp[v]+=mpow(2,n/i-1)-1;
		// if(n==5) gdb(f[n]);
	}
	f[v]=(f[v]%mod+mod)%mod;
	dp[v]%=mod;
	// gdb(n,f[n],g[n],dp[n]);
}
void Solve(){
	scanf("%lld",&n);
	k=sqrt(n);
	Mu::init();
	init();
	calc(n);
	printf("%lld\n",(f[getid(n)]+mpow(2,n-1))%mod);
}
int main(){
	int t=1;
	// scanf("%d",&t);
	while(t--) Solve();
	cerr<<clock()*1.0/CLOCKS_PER_SEC<<'\n';
}