题目描述

小Z在玩一个叫做《淘金者》的游戏。游戏的世界是一个二维坐标。X轴、Y轴坐标范围均为1..N。初始的时候,所有的整数坐标点上均有一块金子,共N*N块。

一阵风吹过,金子的位置发生了一些变化。细心的小Z发现,初始在(i,j)坐标处的金子会变到(f(i),fIj))坐标处。其中f(x)表示x各位数字的乘积,例如f(99)=81,f(12)=2,f(10)=0。

如果金子变化后的坐标不在1..N的范围内,我们认为这块金子已经被移出游戏。同时可以发现,对于变化之后的游戏局面,某些坐标上的金子数量可能不止一块,而另外一些坐标上可能已经没有金子。这次变化之后,游戏将不会再对金子的位置和数量进行改变,玩家可以开始进行采集工作。

小Z很懒,打算只进行K次采集。每次采集可以得到某一个坐标上的所有金子,采集之后,该坐标上的金子数变为0。

现在小Z希望知道,对于变化之后的游戏局面,在采集次数为K的前提下,最多可以采集到多少块金子? 答案可能很大,小Z希望得到对1000000007(10^9+7)取模之后的答案。

输入输出格式

输入格式:

共一行,包含两介正整数N,K。

输出格式:

一个整数,表示最多可以采集到的金子数量。

输入输出样例

输入样例#1: 复制

1 2  5
输出样例#1: 复制

18

说明

N < = 10^12 ,K < = 100000

对于100%的测试数据:K < = N^2

题解

  大概是少数我能看得懂题解的数位dp……

  我们可以先求出每个数能被转移到的次数$f[i]$,那么点$(i,j)$的金子就是$f[i]*f[j]$

  然后这个$f$就用数位dp来求,具体可以看代码

  因为题目中所有数的积的总数也就$10^4$左右,可以直接深搜+离散给弄出来

  然后考虑状态$f[i][j][k]$表示$i$位数,积为$j$(离散后的),$k$表示是否超过$n$

  然后直接递推

  再然后直接用优先队列求答案就行

  然后代码里细节很多,都写在注解里了,自己注意

 //minamoto
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<queue>
#define ll long long
using namespace std;
const int N=4e5+,mod=1e9+;
ll f[][N][],num[N],sz[N];
int tot,K,a[],len;ll n,ans;
struct node{
int x,y;ll v;
node(int x,int y):x(x),y(y){v=sz[x]*sz[y];}
inline bool operator <(const node &b)const
{return v<b.v;}
};
priority_queue<node> q;
inline bool cmp(ll x,ll y){return x>y;}
void dfs(int now,int dep,ll mul){
if(dep==len) num[++tot]=mul;
else{
for(int i=now;i<;++i) dfs(i,dep+,mul*i);
}
}
int main(){
//freopen("testdata.in","r",stdin);
scanf("%lld%d",&n,&K);
ll tmp=n;
while(n) a[++len]=n%,n/=;
num[++tot]=;dfs(,,);
sort(num+,num++tot);
tot=unique(num+,num++tot)-num-;
f[][][]=;//num[2]是1
for(int i=;i<=len;++i) for(int j=;j<=tot;++j) for(int k=;k<;++k) if(f[i][j][k]){
for(int x=i==?:;x<;++x){
//0只能有1个
int r=lower_bound(num+,num+tot+,num[j]*x)-num;
f[i+][r][(k+x)>a[i+]]+=f[i][j][k];
//这里应该写成x+k>a[i+1]来判断是否大于而不是x>k
//首先这里是从低位到高位枚举,所以默认高位全是0
//然后如果x==a[i+1]那就要看后面有没有超来判断当前数是否超了
//所以要+k
}
}
for(int i=;i<=tot;++i){
for(int j=;j<=len-;++j) sz[i]+=f[j][i][]+f[j][i][];
//上面是这样写,因为全都是低位,低位就算超了,高位还是小于限制
sz[i]+=f[len][i][];
//最高位不能超
}
sort(sz+,sz++tot,cmp);
q.push(node(,));
while(!q.empty()&&K){
node u=q.top();q.pop();
ans=(ans+u.v)%mod;
if(!--K) break;
if(u.x!=u.y){
//如果两个数不一样,那么这样的格子有两个
ans=(ans+u.v)%mod;
if(!--K) break;
q.push(node(u.x+,u.y));
}
if(u.x==) q.push(node(u.x,u.y+));
//排序之后肯定从大到小最优
//所以枚举是跑x+1还是y+1
}
printf("%lld\n",ans);
return ;
}

