#include using namespace std; using ll = long long; #define rep(i, s, t) for (int i = (int)(s); i < (int)(t); ++i) #define revrep(i, t, s) for (int i = (int)(t)-1; i >= (int)(s); --i) #define all(x) begin(x), end(x) template bool chmax(T& a, const T& b) { return a < b ? (a = b, 1) : 0; } template bool chmin(T& a, const T& b) { return a > b ? (a = b, 1) : 0; } template class Modint { using mint = Modint; static_assert(mod > 0, "Modulus must be positive"); public: static constexpr int get_mod() noexcept { return mod; } constexpr Modint(long long y = 0) noexcept : x(y >= 0 ? y % mod : (y % mod + mod) % mod) {} constexpr int value() const noexcept { return x; } constexpr mint& operator+=(const mint& r) noexcept { if ((x += r.x) >= mod) x -= mod; return *this; } constexpr mint& operator-=(const mint& r) noexcept { if ((x += mod - r.x) >= mod) x -= mod; return *this; } constexpr mint& operator*=(const mint& r) noexcept { x = static_cast(1LL * x * r.x % mod); return *this; } constexpr mint& operator/=(const mint& r) noexcept { *this *= r.inv(); return *this; } constexpr mint operator-() const noexcept { return mint(-x); } constexpr mint operator+(const mint& r) const noexcept { return mint(*this) += r; } constexpr mint operator-(const mint& r) const noexcept { return mint(*this) -= r; } constexpr mint operator*(const mint& r) const noexcept { return mint(*this) *= r; } constexpr mint operator/(const mint& r) const noexcept { return mint(*this) /= r; } constexpr bool operator==(const mint& r) const noexcept { return x == r.x; } constexpr bool operator!=(const mint& r) const noexcept { return x != r.x; } constexpr mint inv() const noexcept { int a = x, b = mod, u = 1, v = 0; while (b > 0) { int t = a / b; std::swap(a -= t * b, b); std::swap(u -= t * v, v); } return mint(u); } constexpr mint pow(long long n) const noexcept { mint ret(1), mul(x); while (n > 0) { if (n & 1) ret *= mul; mul *= mul; n >>= 1; } return ret; } friend std::ostream& operator<<(std::ostream& os, const mint& r) { return os << r.x; } friend std::istream& operator>>(std::istream& is, mint& r) { long long t; is >> t; r = mint(t); return is; } private: int x; }; using mint = Modint<998244353>; int main() { ios_base::sync_with_stdio(false); cin.tie(nullptr); cout << fixed << setprecision(15); int N; mint X; cin >> N >> X; vector x(N), y(N); rep(i,0,N) cin >> x[i] >> y[i]; vector prodinv(N, 1); rep(i,0,N) { mint p = 1; rep(k,0,N) if (k != i) p *= x[i] - x[k]; prodinv[i] = mint(1)/p; } mint ans = 0; rep(i,0,N) { mint left = 1; rep(k,0,i) left *= X - x[k]; vector right(N, 1); revrep(k,N-1,i) right[k] = right[k+1] * (X - x[k+1]); rep(j,i+1,N) { ans += (x[j] - x[i]) * (y[j] * prodinv[j] - y[i] * prodinv[i]) * left * right[j]; left *= X - x[j]; } } cout << ans << endl; }