結果

問題 No.3027 f-列とh-列
ユーザー ジュ・ビオレ・グレイス
提出日時 2025-02-09 16:37:03
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,116 bytes
コンパイル時間 2,972 ms
コンパイル使用メモリ 81,892 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2025-02-09 20:39:04
合計ジャッジ時間 4,300 ms
ジャッジサーバーID
(参考情報)
judge1 / judge5
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

import std.array, std.conv, std.stdio;
immutable Nmax = 25;
alias Polynomial = long[Nmax+1];
Polynomial sum(Polynomial lhs, Polynomial rhs) {
foreach (i; 0 .. Nmax+1) {
lhs[i] += rhs[i];
}
return lhs;
}
Polynomial mul(Polynomial lhs, Polynomial rhs) {
Polynomial result;
foreach (i; 0 .. Nmax+1) {
foreach (j; 0 .. i+1) {
result[i] += lhs[j] * rhs[i-j];
}
}
return result;
}
Polynomial mul(long lhs, Polynomial rhs) {
foreach (i; 0 .. Nmax+1) {
rhs[i] *= lhs;
}
return rhs;
}
Polynomial solve(Polynomial f) {
Polynomial x_1;
x_1[0] = -1, x_1[1] = 1;
Polynomial x_1_power;
x_1_power[0] = 1;
Polynomial result;
foreach (i; 0 .. Nmax+1) {
result = sum( result, mul(f[i], x_1_power) );
x_1_power = mul(x_1, x_1_power);
}
return result;
}
void main() {
Polynomial f;
auto N = readln[0 .. $-1].to!ulong;
auto coeffs = readln.split.to!(long[]);
assert (coeffs.length == N+1 && N <= Nmax);
foreach (i; 0 .. N+1) {
f[N-i] = coeffs[i];
}
auto ans_pol = solve(f);
string ans;
foreach (i; 0 .. N+1) {
ans ~= ans_pol[N-i].to!string ~ " ";
}
ans[0 .. $-1].writeln;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0