bzoj 3131 [Sdoi2013]淘金(数位dp)的更多相关文章

  1. BZOJ 3131 [SDOI2013]淘金 - 数位DP

    传送门 Solution 这道数位$DP$看的我很懵逼啊... 首先我们肯定要先预处理出 $12$位乘起来的所有的可能情况, 记录入数组 $b$, 发现个数并不多, 仅$1e4$不到. 然后我们考虑算 ...

  2. Bzoj 3131 [Sdoi2013]淘金 题解

    3131: [Sdoi2013]淘金 Time Limit: 30 Sec  Memory Limit: 256 MBSubmit: 733  Solved: 363[Submit][Status][ ...

  3. bzoj 3131 [Sdoi2013]淘金(数位DP+优先队列)

    Description 小Z在玩一个叫做<淘金者>的游戏.游戏的世界是一个二维坐标.X轴.Y轴坐标范围均为1..N.初始的时候,所有的整数坐标点上均有一块金子,共N*N块.    一阵风吹 ...

  4. [SDOI2013]淘金 数位DP

    做了好久.... 大致思路: 求出前k大的方格之和即为答案, 先考虑一维的情况,设f[i]为数位上各个数相乘为i的数的总数,也就是对于数i,有f[i]个数它们各个位相乘为i, 再拓展到二维,根据乘法原 ...

  5. bozoj3131: [Sdoi2013]淘金 数位dp

    链接 https://www.lydsy.com/JudgeOnline/problem.php?id=3131 思路 1. 函数值的素因子只有2.3.5.7 由他们组成的状态不多,爆搜的时候即使搜不 ...

  6. bzoj 3131: [Sdoi2013]淘金

    #include<cstdio> #include<iostream> #include<queue> #include<algorithm> #def ...

  7. BZOJ 3652: 大新闻(数位DP+概率论)

    不得不说数位DP和博弈论根本不熟啊QAQ,首先这道题嘛~~~可以分成两个子问题: 有加密:直接算出0~n中二进制每一位为0或为1分别有多少个,然后分位累加求和就行了= = 无加密:分别算出0~n中二进 ...

  8. BZOJ 1833 数字计数 数位DP

    题目链接 做的第一道数位DP题,听说是最基础的模板题,但还是花了好长时间才写出来..... 想深入了解下数位DP的请点这里 先设dp数组dp[i][j][k]表示数位是i,以j开头的数k出现的次数 有 ...

  9. BZOJ 4521 [CQOI2016]手机号码 - 数位DP

    Description 在$[L, R]$找出有几个数满足两个条件 : 1 : 不同时含有$4$ 和 $8$ 2 : 至少有$3$个相邻的数相同 Solution 非常容易的数位DP, $pos$ 为 ...

随机推荐

  1. java代码继承super

    总结:多态 :. 当重写父类的方法的时,子类对象名可以调用父类的方法,以及不带参的构造方法 package com.addd; public class rr { int a, b; String c ...

  2. Java-API:java.lang百科

    ylbtech-Java-API:java.lang百科 java.lang是提供利用 Java 编程语言进行程序设计的基础类.最重要的类是Object(它是类层次结构的根)和 Class(它的实例表 ...

  3. UE4流关卡

    转自:http://blog.ch-wind.com/ue4-level-streaming/ 流关卡可以使得关卡内容只在玩家“需要”的时候才加载,在很多游戏中都有使用这个技术. 当前UE4版本4.1 ...

  4. Solaris10上如何识别新增加的HDLM LUN

    先在磁盘阵列上将新加LUN映射给主机组,然后在光纤交换机上增加相关zone信息.以下是Solaris10上需要执行的操作步骤. 在Solaris10上重新扫描磁盘 -bash-3.2# cfgadm ...

  5. 阿里云服务器ubuntu安装redis2.8.13

    阿里云服务器ubuntu安装redis2.8.13 2014-09-04 16:14 |  coding云 |  2198次阅读 | 暂无评论   一.下载redis 可以先下载到本地,然后ftp到服 ...

  6. DAY15-web框架本质及第一个Django实例

    Web框架本质 我们可以这样理解:所有的Web应用本质上就是一个socket服务端,而用户的浏览器就是一个socket客户端. 这样我们就可以自己实现Web框架了. 半成品自定义web框架 impor ...

  7. 用原生JS和html5进行视频截图并保存到本地

    <!doctype html> <html> <head> <meta charset="utf-8"> <title> ...

  8. python爬虫实战(3)--图片下载器

    本篇目标 1.输入关键字能够根据关键字爬取百度图片 2.能够将图片保存到本地文件夹 1.URL的格式 进入百度图片搜索apple,这时显示的是瀑布流版本,我们选择传统翻页版本进行爬取.可以看到网址为: ...

  9. 修改Win10默认窗口背景色为护眼色的方法

    按Windows键+R,打开“运行”对话框,输入regedit: 修改[HKEY_CURRENT_USER\Control Panel\Colors]下Windows键值为 199 237 204, ...

  10. Codeforces 58E Expression (搜索)

    题意:给你一个可能不正确的算式a + b = c, 你可以在a,b,c中随意添加数字.输出一个添加数字最少的新等式x + y  = z; 题目链接 思路:来源于这片博客:https://www.cnb ...