結果

問題 No.1495 パンの仕入れ
ユーザー 👑 potato167potato167
提出日時 2022-02-14 21:38:29
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 242 ms / 2,000 ms
コード長 2,946 bytes
コンパイル時間 3,884 ms
コンパイル使用メモリ 207,124 KB
実行使用メモリ 18,320 KB
最終ジャッジ日時 2023-09-11 16:48:13
合計ジャッジ時間 11,332 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 2 ms
4,504 KB
testcase_02 AC 160 ms
4,380 KB
testcase_03 AC 160 ms
4,380 KB
testcase_04 AC 160 ms
4,376 KB
testcase_05 AC 159 ms
4,376 KB
testcase_06 AC 160 ms
4,380 KB
testcase_07 AC 162 ms
4,376 KB
testcase_08 AC 162 ms
4,380 KB
testcase_09 AC 162 ms
4,376 KB
testcase_10 AC 162 ms
4,380 KB
testcase_11 AC 163 ms
4,380 KB
testcase_12 AC 50 ms
4,380 KB
testcase_13 AC 49 ms
4,376 KB
testcase_14 AC 50 ms
4,380 KB
testcase_15 AC 50 ms
4,376 KB
testcase_16 AC 51 ms
4,376 KB
testcase_17 AC 49 ms
4,376 KB
testcase_18 AC 49 ms
4,376 KB
testcase_19 AC 50 ms
4,376 KB
testcase_20 AC 50 ms
4,380 KB
testcase_21 AC 50 ms
4,380 KB
testcase_22 AC 242 ms
18,116 KB
testcase_23 AC 240 ms
18,316 KB
testcase_24 AC 36 ms
5,136 KB
testcase_25 AC 35 ms
5,204 KB
testcase_26 AC 36 ms
5,068 KB
testcase_27 AC 35 ms
5,180 KB
testcase_28 AC 224 ms
18,180 KB
testcase_29 AC 228 ms
18,320 KB
testcase_30 AC 40 ms
4,380 KB
testcase_31 AC 40 ms
4,376 KB
testcase_32 AC 41 ms
4,384 KB
testcase_33 AC 40 ms
4,376 KB
testcase_34 AC 40 ms
4,376 KB
testcase_35 AC 49 ms
4,380 KB
testcase_36 AC 49 ms
4,376 KB
testcase_37 AC 49 ms
4,376 KB
testcase_38 AC 48 ms
4,376 KB
testcase_39 AC 49 ms
4,380 KB
testcase_40 AC 37 ms
4,380 KB
testcase_41 AC 36 ms
4,376 KB
testcase_42 AC 36 ms
4,376 KB
testcase_43 AC 36 ms
4,376 KB
testcase_44 AC 36 ms
4,380 KB
testcase_45 AC 2 ms
4,380 KB
testcase_46 AC 3 ms
4,380 KB
testcase_47 AC 2 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#pragma GCC optimize("Ofast")
#define _GLIBCXX_DEBUG
using namespace std;
using std::cout;
using std::cin;
using std::endl;
using ll=long long;
using ld=long double;
ll ILL=1167167167167167167;
const int INF=2100000000;
const ll mod=1e9+7;
#define rep(i,a) for (int i=0;i<a;i++)
template<class T> using _pq = priority_queue<T, vector<T>, greater<T>>;
template<class T> ll LB(vector<T> &v,T a){return lower_bound(v.begin(),v.end(),a)-v.begin();}
template<class T> ll UB(vector<T> &v,T a){return upper_bound(v.begin(),v.end(),a)-v.begin();}
template<class T> bool chmin(T &a,const T &b){if(a>b){a=b;return 1;}else return 0;}
template<class T> bool chmax(T &a,const T &b){if(a<b){a=b;return 1;}else return 0;}
template<class T> void So(vector<T> &v) {sort(v.begin(),v.end());}
template<class T> void Sore(vector<T> &v) {sort(v.begin(),v.end(),[](T x,T y){return x>y;});}
void yneos(bool a){if(a) cout<<"Yes\n"; else cout<<"No\n";}



void solve();
//  rainy ~ 雨に打たれて ~
int main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	
	int t=1;
	cin>>t;
	rep(i,t) solve();
}

void solve(){
    ll N,M,K;
    cin>>N>>M>>K;
    vector<vector<ll>> G(N);
    rep(i,M){
        ll x,y;
        cin>>x>>y;
        G[x-1].push_back(y);
    }
    int J=0;
    ll bestsum=0;
    vector<ll> best(N),sum(N),sum2(N),size(N);
    rep(i,N){
        size[i]=G[i].size();
        if(size[i]==0){
            J=1;
            continue;
        }
        for(auto x:G[i]){
            sum[i]+=x;
            sum2[i]+=x*x;
        }
        best[i]=(2*sum[i]+size[i]-1)/(2*size[i]);
        bestsum+=best[i];
    }
    if(bestsum<=K){
        if(J){
            ll ans=0;
            rep(i,N){
                ans+=best[i]*(best[i]*size[i]-2*sum[i])+sum2[i];
            }
            cout<<ans<<"\n";
            return ;
        }
        ll l=-1,r=ILL;
        while(r-l>1){
            ll med=(r+l)/2;
            ll tmp=0;
            rep(i,N){
                if(tmp>K) break;
                if(size[i]) tmp+=(med+2*sum[i]-size[i])/(2*size[i])+1;
            }
            if(tmp>K) r=med;
            else l=med;
        }
        ll ans=0;
        ll c=0;
        rep(i,N){
            ll b=(l+2*sum[i]-size[i])/(2*size[i])+1;
            c+=b;
            ans+=b*(b*size[i]-2*sum[i])+sum2[i];
        }
        cout<<ans+r*(K-c)<<"\n";
    }else{
        ll l=-ILL,r=1;
        while(r-l>1){
            ll med=(r+l)/2;
            ll tmp=0;
            rep(i,N){
                if(tmp>K) break;
                if(size[i]) tmp+=max(0ll,(2*sum[i]+med+size[i]-1)/(2*size[i]));
            }
            if(tmp>K) r=med;
            else l=med;
        }
        ll ans=0;
        ll c=0;
        rep(i,N){
            if(size[i]==0) continue;
            ll b=max(0ll,(r+2*sum[i]+size[i]-1)/(2*size[i]));
            c+=b;
            ans+=b*(b*size[i]-2*sum[i])+sum2[i];
        }
        cout<<ans-l*(c-K)<<"\n";
    }
}
0