結果

問題 No.76 回数の期待値で練習
ユーザー femtofemto
提出日時 2015-04-02 02:35:21
言語 C#(csc)
(csc 3.9.0)
結果
AC  
実行時間 84 ms / 5,000 ms
コード長 4,995 bytes
コンパイル時間 2,388 ms
コンパイル使用メモリ 110,156 KB
実行使用メモリ 29,984 KB
最終ジャッジ日時 2023-09-17 00:52:36
合計ジャッジ時間 2,983 ms
ジャッジサーバーID
(参考情報)
judge14 / judge12
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 81 ms
29,984 KB
testcase_01 AC 84 ms
29,964 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc)
Copyright (C) Microsoft Corporation. All rights reserved.

ソースコード

diff #

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

class Solution {
    int T;
    int N = 1000000;
    double[] p, dp;
    double[] ans = new double[]
    { 
        1.0833333333333333, 
        1.2569444444444444, 
        1.5353009259259260, 
        1.6915991512345676,
        2.0513639724794235 
    };

    void calc() {
        p = new double[6];
        dp = new double[N + 10];

        for (int k = 0; k < 5; k++) {
            double min = 0, max = 1, mid = 0;
            for (int i = 0; i < k; i++)
                max -= p[i];

            for (int loop = 0; loop < 100; loop++) {
                mid = (min + max) / 2;
                p[k] = mid;
                p[5] = 1;
                for (int i = 0; i < 5; i++)
                    p[5] -= p[i];
                f(k + 2);
                if (dp[k + 2] > ans[k])
                    max = mid;
                else
                    min = mid;
            }
            p[k] = mid;
            p[5] = 1;
            for (int i = 0; i < 5; i++)
                p[5] -= p[i];
        }

        f(N + 1);

        T = ri();
        var anslist = new List<double>();
        for (int loop = 0; loop < T; loop++) {
            var r = ri();
            anslist.Add(dp[r]);
        }

        foreach (var a in anslist)
            Console.WriteLine(a);
    }

    void f(int n) {
        dp[0] = 0;
        for (int i = 1; i < n + 1; i++) {
            dp[i] = 1;
            for (int j = 0; j < 6; j++) {
                if (i - j - 1 >= 0)
                    dp[i] += dp[i - j - 1] * p[j];
            }
        }
    }

    static void Main(string[] args) {
        new Solution().calc();
    }

    #region
    static int ri() { return int.Parse(Console.ReadLine()); }
    static int[] ria(int n) {
        if (n <= 0) { Console.ReadLine(); return new int[0]; }
        else return Console.ReadLine().Trim().Split(' ').Select(int.Parse).ToArray();
    }
    static void rio(out int p1) { p1 = ri(); }
    static void rio(out int p1, out int p2) { var r = ria(2); p1 = r[0]; p2 = r[1]; }
    static void rio(out int p1, out int p2, out int p3) { var r = ria(3); p1 = r[0]; p2 = r[1]; p3 = r[2]; }
    static void rio(out int p1, out int p2, out int p3, out int p4) { var r = ria(4); p1 = r[0]; p2 = r[1]; p3 = r[2]; p4 = r[3]; }
    static void rio(out int p1, out int p2, out int p3, out int p4, out int p5) { var r = ria(5); p1 = r[0]; p2 = r[1]; p3 = r[2]; p4 = r[3]; p5 = r[4]; }
    static long rl() { return long.Parse(Console.ReadLine()); }
    static long[] rla(int n) {
        if (n <= 0) { Console.ReadLine(); return new long[0]; }
        else return Console.ReadLine().Trim().Split(' ').Select(long.Parse).ToArray();
    }
    static void rlo(out long p1) { p1 = rl(); }
    static void rlo(out long p1, out long p2) { var r = rla(2); p1 = r[0]; p2 = r[1]; }
    static void rlo(out long p1, out long p2, out long p3) { var r = rla(3); p1 = r[0]; p2 = r[1]; p3 = r[2]; }
    static void rlo(out long p1, out long p2, out long p3, out long p4) { var r = rla(4); p1 = r[0]; p2 = r[1]; p3 = r[2]; p4 = r[3]; }
    static void rlo(out long p1, out long p2, out long p3, out long p4, out long p5) { var r = rla(5); p1 = r[0]; p2 = r[1]; p3 = r[2]; p4 = r[3]; p5 = r[4]; }
    static double rd() { return double.Parse(Console.ReadLine()); }
    static double[] rda(int n) {
        if (n <= 0) { Console.ReadLine(); return new double[0]; }
        else return Console.ReadLine().Trim().Split(' ').Select(double.Parse).ToArray();
    }
    static void rdo(out double p1) { p1 = rd(); }
    static void rdo(out double p1, out double p2) { var r = rda(2); p1 = r[0]; p2 = r[1]; }
    static void rdo(out double p1, out double p2, out double p3) { var r = rda(3); p1 = r[0]; p2 = r[1]; p3 = r[2]; }
    static void rdo(out double p1, out double p2, out double p3, out double p4) { var r = rda(4); p1 = r[0]; p2 = r[1]; p3 = r[2]; p4 = r[3]; }
    static void rdo(out double p1, out double p2, out double p3, out double p4, out double p5) { var r = rda(5); p1 = r[0]; p2 = r[1]; p3 = r[2]; p4 = r[3]; p5 = r[4]; }
    static void swap<T>(ref T x, ref T y) { T temp = x; x = y; y = temp; }
    static void wa1<T>(T[] a) { Debug.WriteLine(string.Join(" ", a)); }
    static void wa2<T>(T[][] a) {
        foreach (var row in a) {
            Debug.WriteLine(String.Join(" ", row));
        }
    }
    [DebuggerDisplay("{x} , {y}")]
    class point {
        public int x, y;
        public point(int x, int y) {
            this.x = x; this.y = y;
        }
    }
    #endregion
}

static class Extention {
    public static T[][] ToJagArray<T>(this T[,] a) {
        int n = a.GetLength(0), m = a.GetLength(1);
        var ret = new T[n][];
        for (int i = 0; i < n; i++) {
            ret[i] = new T[m];
            for (int j = 0; j < m; j++) {
                ret[i][j] = a[i, j];
            }
        }
        return ret;
    }
}

0