結果

問題 No.731 等差数列がだいすき
コンテスト
ユーザー te-sh
提出日時 2018-09-20 14:12:07
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 2 ms / 1,500 ms
コード長 824 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,630 ms
コンパイル使用メモリ 154,828 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-13 01:40:15
合計ジャッジ時間 3,445 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 18
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import std.algorithm, std.container, std.conv, std.math, std.range, std.typecons, std.stdio, std.string;

auto rdsp(){return readln.splitter;}
void pick(R,T)(ref R r,ref T t){t=r.front.to!T;r.popFront;}
void readV(T...)(ref T t){auto r=rdsp;foreach(ref v;t)pick(r,v);}
void readA(T)(size_t n,ref T[]t){t=new T[](n);auto r=rdsp;foreach(ref v;t)pick(r,v);}

void main()
{
  int n; readV(n);
  int[] a; readA(n, a);

  auto x = new real[](n), y = new real[](n);
  foreach (i; 0..n) {
    x[i] = i;
    y[i] = a[i];
  }

  auto ax = x.sum/n, ay = y.sum/n;
  x[] -= ax; y[] -= ay;

  auto sxy = 0.0L, sx = 0.0L;
  foreach (i; 0..n) {
    sxy += x[i]*y[i];
    sx += x[i]^^2;
  }

  auto d = sxy/sx, b = ay-ax*d;

  auto c = 0.0L;
  foreach (i; 0..n) c += (d*i+b-a[i])^^2;

  writefln("%.7f %.7f", b, d);
  writefln("%.7f", c);
}
0