結果
問題 | No.2923 Mayor's Job |
ユーザー | Reid |
提出日時 | 2024-10-12 15:43:29 |
言語 | C++23 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 34 ms / 2,000 ms |
コード長 | 4,493 bytes |
コンパイル時間 | 4,563 ms |
コンパイル使用メモリ | 273,208 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-10-12 15:43:35 |
合計ジャッジ時間 | 5,258 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
6,816 KB |
testcase_01 | AC | 2 ms
6,820 KB |
testcase_02 | AC | 2 ms
6,816 KB |
testcase_03 | AC | 16 ms
6,816 KB |
testcase_04 | AC | 17 ms
6,820 KB |
testcase_05 | AC | 16 ms
6,816 KB |
testcase_06 | AC | 15 ms
6,816 KB |
testcase_07 | AC | 18 ms
6,820 KB |
testcase_08 | AC | 19 ms
6,820 KB |
testcase_09 | AC | 17 ms
6,820 KB |
testcase_10 | AC | 19 ms
6,816 KB |
testcase_11 | AC | 26 ms
6,820 KB |
testcase_12 | AC | 34 ms
6,820 KB |
testcase_13 | AC | 6 ms
6,820 KB |
testcase_14 | AC | 2 ms
6,820 KB |
testcase_15 | AC | 1 ms
6,820 KB |
testcase_16 | AC | 2 ms
6,816 KB |
testcase_17 | AC | 2 ms
6,816 KB |
testcase_18 | AC | 2 ms
6,816 KB |
testcase_19 | AC | 2 ms
6,820 KB |
ソースコード
#include <bits/stdc++.h> #include<ext/pb_ds/assoc_container.hpp> #include<ext/pb_ds/tree_policy.hpp> #include<ext/pb_ds/tag_and_trait.hpp> using namespace std; using namespace __gnu_pbds; #define rep(i,n) for(int i = 0;i<(n);i++) #define per(i,n) for(int i = (n-1);i>=0;i--) #define nfor(i,b,e) for(int i = b;i<(e);i++) #define Sort(a) sort((a).begin(), (a).end()) #define Rev(a) reverse((a).begin(), (a).end()) #define all(...) std::begin(__VA_ARGS__), std::end(__VA_ARGS__) #define rall(...) std::rbegin(__VA_ARGS__), std::rend(__VA_ARGS__)//reverseしてからall #define endl '\n' #define YesNo(bool) if(bool){cout<<"Yes"<<endl;}else{cout<<"No"<<endl;} #include <limits> // Tu は最大値のある符号無し整数型 (型チェックに必要な static_assert などは省略) template <typename Tu> static inline constexpr bool is_multiplication_safe(Tu a, Tu b){ return !a || b <= std::numeric_limits<Tu>::max() / a; } template <typename Tu> static inline constexpr bool is_addition_safe(Tu a, Tu b){ if(a*b <= 0) return true; else if(a<0 && b<0) return b >= std::numeric_limits<Tu>::max() - a; else if(a>0 && b>0) return b <= std::numeric_limits<Tu>::max() - a; } //O(NlogN)の計算量が必要 template <typename T> double median(std::vector<T> numbers) { std::sort(numbers.begin(), numbers.end()); size_t size = numbers.size(); if (size % 2 == 0) return (numbers[size / 2 - 1] + numbers[size / 2]) / 2.0; else return numbers[size / 2]; } #ifdef _DEBUG //vectorの要素を列挙 #define debug(x) cout<<#x<<": "<<x<<'\n' template <typename T> void debugv(const vector<T>& vec) { for (const auto& element : vec) cout << element << " "; cout << endl; } //2次元vectorの要素を列挙 template <typename T> void debugvv(const vector<vector<T>>& vec) { for (const auto& e1 : vec) { for(const auto& e2:e1) cout << e2 << " "; cout << endl; } cout << endl; } #else #define debug(x) #define debugv(x) #define debugvv(x) #endif typedef long long ll; typedef unsigned long long ull; //<int>を<使いたい型>に書き変える //*orderd_set.find_by_order(i)で値を取得する。 using ordered_set = tree<char, null_type, less<char>, rb_tree_tag, tree_order_statistics_node_update>; using vi = vector<int>; using vc = vector<char>; using vvi = vector<vector<int>>; using vvc = vector<vector<char>>; using vvvi = vector<vector<vector<int>>>; using vll = vector<ll>; using vvll = vector<vector<ll>>; using Pii = pair<int,int>; #define make(x,y) make_pair(x,y) //vector to string string to_string(vector<int> a){ int n = a.size(); string res = ""; rep(i,n) res += char(a[i]+'0'); return res; } long double euclid_dist(long double xi,long double yi,long double xj,long double yj){ return sqrt((xi-xj)*(xi-xj)+(yi-yj)*(yi-yj)); } ull factorial(int n){ if(n==0||n==1) return 1ULL; else return factorial(n-1)*n; } template<int MOD> struct ModInt { static const int Mod = MOD; unsigned x; ModInt() : x(0) { } ModInt(signed sig) { int sigt = sig % MOD; if (sigt < 0) sigt += MOD; x = sigt; } ModInt(signed long long sig) { int sigt = sig % MOD; if (sigt < 0) sigt += MOD; x = sigt; } int get() const { return (int)x; } ModInt &operator+=(ModInt that) { if ((x += that.x) >= MOD) x -= MOD; return *this; } ModInt &operator-=(ModInt that) { if ((x += MOD - that.x) >= MOD) x -= MOD; return *this; } ModInt &operator*=(ModInt that) { x = (unsigned long long)x * that.x % MOD; return *this; } ModInt operator+(ModInt that) const { return ModInt(*this) += that; } ModInt operator-(ModInt that) const { return ModInt(*this) -= that; } ModInt operator*(ModInt that) const { return ModInt(*this) *= that; } }; using mint = ModInt<1000000007>; int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); //歴史がそれ未満かつ、近く(ユークリッド距離がk以下)なら消滅させる //どうすれば最小?どうやって消す? int n,k; cin >> n >> k; vll h(n);rep(i,n) cin >> h[i]; vector<Pii> xy(n); rep(i,n){ int x,y; cin >> x >> y; xy[i] = make(x,y); } vector<bool> ans(n); rep(i,n){ bool ok = true; rep(j,n){ if(euclid_dist(xy[j].first,xy[j].second,xy[i].first,xy[i].second)<=k && h[j]>h[i]){ ok = false; } } ans[i] = ok; } ll Ans = 0;rep(i,n) if(ans[i]) Ans++; cout << Ans << endl; return 0; }