http://codeforces.com/problemset/problem/1073/E

题意 给定K,L,R,求L~R之间最多不包含超过K个数码的数的和。

显然这是一道数位dp,在做的过程中会发现为了统计数码是否出现过这个状态需要用到状态压缩

因为不同位置出现的数贡献不同,除了传统的dp数组之外还需要记录一个tot来统计这个位置之后出现数字的个数方便后面计算答案。

仔细对比了一下为什么我要开五维数组而网上的题解只需要开二维的数组,发现一是因为网上的题解不对前导0和limit的情况进行记忆化搜索,仔细想了一下这两种情况确实分支的情况比较小,为了他开两倍的空间的确不那么值得。

二是因为dp除了记录位置和是否出现过这个state之外,每个位置是什么数事实上是不需要记录的,对于一个位置和一个状态事实上后续就有固定的答案,不需要额外开一个数字来记录,对于后效答案的更新可以放在进入记忆话搜索之前,因而数字这位数字也可以省略掉。

#include <map>
#include <set>
#include <ctime>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <sstream>
#include <iostream>
#include <algorithm>
#include <functional>
using namespace std;
inline int read(){int now=;register char c=getchar();for(;!isdigit(c);c=getchar());
for(;isdigit(c);now=now*+c-'',c=getchar());return now;}
#define For(i, x, y) for(int i=x;i<=y;i++)
#define _For(i, x, y) for(int i=x;i>=y;i--)
#define Mem(f, x) memset(f,x,sizeof(f))
#define Sca(x) scanf("%d", &x)
#define Sca2(x,y) scanf("%d%d",&x,&y)
#define Sca3(x,y,z) scanf("%d%d%d",&x,&y,&z)
#define Scl(x) scanf("%lld",&x);
#define Pri(x) printf("%d\n", x)
#define Prl(x) printf("%lld\n",x);
#define CLR(u) for(int i=0;i<=N;i++)u[i].clear();
#define LL long long
#define ULL unsigned long long
#define mp make_pair
#define PII pair<int,int>
#define PIL pair<int,long long>
#define PLL pair<long long,long long>
#define pb push_back
#define fi first
#define se second
typedef vector<int> VI;
const double eps = 1e-;
const int maxn = ;
const int INF = 0x3f3f3f3f;
const LL mod = ;
LL L,R;
int K;
LL ten[maxn];
int str[maxn];
LL dp[][][][][];
LL tot[][][][][];
int cnt;
LL dfs(int pos,int num,int use,int k,int limit,bool zero){
if(pos == ){
tot[pos][num][use][limit][zero] = ;
return num;
}
if(~dp[pos][num][use][limit][zero]) return dp[pos][num][use][limit][zero];
LL ed = limit?str[pos - ]:;
LL sum = ;
for(int i = ; i <= ed; i ++){
if(!(use & ( << i)) && (zero || i)){
if(k + > K) continue;
LL x = dfs(pos - ,i,use | ( << i),k + ,limit && i == str[pos - ],zero || i);
sum = (sum + x) % mod;
tot[pos][num][use][limit][zero] += tot[pos - ][i][use | ( << i)][limit && i == (str[pos - ])][zero || i];
}else{
LL x = dfs(pos - ,i,use,k,limit && i == str[pos - ],zero || i);
sum = (sum + x) % mod;
tot[pos][num][use][limit][zero] += tot[pos - ][i][use][limit && i == (str[pos - ])][zero || i];
}
}
tot[pos][num][use][limit][zero] %= mod;
if(zero) sum = (sum + tot[pos][num][use][limit][zero] * num % mod * ten[pos]) % mod;
dp[pos][num][use][limit][zero] = sum;
return sum;
}
LL solve(LL x){
Mem(dp,-); Mem(tot,);
if(x <= ) return ;
cnt = ;
while(x){
str[cnt++] = x % ;
x /= ;
}
return dfs(cnt,,,,,);
}
int main()
{
ten[] = ;
for(int i = ; i <= ; i ++) ten[i] = (ten[i - ] * ) % mod;
scanf("%lld%lld%d",&L,&R,&K);
LL ans1 = solve(R),ans2 = solve(L - );
//cout << ans1 << " " << ans2 << endl;
Prl((ans1 - ans2 + mod) % mod);
#ifdef VSCode
system("pause");
#endif
return ;
}

