#609126 (C++17) No.1361 [Zelkova 4th Tune *] QUADRUPLE-SEQUENCEの詩

提出ソース

結果

問題	No.1361 [Zelkova 4th Tune *] QUADRUPLE-SEQUENCEの詩
ユーザー	HIR180
提出日時	2021-01-22 23:28:58
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	RE
実行時間	-
コード長	3,438 bytes
コンパイル時間	4,437 ms
コンパイル使用メモリ	252,840 KB
実行使用メモリ	24,328 KB
最終ジャッジ日時	2024-06-08 18:28:31
合計ジャッジ時間	39,681 ms
ジャッジサーバーID （参考情報）	judge3 / judge4

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 6,812 KB
testcase_01	AC	2 ms 6,940 KB
testcase_02	AC	2 ms 6,940 KB
testcase_03	AC	2 ms 6,940 KB
testcase_04	AC	2 ms 6,944 KB
testcase_05	AC	2 ms 6,940 KB
testcase_06	AC	2 ms 6,940 KB
testcase_07	AC	2 ms 6,940 KB
testcase_08	AC	4 ms 6,944 KB
testcase_09	AC	2 ms 6,940 KB
testcase_10	AC	5 ms 6,944 KB
testcase_11	AC	4 ms 6,940 KB
testcase_12	AC	3 ms 6,944 KB
testcase_13	AC	6 ms 6,940 KB
testcase_14	AC	2 ms 6,940 KB
testcase_15	AC	7 ms 6,940 KB
testcase_16	AC	4 ms 6,940 KB
testcase_17	AC	6 ms 6,944 KB
testcase_18	AC	73 ms 6,940 KB
testcase_19	AC	28 ms 6,944 KB
testcase_20	AC	80 ms 6,944 KB
testcase_21	AC	20 ms 6,940 KB
testcase_22	AC	35 ms 6,940 KB
testcase_23	AC	53 ms 6,944 KB
testcase_24	AC	41 ms 6,940 KB
testcase_25	AC	49 ms 6,944 KB
testcase_26	AC	64 ms 6,940 KB
testcase_27	AC	40 ms 6,940 KB
testcase_28	RE	-
testcase_29	AC	446 ms 8,044 KB
testcase_30	AC	409 ms 8,168 KB
testcase_31	AC	543 ms 12,524 KB
testcase_32	AC	381 ms 14,220 KB
testcase_33	AC	373 ms 7,920 KB
testcase_34	AC	586 ms 13,620 KB
testcase_35	AC	749 ms 12,772 KB
testcase_36	AC	1,056 ms 15,788 KB
testcase_37	AC	188 ms 6,940 KB
testcase_38	AC	576 ms 9,540 KB
testcase_39	AC	340 ms 7,660 KB
testcase_40	AC	415 ms 8,028 KB
testcase_41	AC	548 ms 10,564 KB
testcase_42	AC	1,063 ms 15,964 KB
testcase_43	AC	563 ms 11,588 KB
testcase_44	AC	770 ms 14,120 KB
testcase_45	AC	651 ms 10,000 KB
testcase_46	AC	535 ms 12,840 KB
testcase_47	AC	296 ms 7,272 KB
testcase_48	AC	1,075 ms 21,888 KB
testcase_49	AC	954 ms 23,044 KB
testcase_50	AC	1,125 ms 24,324 KB
testcase_51	AC	1,073 ms 23,444 KB
testcase_52	AC	877 ms 22,820 KB
testcase_53	AC	1,483 ms 24,248 KB
testcase_54	AC	1,465 ms 23,552 KB
testcase_55	AC	1,580 ms 23,384 KB
testcase_56	AC	1,488 ms 23,596 KB
testcase_57	AC	1,230 ms 23,760 KB
testcase_58	AC	1,106 ms 23,616 KB
testcase_59	AC	1,112 ms 23,360 KB
testcase_60	AC	900 ms 21,872 KB
testcase_61	AC	908 ms 23,152 KB
testcase_62	AC	913 ms 22,092 KB
testcase_63	RE	-
testcase_64	RE	-
testcase_65	RE	-
testcase_66	RE	-
testcase_67	RE	-
testcase_68	RE	-
testcase_69	RE	-
testcase_70	AC	2 ms 6,940 KB
testcase_71	AC	2 ms 6,944 KB
testcase_72	AC	1 ms 6,940 KB
testcase_73	RE	-

権限があれば一括ダウンロードができます

ソースコード

raw source code

