#include #include #define all(a) begin(a),end(a) #define eb emplace_back #define ep emplace #define ef(...) else if(__VA_ARGS__) #define wh(...) while(__VA_ARGS__) #define fin(...) rr pp(__VA_ARGS__) #define RD(T,...) T __VA_ARGS__;li(__VA_ARGS__) #define LL(...) RD(ll,__VA_ARGS__) #define ST(...) RD(str,__VA_ARGS__) #define UL(n,...) u1 __VA_ARGS__;resize({n},__VA_ARGS__);li(__VA_ARGS__) #define UV(n,...) u1 __VA_ARGS__;resize({n},__VA_ARGS__);vi(__VA_ARGS__) #define SS(n,...) strs __VA_ARGS__;resize({n},__VA_ARGS__);li(__VA_ARGS__) #define U2(n,m,...) u2 __VA_ARGS__;resize({n,m},__VA_ARGS__);li(__VA_ARGS__) #define fe(v,e,...) for(au&&__VA_OPT__([)e __VA_OPT__(,__VA_ARGS__)__VA_OPT__(]):v) #define fo(i,...) for(au[i,i##stop,i##step]=rng_fo(0,__VA_ARGS__);i=i##stop;i-=i##step) #define debug(...) pp(#__VA_ARGS__,__VA_ARGS__) #define I(...) template<__VA_ARGS__> #define J typename #define O operator #define au auto #define cs const #define ifc if constexpr #define rr return #define uu using namespace my{uu vo=void;vo main();vo solve();}int main(){my::main();}namespace my{uu namespace std;uu namespace atcoder; uu is=istream;uu os=ostream;uu bo=bool;uu dd=long double;uu ll=long long;uu i128=__int128_t;uu ul=unsigned long long;uu u128=__uint128_t; I(J T)uu v1=vector;I(J T)uu v2=v1>;uu u1=v1;uu u2=v2;uu str=string;uu strs=v1;str sp=" ",nc="",nl="\n"; uu ml=modint998244353;is&O>>(is&i,ml&x){ll t;i>>t;x=t;rr i;}os&O<<(os&o,cs ml&x){rr o<;uu m2=v2; I(J A,J B=A)struct cp{A a={};B b={};cp(){}cp(A a,B b):a(a),b(b){}cp(pairp):a(p.first),b(p.second){} cp&O+=(cs cp&c){a+=c.a,b+=c.b;rr*this;}cp O+(cs cp&c)cs{rr cp{*this}+=c;}bo O==(cs cp&c)cs{rr a==c.a&&b==c.b;} cp&O-=(cs cp&c){a-=c.a,b-=c.b;rr*this;}cp O-(cs cp&c)cs{rr cp{*this}-=c;}cp O-()cs{rr cp(-a,-b);} bo O<(cs cp&c)cs{rr a!=c.a?a(cs cp&c)cs{rr a!=c.a?a>c.a:b>c.b;}friend os&O<<(os&o,cs cp&c){rr o<;uu cls=v1;I(J T)uu qmax=priority_queue;I(J T)uu qmin=priority_queue,greater<>>; I(J T)struct set:std::set{set(v1a={}){fe(a,e)this->ep(e);}};I(J T)struct uset:unordered_set{uset(v1a={}){fe(a,e)this->ep(e);}}; I(J T)struct mset:multiset{mset(v1a={}){fe(a,e)this->ep(e);}vo erase(cs T&x){if(this->contains(x))multiset::erase(this->find(x));}}; I(J T,J U=T)uu map=std::map;I(J T,J U=T)uu umap=unordered_map;I(J T)ll len(cs T&a){rr size(a);} ll oo=3e18;dd ee=1e-12;u1 dx{-1,0,1,0,-1,1,1,-1},dy{0,-1,0,1,-1,-1,1,1};str Yes(bo y=1){rr y?"Yes":"No";}str No(){rr"No";} ll rnd(ll l,ll r=oo){static ul x=76,y=53,z=20,w=41;ul t=x^(x<<11);x=y,y=z,z=w,w=(w^(w>>19))^(t^(t>>8));if(r==oo)r=l,l=0;rr w%(r-l)+l;} ll l2(u128 x){x|=1;ll r=0;wh(x)x>>=1,++r;rr r;}bo odd(ll x){rr x&1;}ll pwm1(ll x){rr 1-2*(x&1);}I(J T)bo at(T x,ll i){rr x>>i&1;} ll l10(u128 x){x|=1;ll r=0;wh(x)x/=10,++r;rr r;}bo eve(ll x){rr~x&1;}ll pw2(ll n){rr 1LL<=0?x/y:(x+1)/y-1;}au rng_fo(bo f,ll a,ll b=oo,ll c=1){if(b==oo)b=a,a=0;rr tuple{a,b-f,c};} I(J T)T cei(T x,T y){rr y<0?x=-x,y=-y:x<=0?x/y:(x-1)/y+1;}I(J T)T abs(T x){rr x<0?-x:x;}u1 io(ll n,ll x=0){u1 a;fo(i,n)a.eb(i+x);rr a;} I(J T)au&O+=(v1&v,cs v1&u){fo(i,len(v))v[i]+=u[i];rr v;}I(J T)au O+(v1v,cs v1&u){rr v+=u;}I(J T)constexpr bo isv=0; I(J T)au&O-=(v1&v,cs v1&u){fo(i,len(v))v[i]-=u[i];rr v;}I(J T)au O-(v1v,cs v1&u){rr v-=u;}I(J T)constexpr bo isv> =1; I(J T,J U)au&O+=(v1&v,U c){fe(v,e)e+=c;rr v;}I(J T,J U)au O+(v1v,U c){rr v+=c;}I(J T)au&O++(v1&v){rr v+=1;} I(J T,J U)au&O-=(v1&v,U c){fe(v,e)e-=c;rr v;}I(J T,J U)au O-(v1v,U c){rr v-=c;}I(J T)au&O--(v1&v){rr v-=1;} I(J T)au&O^=(v1&v,cs v1&u){fe(u,e)v.eb(e);rr v;}I(J T)au O^(v1v,cs v1&u){rr v^=u;} I(J T,J F=less<>)T&srt(T&a,F f={}){sort(all(a),cref(f));rr a;}I(J T)T&rv(T&a){reverse(all(a));rr a;} I(J T)T&rsr(T&a){srt(a),rv(a);rr a;}I(J T)T&uq(T&a){srt(a),a.erase(unique(all(a)),end(a));rr a;} I(J T)decltype(au)bg(T&a){assert(len(a));rr*begin(a);}I(J T)T pof(deque&q){T r=bg(q);q.pop_front();rr r;} I(J T)decltype(au)rb(T&a){assert(len(a));rr*rbegin(a);}I(J T)T pob(deque&q){T r=rb(q);q.pop_back();rr r;} I(J T)au pop(T&a){J T::value_type r;ifc(isv||is_same_v)r=rb(a),a.pop_back();else r=a.top(),a.pop();rr r;} I(J...A)au max(cs A&...a){rr max(initializer_list>{a...});}I(J T,J U=T)vo af(v1&a,U e={}){a.ep(begin(a),e);} I(J...A)au min(cs A&...a){rr min(initializer_list>{a...});}I(J T)vo df(T&a){a.erase(begin(a));} I(J T,J U)ll arglb(cs v1&a,cs U&b){rr lower_bound(all(a),b)-begin(a);}I(J T,J U)bo amax(T&a,cs U&b){rr a&a,cs U&b){rr upper_bound(all(a),b)-begin(a);}I(J T,J U)bo amin(T&a,cs U&b){rr a>b?a=b,1:0;} I(J T=ll,ll n,ll i=0)au vec(cs ll(&l)[n],T x={}){ifc(n==i+1)rr vector(l[i],x);else rr vector(l[i],vec(l,x));} I(J S,bo f(S&,cs S&),J V)au scan(cs V&v){cpr;ifc(!isv)r={v,1};else fe(v,x)if(au t=scan(x);t.b)!r.b?r=t,0:f(r.a,t.a);rr r;} I(J V)au recv(cs V&){ifc(isv)rr recv(V(1)[0]);else rr V();}I(J T)bo add(T&a,cs T&b){rr a+=b,1;} I(J V)au max(cs v1&v){rr scan(v).a;}I(J V)au sum(cs v1&v){rr scan(v).a;} I(J V)au min(cs v1&v){rr scan(v).a;}I(ll n,J...A)vo resize(cs ll(&l)[n],A&...a){((a=vec(l,recv(a))),...);} I(J V)vo fill(v1&v,decltype(recv(V()))x={}){fe(v,e)ifc(isv)fill(e,x);else e=x;} I(J V)au zt(v1a){ifc(isv)fe(a,e)e=zt(e);fo(i,len(a)-1)a[i+1]+=a[i];rr a;} I(J V)au mb(v1a){ifc(isv)fe(a,e)e=mb(e);of(i,len(a)-1)a[i+1]-=a[i];rr a;} struct fio{fio(){ios::sync_with_stdio(0);cin.tie(0);cout<=10)o<&v){str s=nl;ifc(!isv)s=sp;else ifc(isv)s+=nl;fe(v,e)o<&v){fe(v,e)o<&m){fe(m,e)o<&v){fe(v,e)o<>(is&i,v1&v){fe(v,e)i>>e;rr i;} I(J...A)vo li(A&...a){(cin>>...>>a);}I(J...A)vo vi(A&...a){fo(i,(len(a)&...))(cin>>...>>a[i]);} I(J...A)vo pp(cs A&...a){ll i=0;((cout<; graph tog(cs u2&a){graph g;fe(a,e){g.eb();fe(e,v)rb(g).eb(v);}rr g;}u2 tou(cs graph&g){u2 a;fe(g,e){a.eb();fe(e,v,w)rb(a).eb(v);}rr a;} I(J G)au ygi(G&g,ll m,bo f=1){ifc(is_same_v){UV(m,a,b);a-=f,b-=f;fo(i,m)g[a[i]].eb(b[i]);rr tuple{a,b};} else{UV(m,a,b,c);a-=f,b-=f;fo(i,m)g[a[i]].eb(b[i],c[i]);rr tuple{a,b,c};}} I(J G)au mgi(G&g,ll m,bo f=1){ifc(is_same_v){au[a,b]=ygi(g,m,f);fo(i,m)g[b[i]].eb(a[i]);rr tuple{a,b};} else{au[a,b,c]=ygi(g,m,f);fo(i,m)g[b[i]].eb(a[i],c[i]);rr tuple{a,b,c};}}I(J G)au ti(G&g,bo f=1){rr mgi(g,len(g)-1,f);} u1 en(cs str&s,cs str&b="a"){u1 r;fe(s,e)r.eb(len(b)==1?e-b[0]:b.find_first_of(e));rr r;}u1 sien(cs str&b="a"){ST(s);rr en(s,b);} u2 en(cs strs&s,cs str&b="a"){u2 r;fe(s,e)r.eb(en(e,b));rr r;}u2 sien(ll n,cs str&b="a"){SS(n,s);rr en(s,b);} str de(cs u1&a,cs str&b="a"){str r;fe(a,e)r+=(len(b)==1?e+b[0]:0<=e&&e(f)){}I(J...A)decltype(au)O()(A&&...a)cs{rr F::O()(*this,std::forward(a)...);}}; I(J...A)u1 zp(A&...a){u1 v=(a^...);uq(v);([&](u1&u){fe(u,e)e=arglb(v,e);}(a),...);rr v;}I(J...A)vo dec(A&...a){((--a),...);} I(J F,J...A)u1 vs(F f,A&...a){u1 o=io((len(a)&...));srt(o,cref(f));([&](au&e){au c=e;fo(i,len(e))e[i]=c[o[i]];}(a),...);rr o;} I(J F)ll bsl(F f,ll o,ll n){wh(abs(o-n)>1)(f((o+n)/2)?o:n)=(o+n)/2;rr o;}ll sqr(ll x){ll r=max(sqrt(x)-1,0.);wh(r+1<=x/(r+1))++r;rr r;} I(J F)dd bsd(F f,dd o,dd n){fo(t,l2(abs(o-n)/ee)+5)(f((o+n)/2)?o:n)=(o+n)/2;rr o;} cls rle(cs u1&a){cls r;fe(a,e)len(r)&&e==rb(r).a?++rb(r).b:r.eb(e,1).b;rr r;}cls rce(u1 a){rr rle(srt(a));} i128 pow(i128 x,u128 n,u128 m=0){i128 r=1;wh(n)n&1?r*=x:r,x*=x,m?r%=m,x%=m:r,n>>=1;rr r;}I(J T)T sq(T a){rr a*a;} u1 divs(ll n){u1 r;for(ll i=1;i*i<=n;++i)if(n%i==0)r.eb(i),r.eb(n/i);rr uq(r);}I(J T)T cb(T a){rr a*a*a;} cls fact(ll n){umapm;for(ll i=2;i*i<=n;++i)wh(n%i==0)++m[i],n/=i;if(n>1)++m[n];cls r;fe(m,e)r.eb(e);rr srt(r);} I(J T)struct twelvefold{ v1fa,rf; twelvefold(ll n):fa(n+1,1),rf(n+1,1){fo(i,1,n+1)fa[i]=fa[i-1]*i;rf[n]=fa[n].inv();of(i,n)rf[i]=rf[i+1]*(i+1);} T O[](ll n){rr fa[n];}T O()(ll n,ll k){rr c(n,k);} T c(ll n,ll k){rr n<0?pwm1(k)*c(-n+k-1,k):k<0||n; vo main(){ll T=1;fo(t,T)solve();} vo solve(){ LL(N,P); twf tw(N); ml ans=tw[N]; ml su=1; fo(k,N/P+1){ ans-=su*tw.rf[k]; su*=tw(N-P*k,P)*tw[P-1]; } pp(ans); }}