CodeForces1073E 数位dp+状压dp的更多相关文章

  1. 【BZOJ】1076 [SCOI2008]奖励关 期望DP+状压DP

    [题意]n种宝物,k关游戏,每关游戏给出一种宝物,可捡可不捡.每种宝物有一个价值(有负数).每个宝物有前提宝物列表,必须在前面的关卡取得列表宝物才能捡起这个宝物,求期望收益.k<=100,n&l ...

  2. CCF 201312-4 有趣的数 (数位DP, 状压DP, 组合数学+暴力枚举, 推公式, 矩阵快速幂)

    问题描述 我们把一个数称为有趣的,当且仅当: 1. 它的数字只包含0, 1, 2, 3,且这四个数字都出现过至少一次. 2. 所有的0都出现在所有的1之前,而所有的2都出现在所有的3之前. 3. 最高 ...

  3. hdu 4352 "XHXJ's LIS"(数位DP+状压DP+LIS)

    传送门 参考博文: [1]:http://www.voidcn.com/article/p-ehojgauy-ot.html 题解: 将数字num字符串化: 求[L,R]区间最长上升子序列长度为 K ...

  4. [转]状态压缩dp(状压dp)

    状态压缩动态规划(简称状压dp)是另一类非常典型的动态规划,通常使用在NP问题的小规模求解中,虽然是指数级别的复杂度,但速度比搜索快,其思想非常值得借鉴. 为了更好的理解状压dp,首先介绍位运算相关的 ...

  5. 状态压缩dp 状压dp 详解

    说到状压dp,一般和二进制少不了关系(还常和博弈论结合起来考,这个坑我挖了还没填qwq),二进制是个好东西啊,所以二进制的各种运算是前置知识,不了解的话走下面链接进百度百科 https://baike ...

  6. 洛谷 P3343 - [ZJOI2015]地震后的幻想乡(朴素状压 DP/状压 DP+微积分)

    题面传送门 鸽子 tzc 竟然来补题解了,奇迹奇迹( 神仙题 %%%%%%%%%%%% 解法 1: 首先一件很明显的事情是这个最小值可以通过类似 Kruskal 求最小生成树的方法求得.我们将所有边按 ...

  7. 51nod 1673 树有几多愁(链表维护树形DP+状压DP)

    题意 lyk有一棵树,它想给这棵树重标号. 重标号后,这棵树的所有叶子节点的值为它到根的路径上的编号最小的点的编号. 这棵树的烦恼值为所有叶子节点的值的乘积. lyk想让这棵树的烦恼值最大,你只需输出 ...

  8. BZOJ3836 [Poi2014]Tourism 【树形dp +状压dp】

    题目链接 BZOJ3836 题解 显然这是个\(NP\)完全问题,此题的解决全仗任意两点间不存在节点数超过10的简单路径的性质 这意味着什么呢? \(dfs\)树深度不超过\(10\) \(10\)很 ...

  9. SCUT - 254 - 欧洲爆破 - 概率dp - 状压dp

    https://scut.online/p/254 思路很清晰,写起来很恶心. #include<bits/stdc++.h> using namespace std; #define l ...

  10. bzoj1076: [SCOI2008]奖励关(期望dp+状压dp)

    1076: [SCOI2008]奖励关 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 2989  Solved: 1557[Submit][Statu ...

随机推荐

  1. AtCoder Beginner Contest 122 解题报告

    手速选手成功混进rated only里面的前30名,但是总排名就到110+了... A - Double Helix #include <bits/stdc++.h> #define ll ...

  2. MD5进行加密操作

    package com.dyy.test; import java.security.MessageDigest; import java.security.NoSuchAlgorithmExcept ...

  3. Django auth认证系统

    Django自带的用户认证 我们在开发一个网站的时候,无可避免的需要设计实现网站的用户系统.此时我们需要实现包括用户注册.用户登录.用户认证.注销.修改密码等功能,这还真是个麻烦的事情呢. Djang ...

  4. Python面试题练习

    1.实现1--100之和 #解答一 print sum(xrange(101)) #解答二 s=0 for i in xrange(101): s = s + i print s 2.如何在一个函数内 ...

  5. [HDU2065] "红色病毒"问题

    传送门:>Here< 题意:现在有一长度为N的字符串,满足一下条件: (1) 字符串仅由A,B,C,D四个字母组成; (2) A出现偶数次(也可以不出现); (3) C出现偶数次(也可以不 ...

  6. 使用python和selenium写一个百度搜索的case

    今天练习的内容主要写了一个小功能,在百度上搜索某词汇,然后实现web上的back功能 代码如下: import unittest from selenium import webdriver from ...

  7. X问题 HDU - 1573(excrt入门题)

    X问题 Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submiss ...

  8. pytorch搭建简单网络

    pytorch搭建一个简单神经网络 import torch import torch.nn as nn # 定义数据 # x:输入数据 # y:标签 x = torch.Tensor([[0.2, ...

  9. linux常用目录简介

    简单归纳: /bin目录(binary)是二进制执行文件目录,主要用于具体应用 /sbin目录(system binary)是系统管理员专用的二进制代码存放目录,主要用于系统管理 附转载一份: /bi ...

  10. 使用nmon监控Linxu服务器性能

    nmon是IBM开发的Linux软件工具.能够监控多项Linux服务,最常见的如:CPU使用率.内存使用率.硬盘读写.网卡流量等. 并能设置参数,将记录的数据保存到文件,利用Excel做统计分析. 安 ...