結果

問題 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;
}
0