题目传送门(内部题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(数学)的更多相关文章

  1. [考试反思]0729NOIP模拟测试10

    安度因:哇哦. 安度因:谢谢你. 第三个rank1不知为什么就来了.迷之二连?也不知道哪里来的rp 连续两次考试数学都占了比较大的比重,所以我非常幸运的得以发挥我的优势(也许是优势吧,反正数学里基本没 ...

  2. 0823NOIP模拟测试赛后总结

    考了两场感觉虚了... NOIP模拟测试30 分着考的. 就只有T2的美妙的暴力拿分了,60分rank10,挂了. T1是一道sb题,爆零了十分遗憾. 许多人都掉进了输出格式的坑里,C没大写.少个空格 ...

  3. Android单元测试与模拟测试详解

    测试与基本规范 为什么需要测试? 为了稳定性,能够明确的了解是否正确的完成开发. 更加易于维护,能够在修改代码后保证功能不被破坏. 集成一些工具,规范开发规范,使得代码更加稳定( 如通过 phabri ...

  4. [开源]微信在线信息模拟测试工具(基于Senparc.Weixin.MP开发)

    目前为止似乎还没有看到过Web版的普通消息测试工具(除了官方针对高级接口的),现有的一些桌面版的几个测试工具也都是使用XML直接请求,非常不友好,我们来尝试做一个“面向对象”操作的测试工具. 测试工具 ...

  5. 安装nginx python uwsgi环境 以及模拟测试

    uwsgi帮助文档: http://uwsgi-docs-cn.readthedocs.io/zh_CN/latest/WSGIquickstart.html http://uwsgi-docs.re ...

  6. 【模拟】【数学】CSU 1803 2016 (2016湖南省第十二届大学生计算机程序设计竞赛)

    题目链接: http://acm.csu.edu.cn/OnlineJudge/problem.php?id=1803 题目大意: 给定n,m(n,m<=109)1<=i<=n,1& ...

  7. 利用Python中的mock库对Python代码进行模拟测试

    这篇文章主要介绍了利用Python中的mock库对Python代码进行模拟测试,mock库自从Python3.3依赖成为了Python的内置库,本文也等于介绍了该库的用法,需要的朋友可以参考下     ...

  8. 转 C#实现PID控制的模拟测试和曲线绘图

    C#实现PID控制的模拟测试和曲线绘图   本文分两部分,一部分是讲PID算法的实现,另一部分是讲如何用动态的曲线绘制出PID运算的结果. 首先,PID算法的理论模型请参考自动控制理论,最早出现的是模 ...

  9. Mockito:一个强大的用于Java开发的模拟测试框架

    https://blog.csdn.net/zhoudaxia/article/details/33056093 介绍 本文将介绍模拟测试框架Mockito的一些基础概念, 介绍该框架的优点,讲解应用 ...

  10. NOIP模拟测试1(2017081501)

    好,今天是cgg第一次举行模拟测试,希望各位支持. 时间限制:2小时 题目链接: 题目一:水得都没名字了 题目二:车站 题目三:选数 不要觉得2小时太少,我的题目很良心,都很简单. 答案可以在模拟测试 ...

随机推荐

  1. C++中的自定义内存管理

    1,问题: 1,new 关键字创建出来的对象位于什么地方? 1,位于堆空间: 2,有没有可能位于其它地方? 1,有: 2,通过一些方式可以使动态创建的对象位于静态存储区: 3,这个存储区在程序结束后释 ...

  2. css精灵图使用

    1. 精灵技术的使用 CSS 精灵其实是将网页中的一些背景图像整合到一张大图中(精灵图),然而,各个网页元素通常只需要精灵图中不同位置的某个小图,要想精确定位到精灵图中的某个小图,就需要使用CSS的b ...

  3. Spring Boot解决无法访问图片的问题

    找了很多Spring Boot项目访问图片的解决方式,发现都是配置的,有时配置了也没有用.然后自己研究了一种简单操作的方法. 1,在Spring Boot的static目录下创建一个新目录img(或者 ...

  4. Linux安装rabbitmq (解决guest无法登陆问题)

    安装Erlang wget https://packages.erlang-solutions.com/erlang-solutions-1.0-1.noarch.rpm sudo rpm -Uvh ...

  5. AStar 启发函数设计(老物)

    作为我出山的第一篇日志,怎么也得写篇对得起我身份和地位的文章吧? 先容我吐槽一下不小心发的贴图,那个只是我不小心收藏了隔壁兄弟班的课表就别大家这么热情的 BB 我感到很有压力,额,废话不多说,立刻进入 ...

  6. 3-关于ES的几个小疑问和解答

    1.ES如何实现分布式 2.ES如何实现高实时 3.ES如何实现高扩展 4.ES7.x版本为何废弃type 5.搜索原理--知乎es

  7. 初探 -1 JavaScript

    JavaScript 教程 JavaScript 是 Web 的编程语言. 所有现代的 HTML 页面都使用 JavaScript. JavaScript 非常容易学. 本教程将教你学习从初级到高级J ...

  8. 关于android studio从2.3升级到3.0以上可能会遇到的问题

    请参考链接: http://blog.csdn.net/hylczp/article/details/60137958 gradle-3.3-all网盘下载地址: 链接:http://pan.baid ...

  9. AIX中磁带设备的使用

    1.AIX中的磁带设备 在UNIX系统中,磁带一般具有如下共同属性. 磁带密度:指磁带使用时采用的低密度方式和高密度方式   回退操作:磁带一次写入操作完成后,是否立刻倒带,回退到磁带头部.   位置 ...

  10. vue父组件调用子组件方法、父组件向子组件传值、子组件向父组件传值

      一.父组件调用子组件方法 父组件代码  parent.vue <template> <div> <button @click="parentFun" ...