結果
| 問題 |
No.691 E869120 and Constructing Array 5
|
| コンテスト | |
| ユーザー |
vjudge1
|
| 提出日時 | 2025-09-21 17:06:58 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 158 ms / 1,000 ms |
| コード長 | 987 bytes |
| コンパイル時間 | 1,780 ms |
| コンパイル使用メモリ | 167,764 KB |
| 実行使用メモリ | 7,716 KB |
| 最終ジャッジ日時 | 2025-09-21 17:07:05 |
| 合計ジャッジ時間 | 6,465 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 |
| other | AC * 27 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
using ld=long double;
int q;ld x;
int main() {
cin>>q;
for(int i=0;i<q;i++) {
cin>>x;
int sz=10;
int lim=x*x /(4*sz*sz)+1.01L;
for(int j=lim;;j++) {
ld bs=sqrt((ld)(j));
ld rem=sqrt((ld)(j))*(2*sz)-x;
vector<int> delta(sz);
for(int k=0;k<sz;k++) {
int l=0,r=j;
while(r-l>1) {
int m=(l+r)>>1;
if(rem>2*bs-sqrt((ld)(j-m))-sqrt((ld)(j+m))) l=m;
else r=m;
}
rem -= 2*bs-sqrt((ld)(j-l))-sqrt((ld)(j+l));
delta[k]=l;
}
if(rem<1.0e-10L) {
vector<ld> ret(2*sz);
for(int k=0;k<sz;k++) {
ret[2*k]=j-delta[k];
ret[2*k+1]=j+delta[k];
}
ld val=0.0;
for(int k=0;k<ret.size();k++) {
val += sqrt(ret[k]);
}
cerr << fixed << setprecision(15) << j << ' ' << lim << ' ' << val << ' ' << abs(val-x) << endl;
cout << ret.size();
for(int k=0;k<ret.size();k++) {
cout << ' ' << ret[k];
}
cout << '\n';
break;
}
}
}
return 0;
}
vjudge1