提出 #47937860


ソースコード 拡げる

#include <bits/stdc++.h>

using namespace std;

using ll = long long;
using pii = pair<int,int>;
using pll = pair<ll,ll>;
using vi = vector<int>;
using vl = vector<ll>;

using _loop_int = int;
#define REP(i,n) for(_loop_int i=0; i<(_loop_int)(n); i++)
#define FOR(i,a,b) for(_loop_int i=(_loop_int)(a); i<(_loop_int)(b); i++)
#define FORR(i,a,b) for(_loop_int i=(_loop_int)(b)-1; i>=(_loop_int)(a); i--)

#define CHMIN(a,b) (a)=min((a),(b))
#define CHMAX(a,b) (a)=max((a),(b))
#define ALL(v) (v).begin(),(v).end()

#define DEBUG(x) cerr<<#x<<": "<<(x)<<endl
#define DEBUG_VEC(v) cerr<<#v<<": ";REP(__i,(v).size())cerr<<((v)[__i])<<", ";cerr<<endl

const ll MOD = 1000000007ll;

int n;
ll k;
pii pts[252521];

int main(){
    #if 1
    scanf("%d%lld",&n,&k);
    REP(i,n){int x,y;scanf("%d%d",&x,&y);pts[i]=pii(x,y);}
    #else
    srand(830252);
    n = 200000;
    k = 4'000'000'000'000ll;
    REP(i,n){
        int x,y;
        x = rand() % 1'000'000'000;
        y = rand() % 1'000'000'000;
        pts[i] = pii(x,y);
    }
    #endif

    vl xs(n), ys(n);
    REP(i,n)xs[i]=pts[i].first,ys[i]=pts[i].second;
    sort(ALL(xs));sort(ALL(ys));
    vl yys(n);
    yys[0] = 0ll;  // 1*(y-ys[0]) + 0
    FOR(i,1,n){
        yys[i] = yys[i-1] + 1ll*(ys[i]-ys[i-1])*i;
    }
    vl yys2(n);
    yys2[n-1] = 0ll;
    FORR(i,0,n-1){
        yys2[i] = yys2[i+1] + 1ll*(ys[i+1]-ys[i])*(n-1-i);
    }
    vl xxs(n);
    xxs[0] = 0ll;  // 1*(x-xs[0]) + 0
    FOR(i,1,n){
        xxs[i] = xxs[i-1] + 1ll*(xs[i]-xs[i-1])*i;
    }
    vl xxs2(n);
    xxs2[n-1] = 0ll;
    FORR(i,0,n-1){
        xxs2[i] = xxs2[i+1] + 1ll*(xs[i+1]-xs[i])*(n-1-i);
    }

    auto fx = [&](){
        ll ok = 1'000'000'005, ng=-1;
        while(abs(ok-ng)>1) {
            ll len = (ok+ng)/2;

            ll xmin = k+1ll;
            ll ymin = k+1ll;
            {
                const vl& ss = xs;
                const vl& ssa = xxs;
                const vl& ssb = xxs2;
                {
                    int top = 0;
                    REP(bot,n){
                        while(top+1<n && ss[top+1]-ss[bot] <= len) top++;
                        ll v = ssa[bot];
                        if (top+1<n) v += ssb[top+1] + 1ll*(ss[top+1]-ss[bot]-len)*(n-1-top);
                        CHMIN(xmin,v);
                    }
                }
                {
                    int bot = n-1;
                    FORR(top,0,n){
                        while(bot-1>=0 && ss[top]-ss[bot-1] <= len) bot--;
                        ll v = ssb[top];
                        if (bot-1>=0) v += ssa[bot-1] + 1ll*(ss[top]-len-ss[bot-1])*bot;
                        CHMIN(xmin,v);
                    }
                }
            }
            {
                const vl& ss = ys;
                const vl& ssa = yys;
                const vl& ssb = yys2;
                {
                    int top = 0;
                    REP(bot,n){
                        while(top+1<n && ss[top+1]-ss[bot] <= len) top++;
                        ll v = ssa[bot];
                        if (top+1<n) v += ssb[top+1] + 1ll*(ss[top+1]-ss[bot]-len)*(n-1-top);
                        CHMIN(ymin,v);
                    }
                }
                {
                    int bot = n-1;
                    FORR(top,0,n){
                        while(bot-1>=0 && ss[top]-ss[bot-1] <= len) bot--;
                        ll v = ssb[top];
                        if (bot-1>=0) v += ssa[bot-1] + 1ll*(ss[top]-len-ss[bot-1])*bot;
                        CHMIN(ymin,v);
                    }
                }
            }
            (xmin+ymin<=k ? ok : ng) = len;
        }
        return ok;
    };
    ll ans = fx();
    printf("%lld\n",ans);
    return 0;
}