//Let's join Kaede Takagaki Fan Club !!
#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")
#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
using namespace std;
//#define int long long
typedef long long ll;
typedef pair<int,int> P;
typedef pair<int,P> P1;
typedef pair<P,P> P2;
#define pu push
#define pb push_back
#define eb emplace_back
#define mp make_pair
#define eps 1e-7
#define INF 1000000000
#define a first
#define b second
#define fi first
#define sc second
#define rng(i,a,b) for(int i=int(a);i<int(b);i++)
#define rep(i,x) for(int i=0;i<x;i++)
#define repn(i,x) for(int i=1;i<=x;i++)
#define SORT(x) sort(x.begin(),x.end())
#define ERASE(x) x.erase(unique(x.begin(),x.end()),x.end())
#define POSL(x,v) (lower_bound(x.begin(),x.end(),v)-x.begin())
#define POSU(x,v) (upper_bound(x.begin(),x.end(),v)-x.begin())
#define all(x) x.begin(),x.end()
#define si(x) int(x.size())

#ifdef LOCAL
#define dmp(x) cerr<<__LINE__<<" "<<#x<<" "<<x<<endl
#else
#define dmp(x) void(0)
#endif

template<class t,class u> bool chmax(t&a,u b){if(a<b){a=b;return true;}else return false;}
template<class t,class u> bool chmin(t&a,u b){if(b<a){a=b;return true;}else return false;}

template<class t> using vc=vector<t>;

template<class t,class u>
ostream& operator<<(ostream& os,const pair<t,u>& p){
	return os<<"{"<<p.fi<<","<<p.sc<<"}";
}

template<class t> ostream& operator<<(ostream& os,const vc<t>& v){
	os<<"{";
	for(auto e:v)os<<e<<",";
	return os<<"}";
}

template<class T>
void g(T &a){
	cin >> a;
}
template<class T>
void o(const T &a,bool space=false){
	cout << a << (space?' ':'\n');
}
//ios::sync_with_stdio(false);
const ll mod = 1000000007;//998244353
//mt19937_64 mt(chrono::steady_clock::now().time_since_epoch().count());
template<class T>
void add(T&a,T b){
	a+=b;
	if(a >= mod) a-=mod;
}
ll modpow(ll x,ll n){
	ll res=1;
	while(n>0){
		if(n&1) res=res*x%mod;
		x=x*x%mod;
		n>>=1;
	}
	return res;
}
ll divv(ll a, ll b) { // floored division
	if (b < 0) a *= -1, b *= -1;
	if (a >= 0) return a / b;
	return -((-a + b - 1) / b);
}

int n[4];
vc<int>a[4];
ll cnt;
vc<pair<ll,P>>v[2];

void solve(){
	rep(i, 4) cin >> n[i];
	cin >> cnt;
	rep(i, 4){
	    a[i].resize(n[i]);
		rep(j, n[i]) cin >> a[i][j];
	}
	rep(i, 2){
		for(auto at:a[i*2]){
			for(auto at2:a[i*2+1]){
				v[i].pb(mp(1LL*at*at2, mp(at, at2)));
			}
		}
		SORT(v[i]);
	}
	ll lb = -1e18, ub = 1e18;
	//(lb, ub]
	while(ub-lb > 1){
		ll mid = (lb+ub)/2, C = 0;
		for(auto att:v[0]){
			ll at = att.a; 
			ll lv, uv;
			if(at == 0){
				if(mid < 0) continue;
				lv = -1e18, uv = 1e18;
			}
			else if(at > 0){
				lv = -1e18;
				uv = divv(mid, at);
			}
			else{
				uv = 1e18;
				lv = divv(mid, at);
				if(mid%at != 0) lv ++;
			}
			ll X = v[1].size(), Y = 0;
			if(uv != 1e18) X = POSU(v[1], mp(uv+1, mp(-INF, -INF)));
			if(lv != -1e18) Y = POSL(v[1], mp(lv, mp(-INF, -INF)));
			C += X - Y;
		}
		if(C >= cnt) ub = mid;
		else lb = mid;
	}
	for(auto att:v[0]){
		ll at = att.a;
		if(ub%at != 0) continue;
		int x = POSL(v[1], mp(ub/at, mp(-INF, -INF)));
		if(x < v[1].size() && v[1][x].a == ub/at){
		    cout<<ub<<endl;
			cout << att.b.a << " " << att.b.b << " ";
			cout << v[1][x].b.a << " " << v[1][x].b.b << endl;
			return;
		}
	}
	assert(false);
}
signed main(){
	cin.tie(0);
	ios::sync_with_stdio(0);
	cout<<fixed<<setprecision(20);
	int t; t = 1; //cin >> t;
	while(t--) solve();
}

yukicoder

結果

テストケース

ソースコード