結果

問題 No.738 平らな農地
ユーザー jelljell
提出日時 2018-10-27 17:41:49
言語 C++11
(gcc 11.4.0)
結果
CE  
(最新)
AC  
(最初)
実行時間 -
コード長 6,341 bytes
コンパイル時間 905 ms
コンパイル使用メモリ 95,884 KB
最終ジャッジ日時 2023-08-12 07:49:12
合計ジャッジ時間 1,978 ms
ジャッジサーバーID
(参考情報)
judge15 / judge13
このコードへのチャレンジ(β)
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。(リジャッジ後のコンパイルエラーは公開されます)
ただし、clay言語の場合は開発者のデバッグのため、公開されます。

コンパイルメッセージ
main.cpp: In function ‘uint64_t gcd(uint64_t, uint64_t)’:
main.cpp:110:59: error: call of overloaded ‘abs(uint64_t&)’ is ambiguous
 uint64_t gcd(uint64_t a, uint64_t b) { if(!b) return abs(a); if(a % b) return gcd(b, a % b); return b; }
                                                           ^
In file included from /usr/include/c++/8/cstdlib:75,
                 from /usr/include/c++/8/bits/stl_algo.h:59,
                 from /usr/include/c++/8/algorithm:62,
                 from main.cpp:10:
/usr/include/stdlib.h:837:12: note: candidate: ‘int abs(int)’
 extern int abs (int __x) __THROW __attribute__ ((__const__)) __wur;
            ^~~
In file included from /usr/include/c++/8/cstdlib:77,
                 from /usr/include/c++/8/bits/stl_algo.h:59,
                 from /usr/include/c++/8/algorithm:62,
                 from main.cpp:10:
/usr/include/c++/8/bits/std_abs.h:56:3: note: candidate: ‘long int std::abs(long int)’
   abs(long __i) { return __builtin_labs(__i); }
   ^~~
/usr/include/c++/8/bits/std_abs.h:61:3: note: candidate: ‘long long int std::abs(long long int)’
   abs(long long __x) { return __builtin_llabs (__x); }
   ^~~
/usr/include/c++/8/bits/std_abs.h:70:3: note: candidate: ‘constexpr double std::abs(double)’
   abs(double __x)
   ^~~
/usr/include/c++/8/bits/std_abs.h:74:3: note: candidate: ‘constexpr float std::abs(float)’
   abs(float __x)
   ^~~
/usr/include/c++/8/bits/std_abs.h:78:3: note: candidate: ‘constexpr long double std::abs(long double)’
   abs(long double __x)
   ^~~
/usr/include/c++/8/bits/std_abs.h:84:3: note: candidate: ‘constexpr __int128 std::abs(__int128)’
   abs(__GLIBCXX_TYPE_INT_N_0 __x) { return __x >= 0 ? __x : -__x; }
   ^~~
/usr/include/c++/8/bits/std_abs.h:102:3: note: candidate: ‘constexpr __float128 std::abs(__float128)’
   abs(__float128 __x)
   ^~~

ソースコード

diff #

/*
#include <boost/functional/hash.hpp>
#include <boost/multiprecision/cpp_int.hpp>
typedef boost::multiprecision::cpp_int mlint;
*/
// #include <bits/stdc++.h>
#include <deque>
#include <queue>
#include <vector>
#include <algorithm>
#include <iostream>
#include <set>
#include <cmath>
#include <tuple>
#include <string>
#include <functional>
#include <iterator>
#include <random>
#include <unordered_map>
#include <unordered_set>
#include <map>
#include <bitset>
#include <stack>
#include <iomanip>

using namespace std;
typedef int64_t i64;
typedef double db;
typedef vector<vector<int64_t>> mat;
typedef pair<int,int> pii;
typedef pair<int64_t,int> pli;
typedef pair<db,int> pdi;

#define esc(x) cout << (x) << endl,exit(0)
#define inf(T) (numeric_limits<T>::max() / 2 - 1)
#define each(i,v) for(auto i = begin(v),end_copy = end(v); i != end_copy; ++i)
#define reach(i,v) for(auto i = rbegin(v),rend_copy = rend(v); i != rend_copy; ++i)
#define urep(i,s,t) for(int64_t i = (s),t_copy = (t); i <= t_copy; ++i)
#define drep(i,s,t) for(int64_t i = (s),t_copy = (t); i >= t_copy; --i)
#define rep(i,n) urep(i,0,(n) - 1)
#define rep1(i,n) urep(i,1,(n))
#define rrep(i,n) drep(i,(n) - 1,0)
#define all(v) begin(v),end(v)
#define rall(v) rbegin(v),rend(v)
#define vct vector
#define prque priority_queue
#define u_map unordered_map
#define u_set unordered_set
#define l_bnd lower_bound
#define u_bnd upper_bound
#define rsz resize
#define ers erase
#define emp emplace
#define emf emplace_front
#define emb emplace_back
#define pof pop_front
#define pob pop_back
#define mkp make_pair
#define mkt make_tuple
#define fir first
#define sec second

void setup(uint8_t prec = 10) {
    cin.tie(0);
    ios::sync_with_stdio(false);
    cout << fixed << setprecision(prec);

    #ifdef Local

    #define debug 1
    #define print(x) cout << #x << " = " << (x) << " (L" << __LINE__ << ")" << endl

    cout << "-- Execution At Local ---\n" << "\n<Standard Output>\n\n";
    auto print_atexit = []() {
        time_t end_time = time(NULL);
        struct tm *ret = localtime(&end_time);
        cout << "\n----------------\n";
        cout << "\nSuccessfully Executed At " << (ret -> tm_hour) << ":" << (ret -> tm_min) << ":" << (ret -> tm_sec) << "\n\n";
    };
    atexit(print_atexit);

    #endif
}

