【数位dp】bzoj3131: [Sdoi2013]淘金
思路比较自然,但我要是考场上写估计会写挂;好像被什么不得了的细节苟住了?……
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
题目分析
关键问题在于处理乘积为$i$的数的个数。这个还是一个思路比较自然的dp的。
但是好像有什么不得了的细节恶心到了。
1 // luogu-judger-enable-o2 2 #include<bits/stdc++.h> 3 typedef long long ll; 4 const int MO = 1e9+7; 5 6 ll sum[15003]; 7 struct node 8 { 9 int x,y; 10 ll val; 11 node(int a=0, int b=0):x(a),y(b) {} 12 bool operator < (node a) const 13 { 14 return 1ll*sum[x]*sum[y] < 1ll*sum[a.x]*sum[a.y]; 15 } 16 }; 17 int d[103]; 18 ll n,k,cnt,dct,ans; 19 ll f[103][15003][2]; 20 std::map<int, ll> mp1; 21 std::map<ll, int> mp2; 22 std::priority_queue<node> q; 23 24 bool cmp(ll a, ll b){return a > b;} 25 int main() 26 { 27 scanf("%lld%lld",&n,&k); 28 for (ll i=1; i<=n; i*=2ll) 29 for (ll j=1; i*j<=n; j*=3ll) 30 for (ll k=1; i*j*k<=n; k*=5ll) 31 for (ll l=1; i*j*k*l<=n; l*=7ll) 32 mp1[++cnt] = i*j*k*l, mp2[i*j*k*l] = cnt; 33 //MAXcnt==14672 34 for (ll x=n; x; x/=10) d[++dct] = x%10; 35 for (int i=1; i<=9; i++) f[1][mp2[i]][i > d[1]]++; 36 for (int t=2; t<=dct; t++) 37 for (int i=1; i<=cnt; i++) 38 for (int j=1; j<=9; j++) 39 { 40 ll num = mp1[i]; 41 if (num%j) continue; 42 int lb = mp2[num/j]; 43 if (j < d[t]) 44 f[t][i][0] += f[t-1][lb][0]+f[t-1][lb][1]; 45 else if (j > d[t]) 46 f[t][i][1] += f[t-1][lb][0]+f[t-1][lb][1]; 47 else f[t][i][0] += f[t-1][lb][0], f[t][i][1] += f[t-1][lb][1]; 48 } 49 for (int j=1; j<=cnt; j++) 50 for (int i=1; i<=dct; i++) 51 sum[j] += f[i][j][0]+(i==dct?0:f[i][j][1]); 52 // for (int j=1; j<=cnt; j++) 53 // for (int i=1; i<=dct; i++) 54 // sum[j] += f[i][j][0]+(i==dct)?0:f[i][j][1]; //这两个答案是不一样的?…… 55 k = std::min(k, cnt*cnt); 56 std::sort(sum+1, sum+cnt+1, cmp); 57 for (int i=1; i<=cnt; i++) q.push(node(i, 1)); 58 while (k--) 59 { 60 node tt = q.top(); 61 q.pop(), ans = (ans+sum[tt.x]*sum[tt.y]%MO)%MO; 62 if (tt.y==cnt) continue; 63 q.push(node(tt.x, tt.y+1)); 64 } 65 printf("%lld\n",ans); 66 return 0; 67 }
upd:原来$?:$的优先级比$+$低……怪不得很多人写的时候把双目括起来
END
转载于:https://www.cnblogs.com/antiquality/p/9818964.html
【数位dp】bzoj3131: [Sdoi2013]淘金相关推荐
- bzoj千题计划268:bzoj3131: [Sdoi2013]淘金
http://www.lydsy.com/JudgeOnline/problem.php?id=3131 如果已知 s[i]=j 表示有j个<=n数的数码乘积=i 那么就会有 s[a1]*s[a ...
- 【BZOJ3131】淘金(SDOI2013)-数位DP+优先队列
测试地址:淘金 做法: 本题需要用到数位DP+优先队列. 令 F ( x , y ) F(x,y) F(x,y)为坐标 ( x , y ) (x,y) (x,y)上的黄金数目,那么: F ( x , ...
- bzoj 3131 [Sdoi2013]淘金(数位dp)
题目描述 小Z在玩一个叫做<淘金者>的游戏.游戏的世界是一个二维坐标.X轴.Y轴坐标范围均为1..N.初始的时候,所有的整数坐标点上均有一块金子,共N*N块. 一阵风吹过,金子的位置发生了 ...
- 数位dp从会打模板到不会打模板
打了几个数位$dp$,发现自己除了会打模板之外没有任何长进,遇到非模板题依然什么都不会 那么接下来这篇文章将介绍如何打模板(滑稽) 假设我们要处理$l----r$ 采用记忆化搜索的方式,枚举$< ...
- 【bzoj 1833】【codevs 1359】 [ZJOI2010]count 数字计数(数位dp)
1833: [ZJOI2010]count 数字计数 Time Limit: 3 Sec Memory Limit: 64 MB Submit: 2774 Solved: 1230 [Submit ...
- 不要62 ---数位DP
题意:求m到n中不含62和4的数的个数. 题目链接 思路:数位dp模板求满足的数字或不满足的数字,刚学,就求不满足的数. #include<stdio.h> #include<str ...
- bzoj 3598 [ Scoi 2014 ] 方伯伯的商场之旅 ——数位DP
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3598 数位DP...东看西看:http://www.cnblogs.com/Artanis/ ...
- cojs 简单的数位DP 题解报告
首先这道题真的是个数位DP 我们考虑所有的限制: 首先第六个限制和第二个限制是重复的,保留第二个限制即可 第五个限制在转移中可以判断,不用放在状态里 对于第一个限制,我们可以增加一维表示余数即可 对于 ...
- 数位DP 不断学习中。。。。
1, HDU 2089 不要62 :http://acm.hdu.edu.cn/showproblem.php?pid=2089 题意:不能出现4,或者相邻的62, dp[i][0],表示不存在不吉 ...
最新文章
- 干货 | 深度学习检测小目标常用方法
- Java 设计模式 之 中介者模式(Mediator)
- R语言因子转数值类型
- HTTP面试题都在这里
- 技术动态 | 知识图谱构建的研究已走入下半场,但大规模落地应用仍需时间
- Android下 布局加边框 指定背景色 半透明
- .NET开发人员十大必备下载工具
- 漫画:为什么生僻字计算机上打不出来,或者打出来也无法显示呢?
- Security+ 学习笔记53 供应链风险
- NOD32 2.7、3.0最新升级ID 每天实时同步更新
- NetTiers学习笔记12---deepload的递归
- 「 C++ MFC 」“设置线程运行多媒体定时器”教程
- 黑客入侵Windows XP的几个常用方法
- UDS与DoIp整理
- BLE低功耗蓝牙协议栈
- 软件测试场景 例子,软件测试基础(六)用例设计方法之场景法
- Android 视频裁剪 (含裁剪 View)
- 2018年北京信息科技大学第十届程序设计竞赛暨ACM选拔赛 C:计算几何
- 修改windows自带录音机录音时间方法
- 如何在Microsoft Word中转录音频
热门文章
- GDP平减指数整理 (1978-2017年)
- 从两边往中间流水灯c语言程序,单片机控制LED流水灯从中间向两边,从两边向中间...
- 最大公约数、最小公倍数与算术基本定理
- 杰理之爬锡效果不理想【篇】
- 2024年天津农学院专升本拟招生专业及停招酒店管理、市场营销通知
- 主频是反映计算机的什么指标,主频是反映计算机什么的性能指标
- 如何给青少年讲海龟绘图的内容:第一节课
- 纸的大小图解_常用纸张尺寸及示意图(A0,A1A3,A4,A5A8)数据源维基百科.PDF
- thc-ssl-dos攻击https站点
- iMeta | 中国农业大学汪杰组“塑料际”微生物最新综述