結果

問題 No.731 等差数列がだいすき
ユーザー ducktailducktail
提出日時 2018-09-09 14:09:58
言語 Haskell
(9.8.2)
結果
AC  
実行時間 2 ms / 1,500 ms
コード長 1,026 bytes
コンパイル時間 4,214 ms
コンパイル使用メモリ 180,912 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-02 11:01:52
合計ジャッジ時間 5,166 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 1 ms
5,376 KB
testcase_03 AC 1 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 1 ms
5,376 KB
testcase_06 AC 1 ms
5,376 KB
testcase_07 AC 1 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 1 ms
5,376 KB
testcase_10 AC 2 ms
5,376 KB
testcase_11 AC 2 ms
5,376 KB
testcase_12 AC 1 ms
5,376 KB
testcase_13 AC 1 ms
5,376 KB
testcase_14 AC 1 ms
5,376 KB
testcase_15 AC 2 ms
5,376 KB
testcase_16 AC 2 ms
5,376 KB
testcase_17 AC 2 ms
5,376 KB
testcase_18 AC 2 ms
5,376 KB
testcase_19 AC 2 ms
5,376 KB
testcase_20 AC 2 ms
5,376 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
Loaded package environment from /home/judge/.ghc/x86_64-linux-9.8.2/environments/default
[1 of 2] Compiling Main             ( Main.hs, Main.o )
[2 of 2] Linking a.out

ソースコード

diff # 

import Control.Applicative
import Data.ByteString.Char8 (ByteString)
import qualified Data.ByteString.Char8 as B
import Data.List
import Data.Char (isSpace)

main :: IO ()
main = solve <$> readLn <*> f >>= prt
  where
    f = readil B.readInt <$> B.getLine

solve :: Int -> [Int] -> (Double, Double, Double)
solve n as = (b, a, c)
  where
    xb = fromIntegral (n - 1) / 2.0
    yb = fromIntegral (sum as) / fromIntegral n
    dxs = map (\x -> fromIntegral x - xb) [0 .. n-1]
    dys = map (\y -> fromIntegral y - yb) as
    a = (sum $ zipWith (*) dxs dys) / (sum $ map (^2) dxs)
    b = yb - a * xb
    f x = a * x + b
    c = sum $ map (^2) $ zipWith (\u v -> u - fromIntegral v) (map (\x -> f (fromIntegral x)) [0 .. n-1]) as

prt :: (Double, Double, Double) -> IO ()
prt (a, b, c) = do
  putStrLn $ show a ++ " " ++ show b
  print c

readil :: Integral a =>  (ByteString -> Maybe (a, ByteString)) -> ByteString -> [a]
readil f = unfoldr g
  where
    g s = do
      (n, s') <- f s
      return (n, B.dropWhile isSpace s')
0