#pragma GCC optimize("Ofast") #pragma GCC target("avx2") #include using namespace std; #ifdef LOCAL #include "debug.h" #else #define DEBUG(...) #endif template > constexpr T power(T a, long long n, Op op = Op(), T e = {1}) { assert(n >= 0); while (n) { if (n & 1) e = op(e, a); if (n >>= 1) a = op(a, a); } return e; } template struct modular { using m = modular; static constexpr unsigned mod = M; unsigned v; modular(long long x = 0) : v((x %= mod) < 0 ? x + mod : x) {} m operator-() const { return m() -= *this; } m& operator+=(m b) { if ((int)(v += b.v - mod) < 0) v += mod; return *this; } m& operator-=(m b) { if ((int)(v -= b.v) < 0) v += mod; return *this; } m& operator*=(m b) { v = (uint64_t)v * b.v % mod; return *this; } m& operator/=(m b) { return *this *= power(b, mod - 2); } friend m operator+(m a, m b) { return a += b; } friend m operator-(m a, m b) { return a -= b; } friend m operator*(m a, m b) { return a *= b; } friend m operator/(m a, m b) { return a /= b; } friend bool operator==(m a, m b) { return a.v == b.v; } }; using mint = modular; int x[2000], y[2000], a[2000], b[2000]; int main() { cin.tie(nullptr); ios::sync_with_stdio(false); int n; cin >> n; for (int i = 0; i < n; ++i) { cin >> x[i] >> y[i]; } mint res; for (int i = 0; i < n; ++i) { for (int j = 0; j < n; ++j) { a[j] = x[j] - x[i]; b[j] = y[j] - y[i]; } long long cur = 0; for (int j = i + 1; j < n; ++j) { for (int k = j + 1; k < n; ++k) { cur += abs(a[j] * b[k] - b[j] * a[k]); } } res += cur; } cout << res.v << '\n'; }