提出情報

提出日時
問題 F - Minimize Bounding Square
ユーザ rickytheta
言語 C++ 20 (gcc 12.2)
得点 525
コード長 3852 Byte
結果 AC
実行時間 281 ms
メモリ 14436 KiB

コンパイルエラー

Main.cpp: In function ‘int main()’:
Main.cpp:31:10: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   31 |     scanf("%d%lld",&n,&k);
      |     ~~~~~^~~~~~~~~~~~~~~~
Main.cpp:32:27: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   32 |     REP(i,n){int x,y;scanf("%d%d",&x,&y);pts[i]=pii(x,y);}
      |                      ~~~~~^~~~~~~~~~~~~~

ジャッジ結果

セット名 Sample All
得点 / 配点 0 / 0 525 / 525
結果
AC × 3
AC × 104
セット名 テストケース
Sample sample_01.txt, sample_02.txt, sample_03.txt
All hack_01.txt, hack_02.txt, hack_03.txt, sample_01.txt, sample_02.txt, sample_03.txt, test_01.txt, test_02.txt, test_03.txt, test_04.txt, test_05.txt, test_06.txt, test_07.txt, test_08.txt, test_09.txt, test_10.txt, test_11.txt, test_12.txt, test_13.txt, test_14.txt, test_15.txt, test_16.txt, test_17.txt, test_18.txt, test_19.txt, test_20.txt, test_21.txt, test_22.txt, test_23.txt, test_24.txt, test_25.txt, test_26.txt, test_27.txt, test_28.txt, test_29.txt, test_30.txt, test_31.txt, test_32.txt, test_33.txt, test_34.txt, test_35.txt, test_36.txt, test_37.txt, test_38.txt, test_39.txt, test_40.txt, test_41.txt, test_42.txt, test_43.txt, test_44.txt, test_45.txt, test_46.txt, test_47.txt, test_48.txt, test_49.txt, test_50.txt, test_51.txt, test_52.txt, test_53.txt, test_54.txt, test_55.txt, test_56.txt, test_57.txt, test_58.txt, test_59.txt, test_60.txt, test_61.txt, test_62.txt, test_63.txt, test_64.txt, test_65.txt, test_66.txt, test_67.txt, test_68.txt, test_69.txt, test_70.txt, test_71.txt, test_72.txt, test_73.txt, test_74.txt, test_75.txt, test_76.txt, test_77.txt, test_78.txt, test_79.txt, test_80.txt, test_81.txt, test_82.txt, test_83.txt, test_84.txt, test_85.txt, test_86.txt, test_87.txt, test_88.txt, test_89.txt, test_90.txt, test_91.txt, test_92.txt, test_93.txt, test_94.txt, test_95.txt, test_96.txt, test_97.txt, test_98.txt
ケース名 結果 実行時間 メモリ
hack_01.txt AC 79 ms 14308 KiB
hack_02.txt AC 79 ms 14320 KiB
hack_03.txt AC 79 ms 14436 KiB
sample_01.txt AC 1 ms 3832 KiB
sample_02.txt AC 1 ms 3644 KiB
sample_03.txt AC 1 ms 3712 KiB
test_01.txt AC 1 ms 3636 KiB
test_02.txt AC 1 ms 3820 KiB
test_03.txt AC 1 ms 3824 KiB
test_04.txt AC 1 ms 3764 KiB
test_05.txt AC 1 ms 3844 KiB
test_06.txt AC 1 ms 3884 KiB
test_07.txt AC 1 ms 3768 KiB
test_08.txt AC 1 ms 3580 KiB
test_09.txt AC 1 ms 3888 KiB
test_10.txt AC 1 ms 3680 KiB
test_11.txt AC 1 ms 3788 KiB
test_12.txt AC 1 ms 3760 KiB
test_13.txt AC 1 ms 3760 KiB
test_14.txt AC 1 ms 3760 KiB
test_15.txt AC 1 ms 3580 KiB
test_16.txt AC 1 ms 3764 KiB
test_17.txt AC 1 ms 3768 KiB
test_18.txt AC 1 ms 3864 KiB
test_19.txt AC 1 ms 3684 KiB
test_20.txt AC 1 ms 3788 KiB
test_21.txt AC 1 ms 3644 KiB
test_22.txt AC 1 ms 3884 KiB
test_23.txt AC 1 ms 3688 KiB
test_24.txt AC 1 ms 3756 KiB
test_25.txt AC 1 ms 3644 KiB
test_26.txt AC 1 ms 3796 KiB
test_27.txt AC 272 ms 14348 KiB
test_28.txt AC 8 ms 4080 KiB
test_29.txt AC 76 ms 6140 KiB
test_30.txt AC 233 ms 14324 KiB
test_31.txt AC 200 ms 12708 KiB
test_32.txt AC 112 ms 14176 KiB
test_33.txt AC 233 ms 14344 KiB
test_34.txt AC 97 ms 8784 KiB
test_35.txt AC 281 ms 14108 KiB
test_36.txt AC 132 ms 14324 KiB
test_37.txt AC 186 ms 12292 KiB
test_38.txt AC 166 ms 11180 KiB
test_39.txt AC 232 ms 14380 KiB
test_40.txt AC 28 ms 5860 KiB
test_41.txt AC 12 ms 4572 KiB
test_42.txt AC 274 ms 14324 KiB
test_43.txt AC 257 ms 13784 KiB
test_44.txt AC 18 ms 4560 KiB
test_45.txt AC 232 ms 14336 KiB
test_46.txt AC 180 ms 11936 KiB
test_47.txt AC 27 ms 4608 KiB
test_48.txt AC 114 ms 14332 KiB
test_49.txt AC 133 ms 11136 KiB
test_50.txt AC 22 ms 4596 KiB
test_51.txt AC 251 ms 14344 KiB
test_52.txt AC 190 ms 14032 KiB
test_53.txt AC 164 ms 11144 KiB
test_54.txt AC 232 ms 14296 KiB
test_55.txt AC 165 ms 11120 KiB
test_56.txt AC 108 ms 13784 KiB
test_57.txt AC 231 ms 14188 KiB
test_58.txt AC 94 ms 8800 KiB
test_59.txt AC 170 ms 11364 KiB
test_60.txt AC 233 ms 14388 KiB
test_61.txt AC 218 ms 11692 KiB
test_62.txt AC 200 ms 12736 KiB
test_63.txt AC 232 ms 14188 KiB
test_64.txt AC 85 ms 11388 KiB
test_65.txt AC 100 ms 8272 KiB
test_66.txt AC 224 ms 14348 KiB
test_67.txt AC 116 ms 8788 KiB
test_68.txt AC 111 ms 7680 KiB
test_69.txt AC 232 ms 14292 KiB
test_70.txt AC 111 ms 11432 KiB
test_71.txt AC 170 ms 11404 KiB
test_72.txt AC 113 ms 14308 KiB
test_73.txt AC 45 ms 5588 KiB
test_74.txt AC 16 ms 4280 KiB
test_75.txt AC 68 ms 14344 KiB
test_76.txt AC 68 ms 14312 KiB
test_77.txt AC 79 ms 14188 KiB
test_78.txt AC 81 ms 14328 KiB
test_79.txt AC 80 ms 14268 KiB
test_80.txt AC 80 ms 14432 KiB
test_81.txt AC 80 ms 14296 KiB
test_82.txt AC 80 ms 14268 KiB
test_83.txt AC 159 ms 14320 KiB
test_84.txt AC 145 ms 14292 KiB
test_85.txt AC 145 ms 14184 KiB
test_86.txt AC 145 ms 14340 KiB
test_87.txt AC 145 ms 14404 KiB
test_88.txt AC 90 ms 14308 KiB
test_89.txt AC 139 ms 14316 KiB
test_90.txt AC 145 ms 14184 KiB
test_91.txt AC 153 ms 14332 KiB
test_92.txt AC 156 ms 14308 KiB
test_93.txt AC 153 ms 14428 KiB
test_94.txt AC 155 ms 14332 KiB
test_95.txt AC 154 ms 14332 KiB
test_96.txt AC 90 ms 14312 KiB
test_97.txt AC 119 ms 14400 KiB
test_98.txt AC 120 ms 14432 KiB