Article Index

趣题 Entropy 做题笔记

Published on

题意

题解

首先可以以一个对于当前的 $H$ 合适的长度,随机合理的次数,得到一个较接近 $H$ 的解。然后考虑如何微调。

如果加入一个出现次数最多的数,那么熵就会减小。反之则会增大。

而这个增减幅度会随着字符串的长度越来越小,于是最终可以接近 $H$,而 $1000$ 的长度非常充足。

对于一些难处理的数,可以尝试打表。

代码

// :/

double calc(string s) {
  double h = 0;
  vector<tp> cnt(64, 0);
  for (auto& i : s) ++cnt[i];
  for (auto& i : cnt) {
    if (i == 0) continue;
    double p = i * 1. / s.size();
    h -= p * log2(p);
  }
  return h;
}

string random_str(tp n) {
  string s;
  while (n--) s.push_back(lib::rnd(0, 63));
  return s;
}

char mapping[] = { '0', '1', '2', '3', '4', '5', '6', '7', '8', '9', '.', ' ',
  'a', 'b', 'c', 'd', 'e', 'f', 'g', 'h', 'i', 'j', 'k', 'l', 'm',
  'n', 'o', 'p', 'q', 'r', 's', 't', 'u', 'v', 'w', 'x', 'y', 'z',
  'A', 'B', 'C', 'D', 'E', 'F', 'G', 'H', 'I', 'J', 'K', 'L', 'M',
  'N', 'O', 'P', 'Q', 'R', 'S', 'T', 'U', 'V', 'W', 'X', 'Y', 'Z' };

void STRUGGLING([[maybe_unused]] unsigned TEST_NUMBER) {
  double h; bin >> h;
  if (abs(h - 0.01) < 0.0001) { bin << "nnnZnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnn"; return; }
  string s = random_str(h * 99 + 3);
  tp round = 67777 - h * 11000;
  double lst = abs(calc(s) - h);
  while (round--) {
    string t = random_str(h * 99 + 3);
    if (abs(calc(t) - h) < lst) { s = t; lst = abs(calc(s) - h); }
  }
  double tmp = calc(s);
  while (abs(tmp - h) >= 0.0049) {
    vector<tp> cnt(64, 0);
    for (auto& i : s) ++cnt[i];
    tp goal;
    if (tmp > h) goal = *max_element(FULL(cnt));
    else {
      goal = s.size();
      for (auto& i : s) {
        if (cnt[i] == 0) continue;
        if (cnt[i] < goal) goal = cnt[i];
      }
    }
    for (auto& i : s) {
      if (cnt[i] == goal) { s.push_back(i); break; }
    }
    tmp = calc(s);
  }
  for (auto& i : s) bin << mapping[i];
  bin << '\n';
  dbg(s.size());
}

void MIST() {
}

// :\ */

No comments

Post a comment