首页 > 其他分享 >AGC041F Histogram Rooks

AGC041F Histogram Rooks

时间:2024-04-26 18:23:39浏览次数:20  
标签:Rooks int AGC041F void Histogram MAXN using dp MOD

考察一个合法的棋盘,有一些列上有车,其他的列没有车。

有车的列所有格子都被覆盖,没有车的列上的每个格子都需要被同一行的车覆盖。

即考察所有极长的行连续段,如果这个连续段涉及的列里包含没有车的列,那么这个行连续段里必须有车。

考虑枚举有车的列的集合 \(S\),对于每个行连续段,记它涉及的列集合与 \(S\) 的交集大小为 \(c\),如果它包含没有车的列则方案数为 \(2^c - 1\),否则方案数为 \(2^c\)。

容易发现这么计数不能保证 \(S\) 里的列都有车,再枚举 \(T \subseteq S\) 表示钦定 \(T\) 里的列没有车,容斥系数为 \((-1)^T\),把上一段里与 \(S\) 的交集改为与 \(S - T\) 的交集即可。

一个自然的分治过程是找到当前最低的列 \(i\) 然后把最底下的 \(h_i\) 行删掉,这样棋盘就被分成了不相关的两个部分。根据这个分治过程进行 dp 相当于建出笛卡尔树然后自底向上进行 dp。

设 \(dp_{u, i, 0/1}\) 表示考察 \(u\) 的子树,共有 \(i\) 个列属于 \(S - T\),是否有不属于 \(S\)(即没有车)的列。

合并两个儿子是树形背包,然后枚举当前的最低列属于三个部分的哪一个进行转移,再对每个在分治过程中删掉的行连续段乘上 \(2^c\) 或 \(2^c - 1\) 的系数即可。

注意乘上系数的快速幂复杂度可以摊到每个格子上,因此时间复杂度 \(O(n^2)\)。

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using ull = unsigned long long;
using ld = long double;
using i128 = __int128;
using u128 = unsigned __int128;
template <typename T>
void chkmax(T &x, const T &y) {
  if (x < y) x = y;
}
template <typename T>
void chkmin(T &x, const T &y) {
  if (y < x) x = y;
}
constexpr int MOD = 998244353;
void inc(int &x, int y) {
  x += y;
  if (x >= MOD) x -= MOD;
}
int add(int x, int y) {
  x += y;
  return x >= MOD ? x - MOD : x;
}
int neg(int x) { return x ? MOD - x : 0; }
void upd(int &x, ll y) { x = (x + y) % MOD; }
int qpow(int x, int y) {
  int ret = 1;
  for (; y; y >>= 1, x = (ll)x * x % MOD)
    if (y & 1) ret = (ll)ret * x % MOD;
  return ret;
}
constexpr int MAXN = 405;
int n, h[MAXN], stk[MAXN], tp, ls[MAXN], rs[MAXN];
int dp[MAXN][MAXN][2], pw2[MAXN];
void dfs(int u, int pre, int l, int r) {
  if (!u) return;
  dfs(ls[u], u, l, u - 1);
  dfs(rs[u], u, u + 1, r);
  for (int i = 0; i <= u - l; ++i) {
    for (int x = 0; x < 2; ++x) {
      if (!dp[ls[u]][i][x]) continue;
      for (int j = 0; j <= r - u; ++j) {
        for (int y = 0; y < 2; ++y) {
          if (!dp[rs[u]][j][y]) continue;
          upd(dp[u][i + j][x || y], (ll)dp[ls[u]][i][x] * dp[rs[u]][j][y]);
        }
      }
    }
  }
  for (int i = r - l; i >= 0; --i) {
    inc(dp[u][i + 1][0], dp[u][i][0]);
    inc(dp[u][i + 1][1], dp[u][i][1]);
    dp[u][i][1] = dp[u][i][0];
    dp[u][i][0] = neg(dp[u][i][0]);
  }
  for (int i = 0; i <= r - l + 1; ++i) {
    dp[u][i][0] = (ll)dp[u][i][0] * qpow(pw2[i], h[u] - h[pre]) % MOD;
    dp[u][i][1] = (ll)dp[u][i][1] * qpow(pw2[i] - 1, h[u] - h[pre]) % MOD;
  }
}
int main() {
  ios::sync_with_stdio(false);
  cin.tie(nullptr);
  cin >> n;
  for (int i = 1; i <= n; ++i) cin >> h[i];
  for (int i = 1; i <= n; ++i) {
    while (tp && h[i] < h[stk[tp]]) ls[i] = stk[tp--];
    if (tp) rs[stk[tp]] = i;
    stk[++tp] = i;
  }
  for (int i = 1; i < tp; ++i) rs[stk[i]] = stk[i + 1];
  pw2[0] = 1;
  for (int i = 1; i <= n; ++i) pw2[i] = add(pw2[i - 1], pw2[i - 1]);
  dp[0][0][0] = 1;
  dfs(stk[1], 0, 1, n);
  int ans = 0;
  for (int i = 0; i <= n; ++i)
    for (int j = 0; j < 2; ++j) inc(ans, dp[stk[1]][i][j]);
  cout << ans << "\n";
  return 0;
}
/*
g++ F.cpp -o F -std=c++14 -O2 -Wall -Wextra -Wshadow -g
-fsanitize=address,undefined
*/

