3131: [Sdoi2013]淘金


Time Limit: 30 Sec  Memory Limit: 256 MB
Submit: 847  Solved: 423
[Submit][Status][Discuss]

Description


小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)取模之后的答案。

Input


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

Output


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

Sample Input


 

Sample Output



HINT


N < = 10^12 ,K < = 100000

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

题解:


发现题目变成了scoi2012某道题简化版。

可以很轻松愉快的发现如果是各位数字相乘,质因子只有2357,又限制最多只有12位,所以发现预处理出来可以得到数就一万多个。

于是我们就轻松愉快的处理出了,g[i][j],i位数乘积为j的数出现次数。

j很大,但数只有一万多个,hash走起。

好,数位dp转移也很明显了,不想太讲太多,f[i][j]   i位数内乘积为j的数出现次数。。

好了,发现转移完毕后我们要找所有数次数相乘的前k大。

暴力枚举所有数是10000^2的。

我们发现每个数都是从第一个数转移起走,那么我们对于10000多个数把它记录成三元组(i,j,k)当前价值为i,当前数位j,匹配数位k,每次选了当前三元组后

又加入(v[j] * v[k + 1],j,k + 1)放入优先队列即可

AC代码:


# include <iostream>
# include <cstdio>
# include <queue>
# include <algorithm>
using namespace std;
typedef long long LL;
const int N = 5e4 + ;
const int M = 2e5 + ;
const LL mod = 1e9 + ;
LL g[][N],state[N],s[N],n,K;int dt,hs[M],data[],len;
struct node{
LL val;int id;
bool operator <(const node & other)const{return val > other.val;}
}a[N];
struct e{
LL w;int id,Id;
bool operator <(const e & other)const{return w < other.w;}
};
priority_queue<e>st;
void insert(LL s)
{
int p = s % M;
while(hs[p])
{
if(state[hs[p]] == s)return;
p++;
if(p == M)p = ;
}
hs[p] = ++dt;state[hs[p]] = s;
}
int id(LL s)
{
int p = s % M;
while(hs[p])
{
if(state[hs[p]] == s)return hs[p];
p++;
if(p == M)p = ;
}
return ;
}
void init()
{
for(int i = ;i <= ;i++)
{
insert(i);
g[][id(i)] = ;
}
for(int i = ;i <= ;i++)
{
int tmp = dt;
for(int j = ;j <= tmp;j++)
{
LL x = state[j];int v = j,u;
for(LL k = ;k <= ;k++)
{
insert(x * k);
u = id(x * k);
g[i + ][u] += g[i][v];
}
}
}
}
LL dfs(int now,LL K,bool lim,bool first)
{
if(!now)return K == && !first;
if(!lim && !first)return g[now][id(K)];
LL ret = ;int p = lim ? data[now] : ;
if(first)ret += dfs(now - ,K,false,first);
for(LL i = ;i <= p;i++)if(!(K % i))
ret += dfs(now - ,K / i,lim && i == p,false);
return ret;
}
void calc(LL k)
{
len = ;
while(k)data[++len] = k % ,k /= ;
for(int i = ;i <= dt;i++)a[i] = (node){dfs(len,state[i],true,true),i};
sort(a + ,a + dt + );
for(int i = ;i <= dt;i++)if(a[i].val)
st.push((e){a[i].val * a[].val,i,});
}
int main()
{
init();
scanf("%lld %lld",&n,&K);
calc(n);LL ans = ,ret;int u,v;
while(K && !st.empty())
{
(ans += st.top().w) %= mod;u = st.top().id;v = st.top().Id + ;
st.pop();
if(!a[v].val)continue;
st.push((e){a[u].val * a[v].val,u,v});
K--;
}
printf("%lld\n",ans);
}

