BZOJ 4318 OSU! ( 期望DP )

最新推荐文章于 2024-08-15 05:02:54 发布
最新推荐文章于 2024-08-15 05:02:54 发布
阅读量85 收藏 1
点赞数
文章标签： 数据结构与算法
原文链接：http://www.cnblogs.com/LiHior/p/9551843.html
版权

题目链接

题意 : OSU 是一款群众喜闻乐见的休闲软件。 我们可以把 OSU 的规则简化与改编成以下的样子 : 一共有 n 次操作，每次操作只有成功与失败之分，成功对应 1 ，失败对应 0 ，n次操作对应为 1 个长度为 n 的 01 串。在这个串中连续的 X 个 1 可以贡献 X^3 的分数，这 X 个 1 不能被其他连续的 1 所包含（也就是极长的一串 1 ，具体见样例解释） 现在给出n，以及每个操作的成功率，请你输出期望分数，输出四舍五入后保留 1 位小数。 【样例说明】 000分数为0，001分数为1，010分数为1，100分数为1，101分数为2，110分数为8，011分数为8，111分数为27，总和为48，期望为48/8=6.0  ( N<=100000 )

 

分析 :

考虑 期望DP

首先得分是由连续 1 的长度决定的

而且每段连续 1 的贡献是相互独立的

那么考虑这样一个 len[i] = 以 i 位置为极长 1 的结尾的后缀期望长度

假设当前考虑到位置 i + 1

p[i] 为到 i 位置为 1 的概率

dp[i] 为以 i 位置为极长 1 的结尾的后缀期望得分

那么其得分期望 ( 即贡献 ) 有如下计算过程

dp[i + 1] = ( len[i] + 1 )^3 * p[i+1]

 

ans += dp[i+1] - dp[i] + 0 * ( 1 - p[i+1] )

       += dp[i+1] - dp[i] 

后面的 0 * (1 - p[i+1]) 为 i + 1 这个位为 0 的贡献 

前面的 dp[i+1] - dp[i] 为 i + 1 这个位为 1 的贡献

为什么是加 dp[i+1] - dp[i] 而不是 dp[i+1] 呢？

因为 dp[i+1] 本身就是由 dp[i] 递推而来

即 dp[i+1] = dp[i] + X

所以相当于 dp 数组是一个前缀期望得分和的形式

如果要得到每一位的贡献、当然是 dp[i+1] - dp[i]

每一位加起来就是总贡献

那么你可能会想能不能直接递推 dp[i+1] 就行了何必那么麻烦

实际上根据期望的计算公式

你可以得到 dp[i+1] - dp[i] 这个差值的通式 X

那么就可以递推 dp[i+1] = dp[i] + X 

最后答案就是 dp[n]

下面来讲一下其中期望长度怎么递推、即 len[i]

上面说到要计算 dp[i + 1] = ( len[i] + 1 )^3 * p[i+1]

网上很多题解说  E(x^2) != E(x)^2 或 E(x^3) != E(x)^3

指的就是在 ( len[i] + 1 )^3 的计算这里

将这条公式展开有

len[i]^3 + 3*len[i]^2 + 3*len[i] + 1

那么你不能只算出 len[i] 

然后计算 dp[i+1] = ( len[i]^3 + 3*len[i]^2 + 3*len[i] + 1 ) * p[i+1]

你需要另外递推 len[i]^3  和 len[i]^2 的期望

即先算出 len[i] 记为 a

再用公式求 len[i]^2 记为 b

再用公式求 len[i]^3 记为 c

则 dp[i+1] = (c + 3*b + 3*a + 1)*p[i+1]

#include<bits/stdc++.h>
#define LL long long
#define ULL unsigned long long

#define scl(i) scanf("%lld", &i)
#define scll(i, j) scanf("%lld %lld", &i, &j)
#define sclll(i, j, k) scanf("%lld %lld %lld", &i, &j, &k)
#define scllll(i, j, k, l) scanf("%lld %lld %lld %lld", &i, &j, &k, &l)

