hdu4865(概率DP)

2024-09-09 17:18
文章标签 dp 概率 hdu4865

本文主要是介绍hdu4865(概率DP),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题意:已知前一天和今天的天气概率,某天的天气概率和叶子的潮湿程度的概率,n天叶子的湿度,求n天最有可能的天气情况。

思路:概率DP,dp[i][j]表示第i天天气为j的概率,状态转移如下:dp[i][j] = max(dp[i][j, dp[i-1][k]*table2[k][j]*table1[j][col] ) 

代码如下:

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <map>
#include <queue>
#include <vector>
#include <string>
#include <stack>
#include <iostream>
#define N 55
using namespace std;
double dp[N][5];
char s[10];
int pre[N][5];
double table1[3][4] = { {0.6,0.2,0.15,0.05}, { 0.25,0.3,0.2,0.25 }, { 0.05,0.1,0.35,0.5 } };
double table2[3][3] = { {0.5,0.375,0.125}, {0.25,0.125,0.625}, {0.25,0.375,0.375} };
int get_col(char *s)
{if(strcmp(s,"Dry") == 0)    return 0;if(strcmp(s,"Dryish") == 0)  return 1;if(strcmp(s,"Damp") == 0)   return 2;return 3;
}
map<int,string> mp;
int main()
{mp[0] = "Sunny";mp[1] = "Cloudy";mp[2] = "Rainy";int t,cas = 1;scanf("%d",&t);while(t--){int n;scanf("%d",&n);memset(dp,0,sizeof(dp));int i,j,k;scanf("%s",s);int col = get_col(s);dp[1][0] = 0.63*table1[0][col];dp[1][1] = 0.17*table1[1][col];dp[1][2] = 0.20*table1[2][col];for(i = 2; i <= n; i++){scanf("%s",s);col = get_col(s);for(j = 0; j < 3; j++)for(k = 0; k < 3; k++){double temp = dp[i-1][k]*table2[k][j]*table1[j][col];if(dp[i][j] < temp){dp[i][j] = temp;pre[i][j] = k;}}}double temp = 0;for(i = 0; i < 3; i++)if(dp[n][i] > temp){k = i;temp = dp[n][i];}// cout<<k<<endl;//cout<<temp<<endl;//for(i = 1; i <= n; i++){// for(j = 0; j < 3; j++)//       cout<<dp[i][j]<<" ";//     cout<<endl;}//   for(i = 1; i <= n; i++){//         for(j = 0; j < 3; j++)//               cout<<pre[i][j]<<" ";
//           / cout<<endl;}printf("Case #%d:\n",cas++);stack<string> q;while(n){//        printf("%s\n",mp[k]);//  cout<<mp[k]<<endl;q.push(mp[k]);k = pre[n][k];//          cout<<k<<endl;n--;}while(!q.empty())   {cout<<q.top()<<endl;q.pop();}}return 0;
}


这篇关于hdu4865(概率DP)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1151789

相关文章

hdu4826(三维DP)

这是一个百度之星的资格赛第四题 题目链接:http://acm.hdu.edu.cn/contests/contest_showproblem.php?pid=1004&cid=500 题意:从左上角的点到右上角的点,每个点只能走一遍,走的方向有三个:向上,向下,向右,求最大值。 咋一看像搜索题,先暴搜,TLE,然后剪枝,还是TLE.然后我就改方法,用DP来做,这题和普通dp相比,多个个向上

hdu1011(背包树形DP)

没有完全理解这题, m个人,攻打一个map,map的入口是1,在攻打某个结点之前要先攻打其他一个结点 dp[i][j]表示m个人攻打以第i个结点为根节点的子树得到的最优解 状态转移dp[i][ j ] = max(dp[i][j], dp[i][k]+dp[t][j-k]),其中t是i结点的子节点 代码如下: #include<iostream>#include<algorithm

usaco 1.1 Broken Necklace(DP)

直接上代码 接触的第一道dp ps.大概的思路就是 先从左往右用一个数组在每个点记下蓝或黑的个数 再从右到左算一遍 最后取出最大的即可 核心语句在于: 如果 str[i] = 'r'  ,   rl[i]=rl[i-1]+1, bl[i]=0 如果 str[i] = 'b' ,  bl[i]=bl[i-1]+1, rl[i]=0 如果 str[i] = 'w',  bl[i]=b

uva 10154 DP 叠乌龟

题意: 给你几只乌龟,每只乌龟有自身的重量和力量。 每只乌龟的力量可以承受自身体重和在其上的几只乌龟的体重和内。 问最多能叠放几只乌龟。 解析: 先将乌龟按力量从小到大排列。 然后dp的时候从前往后叠,状态转移方程: dp[i][j] = dp[i - 1][j];if (dp[i - 1][j - 1] != inf && dp[i - 1][j - 1] <= t[i]

uva 10118 dP

题意: 给4列篮子,每次从某一列开始无放回拿蜡烛放入篮子里,并且篮子最多只能放5支蜡烛,数字代表蜡烛的颜色。 当拿出当前颜色的蜡烛在篮子里存在时,猪脚可以把蜡烛带回家。 问最多拿多少只蜡烛。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cs

uva 10069 DP + 大数加法

代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <queue>#include <map>#include <cl

uva 10029 HASH + DP

题意: 给一个字典,里面有好多单词。单词可以由增加、删除、变换,变成另一个单词,问能变换的最长单词长度。 解析: HASH+dp 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

XTU 1233 n个硬币连续m个正面个数(dp)

题面: Coins Problem Description: Duoxida buys a bottle of MaiDong from a vending machine and the machine give her n coins back. She places them in a line randomly showing head face or tail face o

dp算法练习题【8】

不同二叉搜索树 96. 不同的二叉搜索树 给你一个整数 n ,求恰由 n 个节点组成且节点值从 1 到 n 互不相同的 二叉搜索树 有多少种?返回满足题意的二叉搜索树的种数。 示例 1: 输入:n = 3输出:5 示例 2: 输入:n = 1输出:1 class Solution {public int numTrees(int n) {int[] dp = new int

状态dp总结

zoj 3631  N 个数中选若干数和(只能选一次)<=M 的最大值 const int Max_N = 38 ;int a[1<<16] , b[1<<16] , x[Max_N] , e[Max_N] ;void GetNum(int g[] , int n , int s[] , int &m){ int i , j , t ;m = 0 ;for(i = 0 ;