链接

https://www.lydsy.com/JudgeOnline/problem.php?id=3131

思路

1.

函数值的素因子只有2、3、5、7

由他们组成的状态不多,爆搜的时候即使搜不对也没关系,我们只是缩小范围而已

所以不要管呢么多,搜到几万就差不多了,包含有可能的就行

2.

\(f[i][j][0/1]\)表示后i位,乘积为j,n的i位>=k(枚举1-9)?1:0

n的i位>k

\(f[i][j][1]+=f[i][j][1]+f[i][j][0]\)

n的i位<k

\(f[i][j][0]+=f[i][j][1]+f[i][j][0]\)

n的i位==k

\(f[i][j][1]+=f[i][j][1]\)

\(f[i][j][0]+=f[i][j][0]\)

3.

统计的时候,sort

之后想象成二维数组,每一行都是单调的,所以只用头做比较,堆维护就可以

单调队列没想过

错误$$感叹

dp调试真快啊,处理起来真恶心曹丹清啊

代码

#include <iostream>
#include <queue>
#include <map>
#include <cstdio>
#include <algorithm>
#include <vector>
#define ll long long
const ll N=5e5+7;
const ll mod=1e9+7;
using namespace std;
ll n,k,val[N],c[N],js;
ll f[15][N][2];
ll A,B,jl[15];
map<ll,ll> hasH;
struct Pair {
ll first,second;
Pair() {}
Pair(ll a,ll b) :
first(a),second(b) {}
};
bool operator < (Pair a,Pair b) {
if(a.second==js+1) return 0;
if(b.second==js+1) return 1;
return c[a.first]*c[a.second]<c[b.first]*c[b.second];
}
bool cmp(ll a,ll b) {
return a>b;
}
priority_queue<Pair> q;
int main() {
// freopen("a.in","r",stdin);
cin>>n>>k;
for(ll a=0,tmp_a=1;a<=39;++a,tmp_a*=2) {
if(tmp_a>n) break;
for(ll b=0,tmp_b=1;b<=25;++b,tmp_b*=3) {
if(tmp_a*tmp_b>n) break;
for(ll c=0,tmp_c=1;c<=17;++c,tmp_c*=5) {
if(tmp_a*tmp_b*tmp_c>n) break;
for(ll d=0,tmp_d=1;d<=14;++d,tmp_d*=7) {
if(tmp_a*tmp_b*tmp_c*tmp_d>n) break;
if(tmp_a*tmp_b*tmp_c*tmp_d==0) continue;
val[++js]=tmp_a*tmp_b*tmp_c*tmp_d;
}
}
}
}
sort(val+1,val+1+js);
js=unique(val+1,val+1+js)-val-1;
for(ll i=1;i<=js;++i) hasH[val[i]]=i;
ll AAA=0;
while(n) jl[++AAA]=n%10,n/=10;
f[0][1][1]=1;
for(ll i=1;i<=AAA;++i) {
for(ll j=1;j<=js;++j) {
for(ll k=1;k<=9;++k) {
if(val[j]%k==0) {
ll w=hasH[val[j]/k];
if(k==jl[i]) {
f[i][j][1]+=f[i-1][w][1];
f[i][j][0]+=f[i-1][w][0];
} else if(k<jl[i]) {
f[i][j][1]+=f[i-1][w][0]+f[i-1][w][1];
} else {
f[i][j][0]+=f[i-1][w][0]+f[i-1][w][1];
}
}
}
}
}
for(ll i=1;i<=js;++i) {
for(ll j=1;j<AAA;++j)
c[i]+=f[j][i][0]+f[j][i][1];
c[i]+=f[AAA][i][1];
}
sort(c+1,c+1+js,cmp);
for(ll i=1;i<=js;++i) q.push(Pair(i,1));
ll ans=0;
for(ll i=1;i<=k;++i) {
Pair x=q.top();
q.pop();
ans+=c[x.first]*c[x.second]%mod;
ans%=mod;
q.push(Pair(x.first,x.second+1));
}
cout<<ans<<"\n";
return 0;
}

bozoj3131: [Sdoi2013]淘金 数位dp的更多相关文章

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

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

  2. [SDOI2013]淘金 数位DP

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

  3. [Bzoj3131][Sdoi2013]淘金(数位dp)(优先队列)

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

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

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

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

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

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

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

  7. 数位DP学习笔记

    数位DP学习笔记 什么是数位DP? 数位DP比较经典的题目是在数字Li和Ri之间求有多少个满足X性质的数,显然对于所有的题目都可以这样得到一些暴力的分数 我们称之为朴素算法: for(int i=l_ ...

  8. Bzoj 3131 [Sdoi2013]淘金 题解

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

  9. 题解-SDOI2013 淘金

    题面 SDOI2013 淘金 有一个 \(X\).\(Y\) 轴坐标范围为 \(1\sim n\) 的范围的方阵,每个点上有块黄金.一阵风来 \((x,y)\) 上的黄金到了 \((f(x),f(y) ...

随机推荐

  1. 关于windows中80端口被占用

    很奇怪,windows7系统中的80端口被pid 为4 的system进程监听. 尝试关闭IIS功能,并在这期间进行过多次重启,均无效. 后来依照这篇文件解决的:https://www.jianshu ...

  2. netcore webapi统一配置跨域问题

    修改startup类中的configure方法

  3. BFS 搜索 蓝桥杯模拟赛

    题目链接:https://nanti.jisuanke.com/t/36117 这个题目想不到用广搜来做,一直在想深搜. 广搜的思路呢,是把最外圈不是黑色(不是0)的数 的位置 i 和 j 进队,赋值 ...

  4. Yii2返回以主键id为键名的数组

    branch.php <?php namespace app\models; use Yii; /** * This is the model class for table "bra ...

  5. 韩松毕业论文笔记-第六章-EFFICIENT METHODS AND HARDWARE FOR DEEP LEARNING

    难得跟了一次热点,从看到论文到现在已经过了快三周了,又安排了其他方向,觉得再不写又像之前读过的N多篇一样被遗忘在角落,还是先写吧,虽然有些地方还没琢磨透,但是paper总是这样吧,毕竟没有亲手实现一下 ...

  6. 从网站上扒网页,保存为file文件格式

    保存下来的页面总是有部分特效缺失,可是文件包里已经有好几个js文件了. 例如想保存易迅的搜索页面,条件筛选栏的按钮全部失效了,按钮-更多.多选等 都没有反应,搜索结果的鼠标悬浮显示完整信息也没有了. ...

  7. Python中*args和**kwargs 的简单使用

    # 在函数定义中使用*args和kwargs传递可变长参数. *args用作传递非命名键值可变长参数列表(位置参数); kwargs用作传递键值可变长参数列表# *args表示任何多个无名参数,它是一 ...

  8. socket聊天的业务逻辑

        一.主要思想:     1.如果用户A想要发消息给用户B,A需要将消息发送到一个服务器上,服务器接收到A发送的消息之后,再把消息发送给B,B接收到消息     2.当用户B断开连接时服务器不会 ...

  9. tomcat 启动时遇到org.apache.jasper.servlet.TldScanner.scanJars At least one JAR was scanned for TLDs yet contained no TLDs

    当发生这样的错误的时候 org.apache.jasper.servlet.TldScanner.scanJars At least one JAR was scanned for TLDs yet ...

  10. Qt QTextEdit根据行号移动光标

    QTextEdit* p = new QTextEdit; QTextBlock block = p->document()->findBlockByNumber(nLineNum); p ...