标签:Rooks,int,AGC041F,void,Histogram,MAXN,using,dp,MOD
From: https://www.cnblogs.com/JCY-std/p/18160639

相关文章

  • CF1411C - Peaceful Rooks | 思维
    links在一个\(n\timesn\)的棋盘上有\(m(m<n)\)个棋子。若棋子处于同一行或同一列便认为他们可以互相攻击。初始时棋子之间均不可互相攻击。你可以进行若干次操作,每次操作可以将棋子纵向移动任意格或横向移动任意格,要求移动之后棋子之间不能互相攻击。求使得棋子均处在主......
  • Atcoder Grand Contest 041 F - Histogram Rooks
    考虑容斥。我们钦定一些格子组成的集合不能被覆盖,设为\(A\)。把与\(A\)中的点同行同列的点抠掉,剩余的点则是可放可不放的,总方案数就是\(2^{\text{剩余点的个数}}\),乘以\((-1)^{|A|}\)并求和即可。这个做法直接优化显然不行。我们考虑设\(A\)中的点所在的列组成的不可重集......
  • [AGC041F] Histogram Rooks 题解
    题目链接点击打开链接题目解法好牛(难)的题!!!所有都被覆盖不难想到容斥暴力的容斥是钦定集合\(S\)中的位置都没被覆盖,然后把不能填棋子的点都去掉,假设剩下的不受限制的点的个数为\(c\),则答案为\(\sum2^c(-1)^{|S|}\)这个暴力是很难直接优化的如果一列被有点被钦定了,那么......
  • Prometheus最佳实践 Summary和Histogram
    本文分享自华为云社区《Prometheus最佳实践Summary和Histogram》,作者:张俭。前言Histogram和Summary都是复杂的指标,不仅仅是因为直方图和summary包含了多个时间序列,而且它们还较难使用正确。观测中的Count和SumHisto和summary都是采样观测,典型的采样维度有 响应大小 和 ......
  • 无涯教程-Seaborn - 直方图(Histogram)
    直方图表示数据分布,方法是沿数据范围形成条形图,然后绘制条形图以显示落入每个条形图的观察数。Seaborn附带了一些数据集,在前几章中只使用了很少的数据集。无涯教程已经学习了如何加载数据集以及如何查找可用数据集列表。importpandasaspdimportseabornassbfrommatplot......
  • 《计算机科学概论 第12版》[美]J.Glenn Brookshear Denis Brylow 译者:刘艺,吴英,毛
    《计算机科学概论》是计算机科学概论课程的经典教材,全书对计算机科学做了百科全书式的精彩阐述,充分展现了计算机科学的历史背景、发展历程和新的技术趋势。《计算机科学概论》首先介绍的是信息编码及计算机体系结构的基本原理,进而讲述操作系统和组网及因特网,接着探讨算法、程序设......
  • 计算机科学概论 (第10版) 作者: [美] J.Glenn Brookshear 译者: 刘艺 / 肖成海 / 马小
    计算机科学概论(第10版)  更新图书信息或封面作者: [美]J.GlennBrookshear出版社: 人民邮电出版社出品方: 图灵教育原作名: ComputerScience:AnOverview译者: 刘艺 / 肖成海 / 马小会出版年: 2009-9页数: 411定价: 59.00元装帧: 平装丛书: 图灵......
  • Grafana学习(5)——Introduction to histograms and heatmaps
    Ahistogramisagraphicalrepresentationofthedistributionofnumericaldata.Itgroupsvaluesintobuckets(sometimesalsocalledbins)andthencountshowmanyvaluesfallintoeachbucket.Insteadofgraphingtheactualvalues,histogramsgraphthe......
  • SPOJ1805 HISTOGRA - Largest Rectangle in a Histogram 题解
    LinkSPOJ1805HISTOGRA-LargestRectangleinaHistogramQuestion在一条水平线上有\(n\)个高为\(a_i\)的矩形,求包含于这些矩形的最大子矩形面积。Solution我们定义\(L_i\)表示有\(a_i\)这个高度的一根悬线,往左最多能平移到什么位置初始化显然,\(a_i=i\)考虑转移......
  • AtCoder Grand Contest 041 F Histogram Rooks
    洛谷传送门AtCoder传送门神题!!!!!!!!!/bx全部格子都被覆盖不好处理,考虑钦定\(k\)个格子不被覆盖,容斥系数就是\((-1)^k\)。发现网格的行不一定连续,但是列是连续的。如果一列有格子被钦定,那么这一列就不能放棋子。由此想到枚举不能放棋子的列(至少有一个棋子被钦定的列)集合\(S\),把......