結果

問題 No.162 8020運動
コンテスト
ユーザー eitaho
提出日時 2015-03-08 19:50:43
言語 C#(csc)
(csc 3.9.0)
結果
AC  
実行時間 1,554 ms / 5,000 ms
コード長 3,638 bytes
コンパイル時間 1,170 ms
コンパイル使用メモリ 116,720 KB
実行使用メモリ 141,048 KB
最終ジャッジ日時 2024-06-24 16:21:11
合計ジャッジ時間 38,809 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26
権限があれば一括ダウンロードができます
コンパイルメッセージ
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.IO;
using System.Linq;
using System.Text;
using System.Diagnostics;
using System.Collections.Generic;
using System.Text.RegularExpressions;
using Enu = System.Linq.Enumerable;

class Program
{
    static readonly int NumTeeth = 14;
    static readonly int Target = 80;

    void Solve()
    {
        int Age = reader.Int();
        var GoodProb = Enu.Range(0, 3).Select(i => 1 - reader.Double() / 100).ToArray();
        var state = new List<Edge>[1 << NumTeeth];

        for (int mask = 0; mask < (1 << NumTeeth); mask++)
        {
            state[mask] = new List<Edge>();
            int A = mask << 1;
            for (int sub = mask; ; sub = (sub - 1) & mask)
            {
                int B = sub << 1;
                double p = 1;
                for (int i = 1; i <= NumTeeth; i++)
                    if ((A >> i & 1) == 1)
                    {
                        int count = (A >> (i - 1) & 1) + (A >> (i + 1) & 1);
                        if ((B >> i & 1) == 1) p *= GoodProb[count];
                        else p *= 1 - GoodProb[count];
                    }
                state[mask].Add(new Edge(sub, p));
                if (sub == 0) break;
            }
        }

        double[] prob = new double[1 << NumTeeth];
        prob[(1 << NumTeeth) - 1] = 1;
        for (int age = Age; age < Target; age++)
        {
            var next = new double[prob.Length];
            for (int mask = 0; mask < (1 << NumTeeth); mask++)
                foreach (var e in state[mask])
                    next[e.next] += prob[mask] * e.prob;
            prob = next;
        }
        double ans = 0;
        for (int mask = 0; mask < prob.Length; mask++)
            ans += BitCount(mask) * prob[mask];
        ans *= 2;

        Console.WriteLine(ans);
    }

    struct Edge
    {
        public int next;
        public double prob;
        public Edge(int next, double p) { this.next = next; prob = p; }
    }

    int BitCount(long x) { return x == 0 ? 0 : 1 + BitCount(x & x - 1); }

    Reader reader = new Reader(Console.In);
    static void Main() { new Program().Solve(); }
}

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 Reader(string file) { reader = new StreamReader(file); }
    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