问题描述
对于一个长度为 nn 的 01 串 S=x1x2x3…xnS=x1x2x3…xn,香农信息熵的定义为 H(S)=H(S)= −Σ1np(xi)log2(p(xi))−Σ1np(xi)log2(p(xi)),其中 p(0),p(1)p(0),p(1) 表示在这个 01 串中 0 和 1 出现的占比。
比如,对于 S=100S=100 来说,信息熵 H(S)=−13log2(13)−23log2(23)−23log2(23)=1.3083H(S)=−31log2(31)−32log2(32)−32log2(32)=1.3083。
对于一个长度为 2333333323333333 的 01 串, 如果其信息熵为 11625907.579811625907.5798, 且 0 出现次数比 1 少, 那么这个 01 串中 0 出现了多少次?
答案提交
这是一道结果填空的题, 你只需要算出结果后提交即可。本题的结果为一 个整数, 在提交答案时只填写这个整数, 填写多余的内容将无法得分。
运行限制
语言 | 最大运行时间 | 最大运行内存 |
---|---|---|
C++ | 1s | 256M |
C | 1s | 256M |
Java | 2s | 256M |
Python3 | 3s | 256M |
PyPy3 | 3s | 256M |
Go | 3s | 256M |
JavaScript | 3s | 256M |
总通过次数: 1812 | 总提交次数: 1955 | 通过率: 92.7%
难度: 简单 标签: 2023, 暴力, 枚举, 省赛
思路:
这道题暴力直接搜肯定没问题,但是个人觉得难点在于对双浮点数的认识
另外我觉得还有人和我一样,忽略了信息熵是减起来的,于是在考场上累加起来了
另外双浮点数有效位数为16位,小数点后15位,所以在和11625907.5798判断的时候,不要直接用==,要减去11625907.5798看差值是不是到了小数点后四位
#include <iostream>
#include <math.h>
#include <stdlib.h>
using namespace std;
int const n = 23333333;
int main()
{
int ans = 0;
double a, b;//分别为01的占比
for (long long i = 0; i < n; i++)
{
double zkb = 0;
a = i * 1.0 / n;
b = (n - i) * 1.0 / n;
zkb = -1*(a * log2(a) * i + b * log2(b) * (n - i));
if (fabs(zkb - 11625907.5798) <= 0.0001)
{
cout << i;
break;
}
}
// 请在此输入您的代码
return 0;
}