template <class T, class U> ostream& operator << (ostream& s, pair<T,U> p) { return s << p.fir << " " << p.sec; }
template <class T> ostream& operator << (ostream& s, vector<T>& v) { for(auto i = v.begin(); i != v.end(); i++) { if(i != begin(v)) s << " "; s << *i; } return s; }

template <class T> void hash_combine(size_t &seed, T const &key) {
    hash<T> hasher;
    seed ^= hasher(key) + 0x9e3779b9 + (seed << 6) + (seed >> 2);
}

namespace std {
    template <class T, class U> struct hash<pair<T,U>> {
        size_t operator()(pair<T,U> const &p) const {
            size_t seed(0);
            hash_combine(seed,p.first);
            hash_combine(seed,p.second);
            return seed;
        }
    };
}

bool odd(i64 x) { return x & 1; }
bool even(i64 x) { return !odd(x); }
bool parity(i64 x, i64 y) { return odd(x) ^ even(y); }
bool bit(i64 n, uint8_t e) { return n >> e & 1; }
int ilog(i64 x, uint64_t b) { if(x) return ilog(x / b,b) + 1; return -1; }
i64 qceil(i64 x, i64 y) { return x > 0 ? (x - 1) / y + 1 :  x / y; }
uint64_t gcd(uint64_t a, uint64_t b) { if(!b) return abs(a); if(a % b) return gcd(b, a % b); return b; }
uint64_t lcm(uint64_t a, uint64_t b) { if(a & b) return a / gcd(a, b) * b; return 0; }
template <class T, class U> bool chmax(T& m, U x) { if(m < x) { m = x; return 1; } return 0; }
template <class T, class U> bool chmin(T& m, U x) { if(m > x) { m = x; return 1; } return 0; }
template <class T> T cmprs(T& v) {
    T tmp = v,ret = v;
    sort(all(tmp));
    tmp.erase(unique(all(tmp)),end(tmp));
    each(i,ret) *i = l_bnd(all(tmp),*i) - begin(tmp) + 1;
    return ret;
}

const int dir[8][2] = { {1,0},{0,1},{-1,0},{0,-1},{1,1},{-1,1},{-1,-1},{1,-1} };
const i64 mod = 1e9 + 7;
const db eps = 1e-15;


template <class Abel = int> struct BIT {
	typedef function<Abel(const Abel&,const Abel&)> opr_t;
	const opr_t opr,invopr;
	const Abel idel;
	const int sz;
	int range = 1;
	vector<Abel> data,val;

	BIT(int sz, Abel init_val, const opr_t &f = plus<Abel>(), const opr_t &g = minus<Abel>(), Abel idel = 0) : sz(sz), opr(f),invopr(g),idel(idel)
	{
		while(sz >= range) range <<= 1;
		data.resize(range);
		init(init_val);
	}

	void init(Abel v) {
		val.assign(range,v);
		val[0] = idel;
		vector<Abel> tmp(range);
		partial_sum(begin(val),end(val),begin(tmp),opr);
		for(int i = 0; i < range; ++i) data[i] = invopr(tmp[i],tmp[i - (i & -i)]);
	}

	Abel query(int n) { return n ? opr(query(n - (n & -n)),data[n]) : idel; }

	Abel query(int l, int r) { return invopr(query(r),query(l - 1)); }

	void add(int idx, Abel diff = 1) { update(idx,val[idx] + diff); }

	void update(int idx, Abel v) {
		Abel prev = val[idx];
		for(val[idx] = v; idx < range; idx += idx & -idx) data[idx] = opr(invopr(data[idx],prev),v);
	}

	int lower_bound(Abel obj, int l = 0, int r = 0) {
		if(!r) r = range;
		if(r - l <= 1) return min(r,sz + 1);
		int mid = (l + r) / 2;
		if(data[mid] >= obj) return lower_bound(obj,l,mid);
		return lower_bound(obj - data[mid],mid,r);
	}

	int upper_bound(Abel obj, int l = 0, int r = 0) {
		if(!r) r = range;
		if(r - l <= 1) return min(r,sz + 1);
		int mid = (l + r) / 2;
		if(data[mid] > obj) return upper_bound(obj,l,mid);
		return upper_bound(obj - data[mid],mid,r);
	}
};

int N,K;
i64 ans = inf(i64);

int main() {
    setup();

    cin >> N >> K;
    BIT<i64> bt(N,0),bt2(N,0);
    vct<pii> v(N);
    vct<int> a(N);
    rep(i,N) {
        cin >> v[i].fir;
        v[i].sec = i;
        a[i] = v[i].fir;
    }
    sort(all(v));
    vct<int> rank(N);
    rep(i,N) rank[v[i].sec] = i + 1;
    rep(i,K - 1) bt.add(rank[i]);
    rep(i,K - 1) bt2.add(rank[i],a[i]);
    urep(i,K,N) {
        bt.add(rank[i - 1]);
        bt2.add(rank[i - 1],a[i - 1]);
        if(odd(K)) {
            int lb = bt.l_bnd((K + 1) / 2);
            chmin(ans,bt2.query(lb,N) - bt2.query(lb));
        } else {
            int lb = bt.l_bnd(K / 2);
            chmin(ans,bt2.query(lb + 1,N) - bt2.query(lb));
        }
        bt.add(rank[i - K],-1);
        bt2.add(rank[i - K],-a[i - K]);
    }
    esc(ans);
}
0