Entries

スポンサーサイト (この記事を編集する[管理者用])

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

コメント

コメントの投稿

コメントの投稿
管理者にだけ表示を許可する

トラックバック

トラックバック URL
http://rsujskf.blog32.fc2.com/tb.php/2431-37adcf8e
この記事にトラックバックする(FC2ブログユーザー)

ARC 016 C - ソーシャルゲーム (この記事を編集する[管理者用])

Source

AtCoder Regular Contest #016
問題文

問題概要

カードが $N$ 種類,ガチャが $M$ 種類ある.
各ガチャを $1$ 回まわすコストと,どのカードが出るかという確率が全て与えられる.
全てのカードを集めるまでにかかるコストの平均値の最小値を求める問題.
$1 \leq N \leq 10$
$1 \leq M \leq 4$
部分点あり(原文を参照のこと)

解法

今どのカードが集まっているかを状態としてビットDPする.
既に持っているカードが出てくることがあるので,
 $X = pX + Y$

 $X = Y / (1-p)$
にするような式変形をする.

Cによるスパゲッティなソースコード
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<string.h>
#define REP(i,a,b) for(i=a;i<b;i++)
#define rep(i,n) REP(i,0,n)

double dp[3000];
int N, M;
int C[10], cost[10];
int idol[10][20], p[10][20];

double solve(int mask){
  int i, j, k;
  double res = 1e100, tmp, prob;

  if(dp[mask] > -1) return dp[mask];
  if(mask==0) return dp[mask] = 0;
  
  rep(k,M){
    tmp = 0;
    prob = 0;

    rep(i,C[k]){
      if(mask&1<<idol[k][i]){
        tmp += (solve(mask^(1<<idol[k][i])) + cost[k]) * p[k][i] / 100.0;
      } else {
        prob += p[k][i] / 100.0;
      }
    }

    if(prob > 1 - 1e-8) continue;
    tmp = (tmp + prob * cost[k]) / (1-prob);

    if(res > tmp) res = tmp;
  }

  return dp[mask] = res;
}

int main(){
  int i,j,k,l;
  double res;

  scanf("%d%d",&N,&M);
  rep(i,M){
    scanf("%d%d",C+i,cost+i);
    rep(j,C[i]){
      scanf("%d%d",idol[i]+j,p[i]+j);
      idol[i][j]--;
    }
  }

  rep(i,1<<N) dp[i] = -100;
  res = solve((1<<N)-1);

  printf("%.10f\n",res);

  return 0;
}

コメント

コメントの投稿

コメントの投稿
管理者にだけ表示を許可する

トラックバック

トラックバック URL
http://rsujskf.blog32.fc2.com/tb.php/2431-37adcf8e
この記事にトラックバックする(FC2ブログユーザー)

Appendix

Recent Articles

ブログ内検索

Ads


(プライバシーポリシー)
上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。