#define scs(i) scanf("%s", i)
#define sci(i) scanf("%d", &i)
#define scd(i) scanf("%lf", &i)
#define scIl(i) scanf("%I64d", &i)
#define scii(i, j) scanf("%d %d", &i, &j)
#define scdd(i, j) scanf("%lf %lf", &i, &j)
#define scIll(i, j) scanf("%I64d %I64d", &i, &j)
#define sciii(i, j, k) scanf("%d %d %d", &i, &j, &k)
#define scddd(i, j, k) scanf("%lf %lf %lf", &i, &j, &k)
#define scIlll(i, j, k) scanf("%I64d %I64d %I64d", &i, &j, &k)
#define sciiii(i, j, k, l) scanf("%d %d %d %d", &i, &j, &k, &l)
#define scdddd(i, j, k, l) scanf("%lf %lf %lf %lf", &i, &j, &k, &l)
#define scIllll(i, j, k, l) scanf("%I64d %I64d %I64d %I64d", &i, &j, &k, &l)

#define lson l, m, rt<<1
#define rson m+1, r, rt<<1|1
#define lowbit(i) (i & (-i))
#define mem(i, j) memset(i, j, sizeof(i))

#define fir first
#define sec second
#define VI vector<int>
#define ins(i) insert(i)
#define pb(i) push_back(i)
#define pii pair<int, int>
#define VL vector<long long>
#define mk(i, j) make_pair(i, j)
#define all(i) i.begin(), i.end()
#define pll pair<long long, long long>

#define _TIME 0
#define _INPUT 0
#define _OUTPUT 0
clock_t START, END;
void __stTIME();
void __enTIME();
void __IOPUT();
using namespace std;

int main(void){__stTIME();__IOPUT();


    int n;

    sci(n);

    double len1 = 0, len2 = 0, len3 = 0, dp = 0, ans = 0;

    for(int i=1; i<=n; i++){

        double p; scd(p);

        dp = len3 * p;

        ans += (len3 + 3*len2 + 3*len1 + 1) * p - dp;

        len3 = (len3 + 3*len2 + 3*len1 + 1) * p;
        len2 = (len2 + 2*len1 + 1) * p;
        len1 = (len1 + 1) * p;
    }

    printf("%.1f\n", len3);


















__enTIME();return 0;}


void __stTIME()
{
    #if _TIME
        START = clock();
    #endif
}

void __enTIME()
{
    #if _TIME
        END = clock();
        cerr<<"execute time = "<<(double)(END-START)/CLOCKS_PER_SEC<<endl;
    #endif
}

void __IOPUT()
{
    #if _INPUT
        freopen("in.txt", "r", stdin);
    #endif
    #if _OUTPUT
        freopen("out.txt", "w", stdout);
    #endif
}
View Code

然后来讲一下网上大部分题解的写法递推 dp[i+1] = dp[i] + X

这个 X 怎么求

期望的得分 E( (len[i]+1)^3 ) = E( len[i]^3 ) + X

那么就可以根据这个求出 X = E( (len[i]+1)^3 ) - E( len[i]^3 )

X = E( (len[i]+1)^3 -  len[i]^3 )

   = E( 3*len[i]^2 + 3*len[i] + 1 )

   = ( 3*len[i]^2 + 3*len[i] + 1 ) * p[i+1]

故得到递推式子 dp[i+1] = dp[i] + ( 3*len[i]^2 + 3*len[i] + 1 ) * p[i+1]

在 len[i] 的计算上注意 上面说的 E(x^2) != E(x)^2

#include<bits/stdc++.h>
#define LL long long
#define ULL unsigned long long

#define scl(i) scanf("%lld", &i)
#define scll(i, j) scanf("%lld %lld", &i, &j)
#define sclll(i, j, k) scanf("%lld %lld %lld", &i, &j, &k)
#define scllll(i, j, k, l) scanf("%lld %lld %lld %lld", &i, &j, &k, &l)

#define scs(i) scanf("%s", i)
#define sci(i) scanf("%d", &i)
#define scd(i) scanf("%lf", &i)
#define scIl(i) scanf("%I64d", &i)
#define scii(i, j) scanf("%d %d", &i, &j)
#define scdd(i, j) scanf("%lf %lf", &i, &j)
#define scIll(i, j) scanf("%I64d %I64d", &i, &j)
#define sciii(i, j, k) scanf("%d %d %d", &i, &j, &k)
#define scddd(i, j, k) scanf("%lf %lf %lf", &i, &j, &k)
#define scIlll(i, j, k) scanf("%I64d %I64d %I64d", &i, &j, &k)
#define sciiii(i, j, k, l) scanf("%d %d %d %d", &i, &j, &k, &l)
#define scdddd(i, j, k, l) scanf("%lf %lf %lf %lf", &i, &j, &k, &l)
#define scIllll(i, j, k, l) scanf("%I64d %I64d %I64d %I64d", &i, &j, &k, &l)

