結果

問題 No.731 等差数列がだいすき
コンテスト
ユーザー bal4u
提出日時 2019-08-19 05:46:53
言語 C(gnu17)
(gcc 15.2.0)
コンパイル:
gcc-15 -O2 -std=gnu17 -Wno-error=implicit-function-declaration -Wno-error=implicit-int -Wno-error=incompatible-pointer-types -Wno-error=int-conversion -DONLINE_JUDGE -o a.out _filename_ -lm
実行:
./a.out
結果
AC  
実行時間 1 ms / 1,500 ms
コード長 895 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 165 ms
コンパイル使用メモリ 40,168 KB
最終ジャッジ日時 2026-02-22 04:22:15
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 18
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

// yukicoder: No.731 等差数列がだいすき
// bal4u 2019.8.19

#include <stdio.h>
#include <stdlib.h>

typedef long long ll;

//// 入出力関係
#if 1
#define gc() getchar_unlocked()
#else
#define gc() getchar()
#endif

int in() { // 非負整数の入力
	int n = 0, c = gc();
	do n = 10 * n + (c & 0xf); while ((c = gc()) >= '0');
	return n;
}

int N;
int y[1003];

int main()
{
	int i;
	ll sx, sx2, sy, sxy;
	ll bo;
	double a, b;   // y = ax + b
	double cost, c;
	
	sx = sx2 = sy = sxy = 0;
	
	N = in();
	sx = (ll)N*(N+1)/2, sx2 = sx*(2*N+1)/3;
	for (i = 1; i <= N; i++) {
		y[i] = in();
		sy += y[i], sxy += i*y[i];
	}
	
	bo = N*sx2 - sx*sx;
	a = (double)(N*sxy - sx*sy) / bo;
	b = (double)(sx2*sy - sxy*sx) / bo;
	
	printf("%.15lf %.15lf\n", a+b, a);
	
	cost = 0, c = a+b;
	for (i = 1; i <= N; i++) {
		cost += (c-y[i])*(c-y[i]);
		c += a;
	}
	printf("%.15lf\n", cost);
	return 0;
}
0