BZOJ4513: [Sdoi2016]储能表(数位dp)
时间:2022-06-16
本文章向大家介绍BZOJ4513: [Sdoi2016]储能表(数位dp),主要内容包括其使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
题意
Sol
一点思路都没有,只会暴力,没想到标算是数位dp??Orz
首先答案可以分成两部分来统计
设
[ f_{i,j}= begin{aligned} ioplus j &left( ioplus j >kright) \ 0 &left( ioplus j <=kright) end{aligned} ]
那么我们要求的就是
[sum_{i=0}^{n - 1} sum_{j = 0}^{m - 1} f(i, j) - k * sum_{i = 0}^{n - 1} sum_{j = 0}^{m - 1} [f(i, j)]]
也就是说,我们要统计出满足条件的数的异或和以及满足条件的数的对数
考虑直接在二进制下数位dp,注意这里我们要记三维状态
(f[len][0/1][0/1][0/1])表示此时到第(len)位,是否顶着(n)的上界,是否顶着(m)的上界,是否顶着(k)的下界
然后直接dp就可以了
// luogu-judger-enable-o2
#include<bits/stdc++.h>
#define Pair pair<LL, LL>
#define MP make_pair
#define fi first
#define se second
#define LL long long
#define int long long
using namespace std;
const int MAXN = 233;
inline LL read() {
char c = getchar(); int x = 0, f = 1;
while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
return x * f;
}
LL N, M, K, mod, Lim, vis[MAXN][2][2][2];
Pair f[MAXN][2][2][2];
void add2(LL &x, LL y) {
if(x + y < 0) x = (x + y + mod);
else x = (x + y >= mod ? x + y - mod : x + y);
}
LL add(LL x, LL y) {
if(x + y < 0) return x + y + mod;
return x + y >= mod ? x + y - mod : x + y;
}
LL mul(LL x, LL y) {
return 1ll * x % mod * y % mod;
}
int Get(LL x) {
int len = 0; while(x) x >>= 1, len++; return len;
}
Pair dfs(int now, int f1, int f2, int f3) {
if(now > Lim) return MP(0, 1);
if(vis[now][f1][f2][f3]) return f[now][f1][f2][f3];
vis[now][f1][f2][f3] = 1;
Pair ans = MP(0, 0);
int L1 = (N >> Lim - now) & 1, L2 = (M >> Lim - now) & 1, L3 = (K >> Lim - now) & 1;
//cout << (f1 &&(!L1)) << endl;
for(int i = 0; i <= (f1 ? L1 : 1); i++) {
for(int j = 0; j <= (f2 ? L2 : 1); j++) {
if(f3 && ((i ^ j) < L3)) continue;
Pair nxt = dfs(now + 1, f1 && (i == L1), f2 && (j == L2), f3 && ((i ^ j) == L3));
add2(ans.se, nxt.se);
add2(ans.fi, add(nxt.fi, mul(nxt.se, mul((i ^ j), (1ll << Lim - now)))));
}
}
return f[now][f1][f2][f3] = ans;
}
int solve() {
memset(vis, 0, sizeof(vis));
memset(f, 0, sizeof(f));
Lim = 0;
N = read(); M = read(); K = read(); mod = read(); N--; M--;
Lim = max(Get(N), max(Get(K), Get(M)));
Pair ans = dfs(1, 1, 1, 1);
return add(ans.fi, -mul(K, ans.se));
}
signed main() {
for(int T = read(); T; T--, printf("%lldn", solve()));
return 0;
}
/*
5000
504363800392059286 554192717354508770 21453916680846604 401134357
*/
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- spring框架应用系列二:component-scan自动扫描注册装配
- 手把手教你,嘴对嘴传达------Apache日志管理日志(rotatelogs分割工具、AWStats日志分析)
- 配合JAVA的AJAX使用
- 手把手教你,嘴对嘴传达------Apache网页优化
- jQuery通过Ajax实现请求后台接口数据
- Git常规操作
- 手把手教你,嘴对嘴传达 ----源码编译安装部署LAMP平台(LAMP平台与编译安装详解,Apache,MySQL与PHP源码编译安装,LAMP平台搭建论坛)
- Vue点击切换样式
- ElementUI引入到vue项目开发
- 手把手教你,嘴对嘴传达------Apache(安全优化防盗链、隐藏版本信息)
- spring框架应用系列三:切面编程(带参数)
- 排障集锦:九九八十一难之第六难!(98)Address already in use: AH00072: make_sock: could not bind to address ::80
- Vue页面中引用自定义组件
- Vue如何引用Vant组件
- js表单验证工具包