#define lson l, m, rt<<1
#define rson m+1, r, rt<<1|1
#define lowbit(i) (i & (-i))
#define mem(i, j) memset(i, j, sizeof(i))

#define fir first
#define sec second
#define VI vector<int>
#define ins(i) insert(i)
#define pb(i) push_back(i)
#define pii pair<int, int>
#define VL vector<long long>
#define mk(i, j) make_pair(i, j)
#define all(i) i.begin(), i.end()
#define pll pair<long long, long long>

#define _TIME 0
#define _INPUT 0
#define _OUTPUT 0
clock_t START, END;
void __stTIME();
void __enTIME();
void __IOPUT();
using namespace std;

int main(void){__stTIME();__IOPUT();


    int n;

    sci(n);

    double len1 = 0, len2 = 0, len3 = 0, dp = 0, ans = 0;

    for(int i=1; i<=n; i++){

        double p; scd(p);

        dp = dp + (3*len2 + 3*len1 + 1) * p;

        len3 = (len3 + 3*len2 + 3*len1 + 1) * p;
        len2 = (len2 + 2*len1 + 1) * p;
        len1 = (len1 + 1) * p;
    }

    printf("%.1f\n", dp);


















__enTIME();return 0;}


void __stTIME()
{
    #if _TIME
        START = clock();
    #endif
}

void __enTIME()
{
    #if _TIME
        END = clock();
        cerr<<"execute time = "<<(double)(END-START)/CLOCKS_PER_SEC<<endl;
    #endif
}

void __IOPUT()
{
    #if _INPUT
        freopen("in.txt", "r", stdin);
    #endif
    #if _OUTPUT
        freopen("out.txt", "w", stdout);
    #endif
}
View Code

 

转载于:https://www.cnblogs.com/LiHior/p/9551843.html

