結果
| 問題 | No.1658 Product / Sum |
| コンテスト | |
| ユーザー |
沙耶花
|
| 提出日時 | 2021-08-27 21:54:13 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 14 ms / 2,000 ms |
| コード長 | 1,146 bytes |
| コンパイル時間 | 4,145 ms |
| コンパイル使用メモリ | 250,228 KB |
| 最終ジャッジ日時 | 2025-01-24 02:56:22 |
|
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 13 |
ソースコード
#include <stdio.h>
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace atcoder;
using mint = modint998244353;
using namespace std;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf 1000000001
pair<long long,long long> NG(-1,-1);
long long get(long long n){
long long ok = 0,ng = 1000000005;
while(ng-ok>1){
long long mid = (ok+ng)/2;
if(mid*mid<=n)ok = mid;
else ng = mid;
}
if(ok*ok!=n)return -1;
return ok;
}
pair<long long,long long> get(long long a,long long b){
long long t = a*a - b*4;
if(t<0)return NG;
t = get(t);
if(t<0)return NG;
pair<long long,long long> ret;
ret.first = a+t;
ret.second = a-t;
if(ret.first%2!=0)return NG;
ret.first /= 2;
ret.second /= 2;
//if(a==36)cout<<ret.first<<','<<ret.second<<endl;
if(min(ret.first,ret.second)<=0)return NG;
return ret;
}
int main(){
long long N,K;
cin>>N>>K;
vector<long long> ans(N,1);
for(int i=N;true;i++){
int S = i;
pair<long long,long long> p = get(S-(N-2),K*S);
if(p==NG)continue;
ans[0] = p.first;
ans[1] = p.second;
break;
}
rep(i,N){
if(i!=0)cout<<' ';
cout<<ans[i];
}
cout<<endl;
return 0;
}
沙耶花