結果

問題 No.75 回数の期待値の問題
コンテスト
ユーザー eitaho
提出日時 2015-01-30 05:31:28
言語 C#(csc)
(csc 3.9.0)
結果
AC  
実行時間 48 ms / 5,000 ms
コード長 3,039 bytes
コンパイル時間 2,317 ms
コンパイル使用メモリ 115,356 KB
実行使用メモリ 28,360 KB
最終ジャッジ日時 2024-06-23 04:25:24
合計ジャッジ時間 2,969 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 16
権限があれば一括ダウンロードができます
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc)
Copyright (C) Microsoft Corporation. All rights reserved.

ソースコード

diff #
raw source code 

using System;
using System.IO;
using System.Text;
using System.Linq;
using System.Diagnostics;
using System.Collections.Generic;
using Enu = System.Linq.Enumerable;

class Program
{
    void Solve()
    {
        int N = reader.Int();
        var mat = Enu.Range(0, N).Select(i => new double[N + 1]).ToArray();
        for (int r = 0; r < N; r++)
        {
            mat[r][r] = mat[r][N] = 6;
            for (int d = 1; d <= 6; d++)
            {
                int c = (r + d > N) ? 0 : r + d;
                if (c < N) mat[r][c]--;
            }
        }
        GaussJordan(mat);
        Console.WriteLine(mat[0][N]);
    }

    static bool GaussJordan(double[][] M)
    {
        double EPS = 1e-10;
        int R = M.Length;
        int C = M[0].Length;

        for (int i = 0; i < R; i++)
        {
            int pivot = i;
            for (int k = i + 1; k < R; k++)
                if (Math.Abs(M[k][i]) > Math.Abs(M[pivot][i]))
                    pivot = k;
            var t = M[pivot]; M[pivot] = M[i]; M[i] = t;
            if (Math.Abs(M[i][i]) < EPS) return false;
            for (int j = C - 1; j >= i; j--)
                M[i][j] /= M[i][i];
            for (int k = 0; k < R; k++)
                if (k != i)
                    for (int j = C - 1; j >= i; j--)
                        M[k][j] -= M[k][i] * M[i][j];
        }
        return true;
    }



    static void Main() { new Program().Solve(); }
    Reader reader = new Reader(Console.In);
    class Reader
    {
        private readonly TextReader reader;
        private readonly char[] separator = { ' ' };
        private readonly StringSplitOptions removeOp = StringSplitOptions.RemoveEmptyEntries;
        private string[] A = new string[0];
        private int i;

        public Reader(TextReader r) { reader = r; }
        public bool HasNext() { return Enqueue(); }
        public string String() { return Dequeue(); }
        public int Int() { return int.Parse(Dequeue()); }
        public long Long() { return long.Parse(Dequeue()); }
        public double Double() { return double.Parse(Dequeue()); }
        public int[] IntLine() { var s = Line(); return s == "" ? new int[0] : Array.ConvertAll(Split(s), int.Parse); }
        public int[] IntArray(int N) { return Enu.Range(0, N).Select(i => Int()).ToArray(); }
        public int[][] IntGrid(int H) { return Enu.Range(0, H).Select(i => IntLine()).ToArray(); }
        public string[] StringArray(int N) { return Enu.Range(0, N).Select(i => Line()).ToArray(); }
        public string Line() { return reader.ReadLine().Trim(); }
        private string[] Split(string s) { return s.Split(separator, removeOp); }
        private bool Enqueue()
        {
            if (i < A.Length) return true;
            string line = reader.ReadLine();
            if (line == null) return false;
            if (line == "") return Enqueue();
            A = Split(line);
            i = 0;
            return true;
        }
        private string Dequeue() { Enqueue(); return A[i++]; }
    }
}
0