ajumn60282
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
BZOJ第一部分
11-03
【BZOJ第一部分】 BZOJ，全称为“北京邮电大学在线评测系统”（Beijing Zhiyuan Online Judge），是中国最早的一批在线编程竞赛平台之一，它为编程爱好者提供了一个练习和检验编程技能的环境。在这个系统中，用户...
BZOJ全代码
02-05
【标题】"BZOJ全代码"所指的是一份包含BZOJ（BestCoder Zhijiang Online Judge）平台所有编程题目解决方案的压缩文件。BZOJ是一个在线编程竞赛平台，它提供了大量的算法题目供参赛者练习和提交代码，以检验编程能力...
bzoj4318 OSU!（期望dp）
Icefox的博客
01-10 414
e1[i]，表示长度为i的串，末尾连续1的个数的期望。 e2[i]，表示长度为i的串，末尾连续1的个数的平方的期望。 e2[i]=(e2[i−1]+2∗e1[i−1]+1)∗pe2[i]=(e2[i-1]+2*e1[i-1]+1)*p 为什么呢？ 假设一共有n种情况，末尾1的个数分别为xi,出现的概率为pi的话 e1[i]=x1∗p1+x2∗p2+...+xn∗pne1[i]=x_1*p_
BZOJ 4318 OSU! 期望DP
世界
10-30 4756
题目大意：给定一个长度为nn的01串，第ii个位置有aia_i的概率为11，最终得分为01串中所有连在一起1的长度的立方和，求得分的期望假如这个01串使确定的，考虑每新增一个位置，如果这个位置是00，则贡献为00，否则贡献为(x+1)3−x3=3x2+3x+1(x+1)^3-x^3=3x^2+3x+1，其中xx为加入之前最长的全1后缀的长度 现在这个问题变成了期望问题，那么我们只需要维护一个xx的
bzoj 4318 OSU! 期望dp
便纵有千种风情
08-15 390
4318: OSU!Descriptionosu 是一款群众喜闻乐见的休闲软件。 我们可以把osu的规则简化与改编成以下的样子: 一共有n次操作，每次操作只有成功与失败之分，成功对应1，失败对应0，n次操作对应为1个长度为n的01串。在这个串中连续的 X个1可以贡献X^3 的分数，这x个1不能被其他连续的1所包含（也就是极长的一串1，具体见样例解释） 现在给出n，以及每个操作的成功率，请
BZOJ 4318 OSU!期望DP
Monica
06-16 812
相信很多人玩儿过这个游戏…我小时候也玩儿过233Descriptionosu 是一款群众喜闻乐见的休闲软件。 我们可以把osu的规则简化与改编成以下的样子: 一共有n次操作，每次操作只有成功与失败之分，成功对应1，失败对应0，n次操作对应为1个长度为n的01串。在这个串中连续的 X个1可以贡献X^3 的分数，这x个1不能被其他连续的1所包含（也就是极长的一串1，具体见样例解释） 现在给出
bzoj4318: OSU!（Dp）
Hanks_o的博客
03-22 349
题目传送门 。 解法： 太强啦。 f[i]表示前i个位置最后的期望得分。 那么只有这个位置填1的时候才能有贡献。 最长全1后缀为x。那么多了一个1就是x+1 用x推出（x+1）^2的期望。 用x，x^2推出（x+1）^3的期望。 %%%Po姐 代码实现： #include&lt;cstdio&gt; #include&lt;cstring&gt; #include&lt;cs...
2018.09.28【BZOJ4318】OSU!（期望DP）
zxyoi_dreamer的博客（不定期诈尸）
09-28 258
传送门 解析： 不错的期望DP入门题。 思路： 首先，要是没有明白期望的线性性。这道题是无法理解的。 一个常识，要维护高次，我们可以将高次展开，利用低次来维护。 所以我们同时维护长度的期望，长度平方的期望，长度立方的期望。 注意这里可能就有人不懂了，长度平方的期望不等于长度的期望的平方。 同样，长度立方的期望不等于长度的期望的立方。 我们用len1i,len2i,len3ilen1_i,len2...
【BZOJ】4318: OSU! 期望DP
weixin_33866037的博客
03-06 68
【题意】有一个长度为n的01序列，每一段极大的连续1的价值是L^3（长度L）。现在给定n个实数表示该位为1的概率，求期望总价值。n<=10^5。 【算法】期望DP 【题解】后缀长度是一个很关键的量，设g[i]表示前i个的期望后缀长度。根据全期望公式，依赖于第i-1位为0或1：（以下所有公式最后省略+(1-ai)*0） $$g[i]=a_i*(g[i-1]+1)$$ 设f[i]表示前i...
bzoj4318 OSU!
Elijahqi
01-10 239
http://www.elijahqi.win/2018/01/10/bzoj4318-osu/ ‎ Description osu 是一款群众喜闻乐见的休闲软件。 我们可以把osu的规则简化与改编成以下的样子: 一共有n次操作，每次操作只有成功与失败之分，成功对应1，失败对应0，n次操作对应为1个长度为n的01串。在这个串中连续的 X个1可以贡献X^3 的分数，这x个1不能被其他连续的1
BZOJ 4318 OSU!（概率DP）
ZSQ
09-05 285
Description osuosuosu是一款群众喜闻乐见的休闲软件。 我们可以把ososos的规则简化与改编成以下的样子: 一共有nnn次操作，每次操作只有成功与失败之分，成功对应111，失败对应000，nnn次操作对应为111个长度为nnn的010101串。在这个串中连续的XXX个111可以贡献X3X3X^3的分数，这XXX个111不能被其他连续的111所包含（也就是极长的一串11...
【BZOJ4318】OSU!【期望DP】
周周的博客
08-17 173
题意： 有 n 个操作，有 Ai 的概率成功，连续成功 x 次 对答案的贡献是 x ^ 3，断掉就要重新记，求期望的答案  ( 期望DP )这道题和 bzoj3450一样，题意也一样，不过要注意的是平方的期望不等于期望的平方！既然E(x^2)!=(E(x))^2。 分析：第i次点击可能成功也可能失败； 第i次点击时的长度，设第i-1次点击时的长度为x; len=(1-pi)*0+pi*(x...
BZOJ3230相似子串测试数据
08-19
标题"BZOJ3230相似子串测试数据"指的是一个与在线判题平台BZOJ（BestCoder ZOJ，中国大学算法竞赛在线评测系统）相关的编程问题，具体问题编号为3230，涉及的主要知识点是"相似子串"的搜索和处理。在编程竞赛中，...
bzoj problem
08-21
【标题】"bzoj problem" 指的是“八中OJ（Online Judge）问题集”，这是一个在线编程竞赛平台的题目合集。八中OJ可能是某所中学的在线编程训练系统，它提供了大量的编程题目供用户练习和挑战，以提升编程能力和算法...
BZOJ题目镜像
09-14
【BZOJ题目镜像】是为了解决BZOJ（Blue Zoo Online Judge）系统因服务器问题时常出现访问不稳定或无法访问的情况而设立的一个备份资源。BZOJ是中国早期的在线编程竞赛平台，提供了大量的算法题目供编程爱好者进行...
C语言 | Leetcode C语言题解之第337题打家劫舍III
最新发布
m0_59237910的博客
08-15 205
C语言 | Leetcode C语言题解之第337题打家劫舍III
543 二叉树的直径
Noric! 的博客
08-13 373
说明：如果有多个转折点，那么路径上会产生分支，该路径不为合法路径；若没有转折点，可将末尾节点的另外一个子节点加入路径，该路径不是最长的；每一条最长路径中一定存在一个转折点 + 若干非转折点（见说明），只需要把每一个点看做转折点，然后开始构建最长路径，从这些路径长度中选取最大值，即为树的直径。那么，如何构建最长的路径呢，很简单，可以通过递归来解决，在递归之中每一步，情况一：该节点为非转折点，最长路径经过其。作为非转折点能构建的最长的路径长度。情况二：该节点为转折点，最长路径经过其。为转折点构建的路径长度。
Java | Leetcode Java题解之第337题打家劫舍III
m0_57195758的博客
08-15 287
Java | Leetcode Java题解之第337题打家劫舍III
C++ | Leetcode C++题解之第337题打家劫舍III
Ddddddd_158的博客
08-15 230
C++ | Leetcode C++题解之第337题打家劫舍III
bzoj 1040 zjoi2008 骑士 树形dp
08-08
题目描述 有一个 $n$ 个点的棋盘，每个点上有一个数字 $a_i$，你需要从 $(1,1)$ 走到 $(n,n)$，每次只能往右或往下走，每个格子只能经过一次，路径上的数字和为 $S$。定义一个点 $(x,y)$ 的权值为 $a_x+a_y$，求所有满足条件的路径中，所有点的权值和的最小值。 输入格式 第一行一个整数 $n$。 接下来 $n$ 行，每行 $n$ 个整数，表示棋盘上每个点的数字。 输出格式 输出一个整数，表示所有满足条件的路径中，所有点的权值和的最小值。 数据范围 $1\leq n\leq 300$ 输入样例 3 1 2 3 4 5 6 7 8 9 输出样例 25 算法1 (树形dp) $O(n^3)$ 我们可以先将所有点的权值求出来，然后将其看作是一个有权值的图，问题就转化为了在这个图中求从 $(1,1)$ 到 $(n,n)$ 的所有路径中，所有点的权值和的最小值。 我们可以使用树形dp来解决这个问题，具体来说，我们可以将这个图看作是一棵树，每个点的父节点是它的前驱或者后继，然后我们从根节点开始，依次向下遍历，对于每个节点，我们可以考虑它的两个儿子，如果它的两个儿子都被遍历过了，那么我们就可以计算出从它的左儿子到它的右儿子的路径中，所有点的权值和的最小值，然后再将这个值加上当前节点的权值，就可以得到从根节点到当前节点的路径中，所有点的权值和的最小值。 时间复杂度 树形dp的时间复杂度是 $O(n^3)$。 C++ 代码 算法2 (动态规划) $O(n^3)$ 我们可以使用动态规划来解决这个问题，具体来说，我们可以定义 $f(i,j,s)$ 表示从 $(1,1)$ 到 $(i,j)$ 的所有路径中，所有点的权值和为 $s$ 的最小值，那么我们就可以得到如下的状态转移方程： $$ f(i,j,s)=\min\{f(i-1,j,s-a_{i,j}),f(i,j-1,s-a_{i,j})\} $$ 其中 $a_{i,j}$ 表示点 $(i,j)$ 的权值。 时间复杂度 动态规划的时间复杂度是 $O(n^3)$。 C++ 代码
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值