#include #include #define rep(i,b) for(int i=0;i=0;i--) #define rep1(i,b) for(int i=1;i using mpq = priority_queue, greater>; template bool chmax(T &a, const T &b) { if (a bool chmin(T &a, const T &b) { if (b ll sumv(const vector&a){ll res(0);for(auto&&x:a)res+=x;return res;} bool yn(bool a) { if(a) {cout << "Yes" << endl; return 1;} else {cout << "No" << endl; return 0;}} #define dame { cout << "No" << endl; return 0;} #define dame1 { cout << -1 << endl; return 0;} #define test(x) cout << "test" << x << endl; #define deb(x,y) cout << x << " " << y << endl; #define deb3(x,y,z) cout << x << " " << y << " " << z << endl; #define deb4(x,y,z,x2) cout << x << " " << y << " " << z << " " << x2 << endl; #define out cout << ans << endl; #define outv fore(yans , ans) cout << yans << "\n"; #define show(x) cerr<<#x<<" = "<; using pil = pair; using pli = pair; using pii = pair; using tp = tuple; using vi = vector; using vl = vector; using vs = vector; using vb = vector; using vpii = vector; using vpli = vector; using vpll = vector; using vpil = vector; using vvi = vector>; using vvl = vector>; using vvs = vector>; using vvb = vector>; using vvpii = vector>; using vvpli = vector>; using vvpll = vector; using vvpil = vector; using mint = modint998244353; using vm = vector; using vvm = vector>; vector dx={1,0,-1,0,1,1,-1,-1},dy={0,1,0,-1,1,-1,1,-1}; ll gcd(ll a, ll b) { return a?gcd(b%a,a):b;} ll lcm(ll a, ll b) { return a/gcd(a,b)*b;} const double eps = 1e-10; const ll LINF = 1001002003004005006ll; const int INF = 1001001001; int main(){ ll z,n; cin>>n>>z; mint ans = 0; vl x(n),y(n); rep(i,n) cin>>x[i]>>y[i]; vm vsum(n,1); mint sumz = 1; int zero = -1; rep(i,n){ rep(j,n){ if (i == j) continue; vsum[i] *= (x[i] - x[j]); } if (z == x[i]) zero = i; sumz *= (z - x[i]); } if (zero >= 0){ ans += mint(y[zero])*(n-1); rep(i,n){ if (i == zero) continue; mint tmp1 = y[i],tmp2 = 1; rep(j,n){ if (j == i || j == zero) continue; tmp1 *= (z - x[j]); tmp2 *= (x[i] - x[j]); } ans += tmp1 / tmp2; } cout << ans.val() << endl; return 0; } mint a=0,b=0; rep(i,n){ a += mint(1) / (z - x[i]); b += mint(x[i]) / (z - x[i]); } rep(i,n){ mint tmp = 1; tmp *= sumz; tmp *= (a*x[i] - b); tmp /= (vsum[i] * (z - x[i])); tmp *= y[i]; ans += tmp; } cout << ans.val() << endl; return 0; }