Solution -「NOI 2021」「洛谷 P7740」机器人游戏
\(\mathcal{Description}\)
Link.
自己去读题面叭~
\(\mathcal{Solution}\)
首先,参悟【样例解释 #2】。一种暴力的思路即为钦定集合 \(S\) 内的位置都合法,容斥计数。其中对于每条纸带的每个位置,有三种情况(令 _
为“保持不变”,注意没有被机器人经过的位置都有这种修改):
- 同时存在
_
和*
;或者同时存在0
和1
:只能为空,方案数为 \(1\); - 否则,存在(
_
或*
)且存在(0
或1
):只能为空或01
中的一个,方案数为 \(2\); - 否则,方案数为 \(3\)。
暴力做的复杂度是 \(\mathcal O(2^nn^am^b)\),\(a,b\) 看细节实现,不太可观。
尝试钦定最后一个在 \(S\) 中的位置为 \(r\),那么所有修改长度大于 \(n-r+1\) 的机器人所在的纸带都必须全空。在此限制下带着容斥系数状压 DP,令 \(f(i,S)\) 表示考虑了纸带的位置 \(i\),前若干个位置的选择情况为 \(S\),注意此处 \(|S|\le2^{\min\{i,n-r\}}=\mathcal O(2^{n/2})\),转移时可以 \(\mathcal O(nm)\) 暴力计算当前位置的贡献,得到了 \(\mathcal O(2^{n/2}nm)\) 的算法。
顺带,注意到统计当前位置贡献时可以描述为形如“有多少个机器人的修改序列中,下标属于 \(S\) 的位置上存在 _
/*
/0
/1
”,可以用 std::bitset
优化。最终复杂度为 \(\mathcal O(\frac{2^{n/2}nm}{\omega})\)。
\(\mathcal{Code}\)
/*~Rainybunny~*/
#include <bitset>
#include <cstdio>
#include <cassert>
#include <cstring>
#define rep( i, l, r ) for ( int i = l, rep##i = r; i <= rep##i; ++i )
#define per( i, r, l ) for ( int i = r, per##i = l; i >= per##i; --i )
typedef unsigned long long ULL;
const int MAXN = 32, MAXM = 1e3, MOD = 1e9 + 7;
int n, m, pw2[MAXM + 5], pw3[MAXM + 5];
inline int imin( const int a, const int b ) { return a < b ? a : b; }
inline void subeq( int& a, const int b ) { ( a -= b ) < 0 && ( a += MOD ); }
inline int mul( const int a, const int b ) { return int( 1ll * a * b % MOD ); }
inline int add( int a, const int b ) { return ( a += b ) < MOD ? a : a - MOD; }
inline void addeq( int& a, const int b ) { ( a += b ) >= MOD && ( a -= MOD ); }
struct Robot {
int len, val[MAXN + 5];
inline int operator [] ( const int k ) const {
return k <= len ? val[k] : 0;
}
inline void read() {
static char str[105]; scanf( "%s", str );
val[len = 1] = 0;
for ( int i = 0; str[i]; ++i ) {
if ( str[i] == 'R' ) val[++len] = 0;
else if ( str[i] == '0' ) val[len] = 1;
else if ( str[i] == '1' ) val[len] = 2;
else val[len] = 3 - val[len];
}
}
} rbt[MAXM + 5];
inline int solve( const int r ) {
static int f[2][1 << 17][2];
static std::bitset<MAXM + 5> sig[MAXN + 5][4], g[1 << 17][4];
int allc = 0;
rep ( i, 0, n - 1 ) rep ( j, 0, 3 ) sig[i][j].reset();
rep ( i, 1, m ) if ( rbt[i].len <= n - r + 1 ) {
++allc;
rep ( j, 0, n - 1 )
sig[j][rbt[i][j + 1]].set( i );
}
rep ( i, 0, imin( n - r, r ) ) rep ( j, 0, 3 ) {
g[1 << i][j] = sig[i][j];
}
rep ( S, 1, ( 1 << imin( n - r, r ) << 1 ) - 1 ) if ( S & ( S - 1 ) ) {
rep ( i, 0, 3 ) {
g[S][i] = g[S ^ ( S & -S )][i] | g[S & -S][i];
}
}
auto calc = [&]( const bool sty, const ULL S )->int {
static std::bitset<MAXM + 5> x, y;
x = ( sty ? g[S][3] : g[S][0] & g[S][3] ) | ( g[S][1] & g[S][2] );
y = ~x & ( sty ? g[S][1] | g[S][2] :
( g[S][0] | g[S][3] ) & ( g[S][1] | g[S][2] ) );
int t = y.count();
return mul( pw2[t], pw3[allc - t - x.count()] );
};
memset( f[0], 0, sizeof f[0] ), f[0][0][0] = 1;
for ( int i = 0, sta = 0; i < r; ++i, sta ^= 1 ) {
memset( f[!sta], 0, sizeof f[!sta] );
rep ( S, 0, ( 1 << imin( i, n - r ) ) - 1 ) {
ULL T = ULL( S ) << 1ull & ( ( 1ull << n >> r ) - 1 );
bool flg = ULL( S ) << 1ull << r >> n & 1ull;
rep ( t, 0, 1 ) {
int cur = f[sta][S][t]; f[sta][S][t] = 0;
if ( !cur ) continue;
subeq( f[!sta][r == n ? 0 : T | 1][t || flg || r == n],
mul( cur, calc( t || i + 1 < r, S << 1 | 1 ) ) );
if ( i + 1 < r ) {
addeq( f[!sta][T][t || flg],
mul( cur, calc( 1, S << 1 ) ) );
}
}
}
}
int ret = 0;
rep ( t, 0, 1 ) {
for ( int d = 1, sta = r & 1; d <= n - r; ++d, sta ^= 1 ) {
rep ( i, 0, imin( r, n - r ) - 1 ) {
rep ( j, t, 3 ) g[1 << i][j] = sig[i + d][j];
}
rep ( S, 0, ( 1 << imin( r, n - r ) ) - 1 ) {
if ( S & ( S - 1 ) ) {
rep ( j, t, 3 ) {
g[S][j] = g[S ^ ( S & -S )][j] | g[S & -S][j];
}
}
f[!sta][S][t] = mul( f[sta][S][t], calc( t, S ) );
f[sta][S][t] = 0;
}
}
rep ( S, 0, ( 1 << imin( r, n - r ) ) - 1 ) {
addeq( ret, f[n & 1][S][t] );
}
}
return ret;
}
int main() {
// freopen( "robot.in", "r", stdin );
// freopen( "robot.out", "w", stdout );
scanf( "%d %d", &n, &m );
pw2[0] = pw3[0] = 1;
rep ( i, 1, m ) {
rbt[i].read();
pw2[i] = mul( 2, pw2[i - 1] ), pw3[i] = mul( 3, pw3[i - 1] );
}
int ans = 0;
rep ( i, 1, n ) subeq( ans, solve( i ) );
printf( "%d\n", ans );
return 0;
}
Solution -「NOI 2021」「洛谷 P7740」机器人游戏的更多相关文章
- 「区间DP」「洛谷P1043」数字游戏
「洛谷P1043」数字游戏 日后再写 代码 /*#!/bin/sh dir=$GEDIT_CURRENT_DOCUMENT_DIR name=$GEDIT_CURRENT_DOCUMENT_NAME ...
- 洛谷P1118 数字三角形游戏
洛谷1118 数字三角形游戏 题目描述 有这么一个游戏: 写出一个1-N的排列a[i],然后每次将相邻两个数相加,构成新的序列,再对新序列进行这样的操作,显然每次构成的序列都比上一次的序列长度少1,直 ...
- 洛谷P1274-魔术数字游戏
Problem 洛谷P1274-魔术数字游戏 Accept: 118 Submit: 243Time Limit: 1000 mSec Memory Limit : 128MB Probl ...
- 洛谷 P4012 深海机器人问题【费用流】
题目链接:https://www.luogu.org/problemnew/show/P4012 洛谷 P4012 深海机器人问题 输入输出样例 输入样例#1: 1 1 2 2 1 2 3 4 5 6 ...
- 洛谷P1288 取数游戏II(博弈)
洛谷P1288 取数游戏II 先手必胜的条件需要满足如下中至少 \(1\) 条: 从初始位置向左走到第一个 \(0\) 的位置,经过边的数目为偶数(包含 \(0\) 这条边). 从初始位置向右走到第一 ...
- Solution -「JSOI 2019」「洛谷 P5334」节日庆典
\(\mathscr{Description}\) Link. 给定字符串 \(S\),求 \(S\) 的每个前缀的最小表示法起始下标(若有多个,取最小的). \(|S|\le3\time ...
- Solution -「洛谷 P4372」Out of Sorts P
\(\mathcal{Description}\) OurOJ & 洛谷 P4372(几乎一致) 设计一个排序算法,设现在对 \(\{a_n\}\) 中 \([l,r]\) 内的元素排 ...
- Solution -「POI 2010」「洛谷 P3511」MOS-Bridges
\(\mathcal{Description}\) Link.(洛谷上这翻译真的一言难尽呐. 给定一个 \(n\) 个点 \(m\) 条边的无向图,一条边 \((u,v,a,b)\) 表示从 ...
- Solution -「APIO 2016」「洛谷 P3643」划艇
\(\mathcal{Description}\) Link & 双倍经验. 给定 \(n\) 个区间 \([a_i,b_i)\)(注意原题是闭区间,这里只为方便后文描述),求 \(\ ...
随机推荐
- [Beyond Compare] 排除/忽略 .svn 文件夹
[Beyond Compare] Exclude .svn folders Beyond Compare 3 Session >> Session Settings... >> ...
- vue爬坑之路(插件安装)
npm install vue-table-with-tree-grid --save import ZkTable from 'vue-table-with-tree-grid' Vue.use(Z ...
- .net core 和 WPF 开发升讯威在线客服系统:调用百度翻译接口实现实时自动翻译
业余时间用 .net core 写了一个在线客服系统.并在博客园写了一个系列的文章,写介绍这个开发过程. 我把这款业余时间写的小系统丢在网上,陆续有人找我要私有化版本,我都给了,毕竟软件业的初衷就是免 ...
- Bootstrap实战 - 注册和登录
一.介绍 注册和登录在社交和商业网站中是必不可少的一个部分. 二.知识点 2.1 标签页 2.1.1 基础标签页 标签页的使用与导航栏类似,同时都依赖于基础样式 nav,不同的是附加样式变成了 nav ...
- LCT小记
不用说了,直接上怎么 die( 千万不要和 Treap 一样写左旋 zig 和右旋 zag,莫名死亡.Splay 只支持一个 rotate 上旋一个节点即可. splay() 之前记得弄一个栈存储 u ...
- 【Android】安卓中的存储
[Android]安卓中的存储 1.存储在App内部 最简单的一种.在尝试过程中发现,手机中很多文件夹都没有权限读写.我们可以将我们需要写的文件存放到App中的files文件夹中,当然我们有权限在整个 ...
- spring 事务失效的几种场景
以下场景是基于mysql数据库,InnoDB的存储引擎. 一.没有添加@Transactional注解 二.方法声明是private或者static 三.没有抛出异常而是try catch了异常 下面 ...
- unity3d之sokect通信
using System.Collections; using System.Collections.Generic; using UnityEngine; using System; using S ...
- 从零开始, 开发一个 Web Office 套件 (1): 富文本编辑器
这是一个系列博客, 最终目的是要做一个基于HTML Canvas 的, 类似于微软 Office 的 Web Office 套件, 包括: 文档, 表格, 幻灯片... 等等. 富文本编辑器 万里长征 ...
- JUC之阻塞队列(BlockingQueue)基础
阻塞队列 阻塞队列(BlockingQueue)是一个支持两个附加操作的队列.这两个附加的操作支持阻塞的插入和移除方法. 阻塞队列常用于生产者和消费者的场景: 当队列为空的,从队列中获取元素的操作将会 ...