#line 2 "/Users/kaage/Desktop/ProgrammingWorkspace/library/other/template.hpp" #define _CRT_SECURE_NO_WARNINGS #pragma target("avx2") #pragma optimize("O3") #pragma optimize("unroll-loops") #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define rep(i,n) for(int i=0;i<(n);i++) #define REP(i,n) for(int i=1;i<=(n);i++) #define all(V) V.begin(),V.end() typedef long long lint; typedef unsigned long long ulint; typedef std::pair P; typedef std::pair LP; constexpr int INF = INT_MAX/2; constexpr lint LINF = LLONG_MAX/2; constexpr double eps = DBL_EPSILON; constexpr double PI=3.141592653589793238462643383279; template class prique :public std::priority_queue, std::greater> {}; template inline bool chmax(T& lhs, const U& rhs) { if (lhs < rhs) { lhs = rhs; return 1; } return 0; } template inline bool chmin(T& lhs, const U& rhs) { if (lhs > rhs) { lhs = rhs; return 1; } return 0; } inline lint gcd(lint a, lint b) { while (b) { lint c = a; a = b; b = c % b; } return a; } inline lint lcm(lint a, lint b) { return a / gcd(a, b) * b; } bool isprime(lint n) { if (n == 1)return false; for (int i = 2; i * i <= n; i++) { if (n % i == 0)return false; } return true; } template T mypow(T a, lint b) { T res(1); while(b){ if(b&1)res*=a; a*=a; b>>=1; } return res; } lint modpow(lint a, lint b, lint m) { lint res(1); while(b){ if(b&1){ res*=a;res/=m; } a*=a;a/=m; b>>=1; } return res; } template void printArray(std::vector& vec) { rep(i, vec.size()){ std::cout << vec[i]; std::cout<<(i==(int)vec.size()-1?"\n":" "); } } template void printArray(T l, T r) { T rprev = std::prev(r); for (T i = l; i != rprev; i++) { std::cout << *i << " "; } std::cout << *rprev << std::endl; } #line 3 "/Users/kaage/Desktop/ProgrammingWorkspace/library/algebraic/ModInt.hpp" class ModInt { lint value; public: static const unsigned int modulo; ModInt() : value(0) {} template ModInt(T value = 0) : value(value) { if (value < 0)value = -(lint)(-value % modulo) + modulo; this->value = value % modulo; } inline ModInt inv()const{return mypow(*this,modulo-2);} inline operator int()const { return value; } inline ModInt& operator+=(const ModInt& x) { value += x.value; if (value >= modulo)value -= modulo; return *this; } inline ModInt& operator++() { if (value == modulo - 1)value = 0; else value++; return *this; } inline ModInt operator-()const { return ModInt(0) -= *this; } inline ModInt& operator-=(const ModInt& x) { value -= x.value; if (value < 0)value += modulo; return *this; } inline ModInt& operator--() { if (value == 0)value = modulo - 1; else value--; return *this; } inline ModInt& operator*=(const ModInt& x) { value = value * x.value % modulo; return *this; } inline ModInt& operator/=(const ModInt& rhs) { return *this*=rhs.inv(); } template ModInt operator+(const T& rhs)const { return ModInt(*this) += rhs; } template ModInt& operator+=(const T& rhs) { return operator+=(ModInt(rhs)); } template ModInt operator-(const T& rhs)const { return ModInt(*this) -= rhs; } template ModInt& operator-=(const T& rhs) { return operator-=(ModInt(rhs)); } template ModInt operator*(const T& rhs)const { return ModInt(*this) *= rhs; } template ModInt& operator*=(const T& rhs) { return operator*=(ModInt(rhs)); } template ModInt operator/(const T& rhs)const { return ModInt(*this) /= rhs; } template ModInt& operator/=(const T& rhs) { return operator/=(ModInt(rhs)); } }; std::istream& operator>>(std::istream& ist, ModInt& x) { lint a; ist >> a; x = a; return ist; } #line 3 "main.cpp" const unsigned int ModInt::modulo=998244353; int main(){ lint n,k,l; std::cin>>n>>k>>l; if(n