結果

問題 No.2769 Number of Rhombi
ユーザー KeebyKeeby
提出日時 2024-06-13 04:18:06
言語 C++23
(gcc 12.3.0 + boost 1.83.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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 3 ms
6,816 KB
testcase_02 AC 2 ms
6,944 KB
testcase_03 AC 124 ms
19,280 KB
testcase_04 AC 118 ms
17,924 KB
testcase_05 AC 102 ms
17,620 KB
testcase_06 AC 128 ms
16,968 KB
testcase_07 AC 158 ms
19,164 KB
testcase_08 AC 159 ms
18,932 KB
testcase_09 AC 157 ms
19,052 KB
testcase_10 AC 184 ms
19,000 KB
testcase_11 AC 160 ms
19,068 KB
testcase_12 AC 160 ms
18,916 KB
testcase_13 AC 161 ms
19,040 KB
testcase_14 AC 162 ms
19,112 KB
testcase_15 AC 164 ms
18,976 KB
testcase_16 AC 165 ms
19,104 KB
testcase_17 AC 186 ms
19,008 KB
testcase_18 AC 161 ms
19,060 KB
testcase_19 AC 164 ms
19,004 KB
testcase_20 AC 160 ms
19,240 KB
testcase_21 AC 154 ms
18,992 KB
testcase_22 AC 158 ms
19,068 KB
testcase_23 AC 163 ms
19,036 KB
testcase_24 AC 2 ms
6,940 KB
testcase_25 AC 187 ms
18,940 KB
testcase_26 AC 166 ms
19,072 KB
testcase_27 AC 164 ms
19,044 KB
testcase_28 AC 167 ms
18,964 KB
testcase_29 AC 163 ms
18,928 KB
testcase_30 AC 160 ms
19,140 KB
testcase_31 AC 161 ms
18,932 KB
testcase_32 AC 183 ms
18,960 KB
testcase_33 AC 159 ms
19,144 KB
testcase_34 AC 163 ms
18,936 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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