結果

問題 No.2769 Number of Rhombi
コンテスト
ユーザー Keeby
提出日時 2024-06-13 04:18:06
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 187 ms / 5,000 ms
コード長 3,363 bytes
コンパイル時間 3,251 ms
コンパイル使用メモリ 259,540 KB
実行使用メモリ 19,280 KB
最終ジャッジ日時 2024-06-13 04:18:18
合計ジャッジ時間 9,714 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 32
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
//#include <atcoder/modint>
using namespace std;
//using namespace atcoder;
//using mint = modint998244353;
#define rep(i,n) for(long long i=0;i<n;i++)
#define rep1(i,n) for(long long i=1;i<=n;i++)
#define Rep(i,a,b) for(long long i=a;i<=b;i++)
#define all(a) (a).begin(), (a).end()
#define fst first
#define snd second
typedef unsigned long long ull;
typedef long long ll;
typedef vector<int> vec;
typedef vector<vector<int>> vvec;
typedef vector<long long> vecll;
typedef vector<vector<long long>> vvecll;

//#include<boost/multiprecision/cpp_int.hpp>
//using namespace boost::multiprecision;

const ll mod=998244353;

class mint {
public:
    long long x;
    constexpr mint(long long x=0) : x((x%mod+mod)%mod) {}
    constexpr mint operator-() const { 
      return mint(-x);
    }
    constexpr mint& operator+=(const mint& a) {
        if ((x += a.x) >= mod) x -= mod;
        return *this;
    }
    constexpr mint& operator-=(const mint& a) {
        if ((x += mod-a.x) >= mod) x -= mod;
        return *this;
    }
    constexpr mint& operator*=(const mint& a) {
        (x *= a.x) %= mod;
        return *this;
    }
    constexpr mint operator+(const mint& a) const {
        mint res(*this);
        return res+=a;
    }
    constexpr mint operator-(const mint& a) const {
        mint res(*this);
        return res-=a;
    }
    constexpr mint operator*(const mint& a) const {
        mint res(*this);
        return res*=a;
    }
    constexpr mint pow(long long t) const {
        if (!t) return 1;
        mint a = pow(t>>1);
        a *= a;
        if (t&1) a *= *this;
        return a;
    }

    // for prime mod
    constexpr mint inv() const {
        return pow(mod-2);
    }
    constexpr mint& operator/=(const mint& a) {
        return (*this) *= a.inv();
    }
    constexpr mint operator/(const mint& a) const {
        mint res(*this);
        return res/=a;
    }
};
ostream& operator<<(ostream& os, const mint& m){
    os << m.x;
    return os;
}


//xのn乗をpで割ったあまりを求める
ll modpow(ll x,ll n,ll p=mod){
  if(n==0) return 1;
  if(n%2==1) return (x*modpow(x,n-1,p))%p;
  else{
    ll t=modpow(x,n/2,p);
    assert(t>0);
    return (t*t)%p;
  }
}

//xの逆元(p:素数)
ll inv(ll x,ll p=mod){
  return modpow(x,p-2,p);
}

template <typename T>
int index(vector<T> &a, T x){
	return lower_bound(all(a),x)-a.begin();
}

const ll mod2=1000000009;

const ll INF=5e18;
const int INT_INF=1e9;
const double pi=3.14159265358979;

void slope(pair<ll,ll> &k){
  ll g=gcd(k.fst,k.snd);
  k.fst/=g;k.snd/=g;
  if(k.fst<0){
    k.fst/=-1;k.snd/=-1;
  }else if(k.fst==0 && k.snd<0){
    k.snd*=-1;
  }
  return;
}

void rotate(pair<ll,ll> &k){
  swap(k.fst,k.snd);
  k.fst*=-1;
  slope(k);
  return;
}


int main(){
	int n;
  cin >> n;
  vecll x(n),y(n);
  rep(i,n){
    cin >> x[i] >> y[i];
    x[i]*=2;y[i]*=2;
  }
  vector<pair<pair<ll,ll>,pair<ll,ll>>> a,b;
  rep(i,n){
    Rep(j,i+1,n-1){
      pair<ll,ll> p,k;
      p=make_pair((x[i]+x[j])/2,(y[i]+y[j])/2);
      k=make_pair(x[i]-x[j],y[i]-y[j]);
      slope(k);
      if(k.snd>0){
        a.push_back({p,k});
      }else{
        rotate(k);
        b.push_back({p,k});
      }
    }
  }
  sort(all(a));sort(all(b));
  ll ans=0;
  for(auto p:a){
    ans+=upper_bound(all(b),p)-lower_bound(all(b),p);
  }
  cout << ans << endl;
	return 0;
}
0