結果

問題 No.1597 Matrix Sort
ユーザー bayashiko
提出日時 2021-07-09 21:54:13
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 386 ms / 1,500 ms
コード長 1,821 bytes
コンパイル時間 2,212 ms
コンパイル使用メモリ 198,104 KB
最終ジャッジ日時 2025-01-22 21:36:26
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 27
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#pragma GCC optimize("Ofast")
//#pragma GCC target("avx2")
#pragma GCC optimize("unroll-loops")
using namespace std;
//#include<boost/multiprecision/cpp_int.hpp>
//#include<boost/multiprecision/cpp_dec_float.hpp>
//namespace mp=boost::multiprecision;
//#define mulint mp::cpp_int
//#define mulfloat mp::cpp_dec_float_100
struct __INIT{__INIT(){cin.tie(0);ios::sync_with_stdio(false);cout<<fixed<<setprecision(15);}} __init;
#define INF (1<<30)
#define LINF (lint)(1LL<<56)
#define endl "\n"
#define rep(i,n) for(lint (i)=0;(i)<(n);(i)++)
#define reprev(i,n) for(lint (i)=(n-1);(i)>=0;(i)--)
#define flc(x) __builtin_popcountll(x)
#define pint pair<int,int>
#define pdouble pair<double,double>
#define plint pair<lint,lint>
#define fi first
#define se second
#define all(x) x.begin(),x.end()
#define vec vector<lint>
#define nep(x) next_permutation(all(x))
typedef long long lint;
int dx[8]={1,1,0,-1,-1,-1,0,1};
int dy[8]={0,1,1,1,0,-1,-1,-1};
const int MAX_N=3e5+5;
template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return 1;}return 0;}
template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return 1;}return 0;}
//vector<int> bucket[MAX_N/1000];
constexpr int MOD=1000000007;
//constexpr int MOD=998244353;
/*#include<atcoder/all>
using namespace atcoder;
typedef __int128_t llint;*/

int main(void){
    lint N,K,P;
    cin >> N >> K >> P;
    int A[N],B[N];
    rep(i,N) cin >> A[i];
    rep(i,N) cin >> B[i];
    sort(A,A+N);
    sort(B,B+N);
    int lo=-1,hi=P;
    while(hi-lo>1){
        int mid=(hi+lo)/2;
        lint cnt=0;
        rep(i,N){
            cnt+=upper_bound(B,B+N,mid-A[i])-lower_bound(B,B+N,0-A[i]);
            cnt+=upper_bound(B,B+N,mid+P-A[i])-lower_bound(B,B+N,P-A[i]);
        }
        if(cnt>=K) hi=mid;
        else lo=mid;
    }
    cout << hi << endl;
}
0