[CSP-S模拟测试]:Six(数学)
题目传送门(内部题85)
输入格式
一个正整数$N$。
输出格式
一个数表示答案对$1000000007$取模后的结果
样例
样例输入1:
6
样例输出1:
28
样例输入2:
203021
样例输出2:
33628
样例输入3:
60357056536
样例输出3:
907882
样例输入4:
12156144
样例输出4:
104757552
数据范围与提示
样例解释:
第一组样例:$\{(2),(2,2),(2,2,3),(2,2,3,3),(2,3),(2,3,2),(2,3,2,3),(2,3,3),(2,3,3,2),(2,6),(2,6,3),(3),(3,2),(3,2,2),(3,2,2,3),(3,2,3),(3,2,3,2),(3,3),(3,3,2),(3,3,2,2),(3,6),(3,6,2),(6),(6,2),(6,2,3),(6,3),(6,3,2),(6,6)\}$
数据范围:
对于$32\%$的数据,$N$最多有$2$个不同的质因数
对于$56\%$的数据,$N$最多有$4$个不同的质因数
对于$100\%$的数据,$N\leqslant 10^{15}$且$N$最多有$6$个不同的质因数
题解
爆搜能有一半的分,在加上记忆化就$A$了……
记录一下二进制位下每一位$1$的个数状态,但是注意$1$可以是两个,所以还需要再开一个数记录一下是否已经达到了两个。
队列里的数最多只有$12$个,所以状态数也不多,$50,000$以内,用$map$记录就好了。
时间复杂度:$\Theta(\sqrt{N}+50,000)$。
期望得分:$100$分。
实际得分:$100$分。
代码时刻
#include<bits/stdc++.h>
using namespace std;
const int mod=1000000007;
map<pair<long long,long long>,int>dp;
long long N;
long long d[1001],fla[1001],top;
pair<long long,int> pri[10];
long long dfs(int x,long long res1,long long res2)
{
if(dp[make_pair(res1,res2)])return dp[make_pair(res1,res2)];
dp[make_pair(res1,res2)]=1;
for(int i=1;i<(1<<top);i++)
{
int sum=0;bool flag=0;
for(int j=1;j<(1<<top);j++)
{
if(!(i&j))continue;
if((res1>>j)&1)sum++;
if((res2>>j)&1)flag=1;
if(flag||sum>1)goto nxt;
}
if((res1>>i)&1)dp[make_pair(res1,res2)]=(dp[make_pair(res1,res2)]+d[i]*dfs(x+1,res1^(1LL<<i),res2|(1LL<<i))%mod)%mod;
else dp[make_pair(res1,res2)]=(dp[make_pair(res1,res2)]+d[i]*dfs(x+1,res1|(1LL<<i),res2)%mod)%mod;
nxt:;
}
return dp[make_pair(res1,res2)];
}
int main()
{
scanf("%lld",&N);
for(long long i=2;i*i<=N;i++)
if(!(N%i))
{
pri[++top].first=i;
while(!(N%i))
{
pri[top].second++;
N/=i;
}
}
if(N!=1)pri[++top]=make_pair(N,1);
d[0]=1;
for(int i=1;i<=top;i++)fla[1<<i-1]=i;
for(int i=1;i<(1<<top);i++)d[i]=d[i^(i&-i)]*pri[fla[i&-i]].second;
printf("%lld\n",dfs(0,0,0)-1);
return 0;
}
rp++
[CSP-S模拟测试]:Six(数学)的更多相关文章
- [考试反思]0729NOIP模拟测试10
安度因:哇哦. 安度因:谢谢你. 第三个rank1不知为什么就来了.迷之二连?也不知道哪里来的rp 连续两次考试数学都占了比较大的比重,所以我非常幸运的得以发挥我的优势(也许是优势吧,反正数学里基本没 ...
- 0823NOIP模拟测试赛后总结
考了两场感觉虚了... NOIP模拟测试30 分着考的. 就只有T2的美妙的暴力拿分了,60分rank10,挂了. T1是一道sb题,爆零了十分遗憾. 许多人都掉进了输出格式的坑里,C没大写.少个空格 ...
- Android单元测试与模拟测试详解
测试与基本规范 为什么需要测试? 为了稳定性,能够明确的了解是否正确的完成开发. 更加易于维护,能够在修改代码后保证功能不被破坏. 集成一些工具,规范开发规范,使得代码更加稳定( 如通过 phabri ...
- [开源]微信在线信息模拟测试工具(基于Senparc.Weixin.MP开发)
目前为止似乎还没有看到过Web版的普通消息测试工具(除了官方针对高级接口的),现有的一些桌面版的几个测试工具也都是使用XML直接请求,非常不友好,我们来尝试做一个“面向对象”操作的测试工具. 测试工具 ...
- 安装nginx python uwsgi环境 以及模拟测试
uwsgi帮助文档: http://uwsgi-docs-cn.readthedocs.io/zh_CN/latest/WSGIquickstart.html http://uwsgi-docs.re ...
- 【模拟】【数学】CSU 1803 2016 (2016湖南省第十二届大学生计算机程序设计竞赛)
题目链接: http://acm.csu.edu.cn/OnlineJudge/problem.php?id=1803 题目大意: 给定n,m(n,m<=109)1<=i<=n,1& ...
- 利用Python中的mock库对Python代码进行模拟测试
这篇文章主要介绍了利用Python中的mock库对Python代码进行模拟测试,mock库自从Python3.3依赖成为了Python的内置库,本文也等于介绍了该库的用法,需要的朋友可以参考下 ...
- 转 C#实现PID控制的模拟测试和曲线绘图
C#实现PID控制的模拟测试和曲线绘图 本文分两部分,一部分是讲PID算法的实现,另一部分是讲如何用动态的曲线绘制出PID运算的结果. 首先,PID算法的理论模型请参考自动控制理论,最早出现的是模 ...
- Mockito:一个强大的用于Java开发的模拟测试框架
https://blog.csdn.net/zhoudaxia/article/details/33056093 介绍 本文将介绍模拟测试框架Mockito的一些基础概念, 介绍该框架的优点,讲解应用 ...
- NOIP模拟测试1(2017081501)
好,今天是cgg第一次举行模拟测试,希望各位支持. 时间限制:2小时 题目链接: 题目一:水得都没名字了 题目二:车站 题目三:选数 不要觉得2小时太少,我的题目很良心,都很简单. 答案可以在模拟测试 ...
随机推荐
- C++中的自定义内存管理
1,问题: 1,new 关键字创建出来的对象位于什么地方? 1,位于堆空间: 2,有没有可能位于其它地方? 1,有: 2,通过一些方式可以使动态创建的对象位于静态存储区: 3,这个存储区在程序结束后释 ...
- css精灵图使用
1. 精灵技术的使用 CSS 精灵其实是将网页中的一些背景图像整合到一张大图中(精灵图),然而,各个网页元素通常只需要精灵图中不同位置的某个小图,要想精确定位到精灵图中的某个小图,就需要使用CSS的b ...
- Spring Boot解决无法访问图片的问题
找了很多Spring Boot项目访问图片的解决方式,发现都是配置的,有时配置了也没有用.然后自己研究了一种简单操作的方法. 1,在Spring Boot的static目录下创建一个新目录img(或者 ...
- Linux安装rabbitmq (解决guest无法登陆问题)
安装Erlang wget https://packages.erlang-solutions.com/erlang-solutions-1.0-1.noarch.rpm sudo rpm -Uvh ...
- AStar 启发函数设计(老物)
作为我出山的第一篇日志,怎么也得写篇对得起我身份和地位的文章吧? 先容我吐槽一下不小心发的贴图,那个只是我不小心收藏了隔壁兄弟班的课表就别大家这么热情的 BB 我感到很有压力,额,废话不多说,立刻进入 ...
- 3-关于ES的几个小疑问和解答
1.ES如何实现分布式 2.ES如何实现高实时 3.ES如何实现高扩展 4.ES7.x版本为何废弃type 5.搜索原理--知乎es
- 初探 -1 JavaScript
JavaScript 教程 JavaScript 是 Web 的编程语言. 所有现代的 HTML 页面都使用 JavaScript. JavaScript 非常容易学. 本教程将教你学习从初级到高级J ...
- 关于android studio从2.3升级到3.0以上可能会遇到的问题
请参考链接: http://blog.csdn.net/hylczp/article/details/60137958 gradle-3.3-all网盘下载地址: 链接:http://pan.baid ...
- AIX中磁带设备的使用
1.AIX中的磁带设备 在UNIX系统中,磁带一般具有如下共同属性. 磁带密度:指磁带使用时采用的低密度方式和高密度方式 回退操作:磁带一次写入操作完成后,是否立刻倒带,回退到磁带头部. 位置 ...
- vue父组件调用子组件方法、父组件向子组件传值、子组件向父组件传值
一.父组件调用子组件方法 父组件代码 parent.vue <template> <div> <button @click="parentFun" ...