[Bzoj3131][Sdoi2013]淘金(数位dp)(优先队列)的更多相关文章

  1. [SDOI2013]淘金 数位DP

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

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

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

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

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

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

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

  5. [BZOJ3131] [Sdoi2013]淘金

    [BZOJ3131] [Sdoi2013]淘金 Description 小Z在玩一个叫做<淘金者>的游戏.游戏的世界是一个二维坐标.X轴.Y轴坐标范围均为1..N.初始的时候,所有的整数坐 ...

  6. bzoj千题计划268:bzoj3131: [Sdoi2013]淘金

    http://www.lydsy.com/JudgeOnline/problem.php?id=3131 如果已知 s[i]=j 表示有j个<=n数的数码乘积=i 那么就会有 s[a1]*s[a ...

  7. 【数位dp】bzoj3131: [Sdoi2013]淘金

    思路比较自然,但我要是考场上写估计会写挂:好像被什么不得了的细节苟住了?…… Description 小Z在玩一个叫做<淘金者>的游戏.游戏的世界是一个二维坐标.X轴.Y轴坐标范围均为1. ...

  8. [您有新的未分配科技点]数位DP:从板子到基础(例题 bzoj1026 windy数 bzoj3131 淘金)

    只会统计数位个数或者某种”符合简单规律”的数并不够……我们需要更多的套路和应用 数位dp中常用的思想是“分类讨论”思想.下面我们就看一道典型的分类讨论例题 1026: [SCOI2009]windy数 ...

  9. bzoj 3131 [Sdoi2013]淘金(数位dp)

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

随机推荐

  1. Java编译时根据调用该方法的类或对象所属的类决定

    class Base{     int x = 1;     static int y = 2; } class Subclass extends Base{     int x = 4;     i ...

  2. 洛谷 P1918 保龄球

    题目描述 DL 算缘分算得很烦闷,所以常常到体育馆去打保龄球解闷.因为他保龄球已经打了几十年了,所以技术上不成问题,于是他就想玩点新花招. DL 的视力真的很不错,竟然能够数清楚在他前方十米左右每个位 ...

  3. H3C S5024P交换机 vlan实验

    H3C S5024P交换机第二次vlan实验 实验1 与交换机端口G0/1和G0/2相连的PC1与PC2属于VLAN 1,与G0/3和G0/4相连的PC3和PC4属于VLAN 2,PC1.PC2.PC ...

  4. laravel扩展包服务提供者的注册的两种方式

    一. 包自动发现 在 Laravel 应用的配置文件 config/app.php 中,providers 配置项定义了一个会被 Laravel 加载的服务提供者列表.当安装完新的扩展包后,在老版本中 ...

  5. Python 中print 和return 的区别

    1.print() print()函数的作用是输出数据到控制台,就是打印在你能看到的界面上. 2.return return语句[表达式]退出函数,选择性地向调用方返回一个表达式.不带参数值的retu ...

  6. Windows下使用ffmpeg与java实现截取视频缩略图

    [ffmpeg.exe可执行文件获取]: 网上搜索后得到编译好的ffmpeg文件,下载地址:http://download.csdn.net/source/453719 [安装]: 将下载的文件解压, ...

  7. Springboot 配置文件与对象之间进行映射之@ConfigurationProperties

    一.将配置文件与实体类绑定1.1.将yaml配置文件的属性映射到Javabean中1.1.1.yaml配置文件注意:键值对的语法,键:之后必须要有空格 1.1.2.Javabean 定义注意:java ...

  8. activiti工作流学习链接

    首页: http://www.activiti.org/书籍: activiti in action  入门demo: kft-activiti-demo   http://www.oschina.n ...

  9. Multiplication Puzzle POJ - 1651

    解法 区间dp例题,长度从2开始到n结束起点从1到n,中间枚举的时候是看着左端点右端点与中点的乘积 代码 #include <iostream> #include <cstring& ...

  10. 洛谷P1035 级数求和

    #include <iostream> using namespace std; int main(){ long k,i; cin >> k; double s=0.0; ; ...