题目

Description

终于,在众亲们的奋斗下,最终boss 恭介被关进了库特设计的密室。正当她们松了一口气时,这个世界却发生了天翻覆地的变化:地面开始下沉,天空开始变成血红色,海水沸腾……一幅世界末日的图景。美鱼从她手中的古籍《若山牧水诗歌集》中发现了原因:白鸟は かなしからずや 空の青 海のあをにも 染まずただよふ 。大(xia)意(shuo)就是狡猾的恭介在创造这个世界的时候就篡改了法则。而这个法则的起源,就是一只生死之间的猫。这个猫被关在一个黑盒子里,盒子里有两个毒气罐,如果有任意一个毒气罐被打开那么猫将会被杀死,法则也能得到纠正。然而外界能控制的仅仅是这两个毒气罐被打开的概率。假设第一个毒气罐被打开的概率为1/x,第二个毒气罐为1/y(x,y 为正整数),那么当两个概率和为1/(n!)时,猫将会被莫名其妙地杀死。现在美鱼想知道,有多少对(x,y)可以让猫被莫名其妙杀死。

Input

一行,一个正整数n

Output

一行,满足题意的(x,y)对数。

Sample Input

6

Sample Output

135

Data Constraint

对于30%的数据 n<=6

对于60%的数据 n<=50

对于100%的数据 n<=700000

分析

  • 来手推一下
  • T2:先移项:
    x,y,n!为正整数



    则 
    所以k|(n!)^2+n!k
    所以k|(n!)^2

  • 然后我们求约数个数就是
  • 首先如果求n的因子个数
  • 那么我们就是分解质因数
  • 至于怎么算每个质数p出现的次数,p的次数等于[n/p]+[n/p^2]+[n/p^3]+….这个可以看成是把1~n里p的倍数的数字都除掉一个p,这样的数有[n/p]个

  • 除完之后这[n/p]个数会变成1~[n/p],然后再除掉一个p,就是[[n/p]/p]个,然后变成1~[n/p^2],一直做下去。

代码

 1 #include<iostream>
 2 #include<cstdio>
 3 using namespace std;
 4 const long long md=1000000000;
 5 int n,sum,num=1,vis[700002];
 6 unsigned long long s[1000002],f[10000002];
 7 void calc(long long x)
 8 {
 9     int p=0;
10     for (int i=1;i<=num;i++)
11     {
12         s[i]=s[i]*x+p;
13         p=s[i]/md;
14         s[i]%=md;
15     }
16     while (p) s[++num]=p%md,p/=md;
17 }
18 int main ()
19 {
20     int n;
21     cin>>n;
22     s[1]=1;
23     long long k=1;
24     for (int i=2;i<=n;i++)
25     if (!vis[i])
26     {
27         for (int j=1;j<=n/i;j++) vis[i*j]=1;
28         int t=n/i,t1=0;
29         while (t) t1+=t,t/=i;
30         t1=t1*2+1;
31         if (k*t1>md)
32         {
33             calc(k);
34             k=1;
35         }
36         k*=t1;
37     }
38     if (k>1) calc(k);
39     printf("%lld",s[num]);
40     for (int i=num-1;i>=1;i--)
41       printf("%09lld",s[i]);
42 } 

转载于:https://www.cnblogs.com/zjzjzj/p/11360876.html

JZOJ 3456. 【NOIP2013模拟联考3】恭介的法则(rule)相关推荐

  1. JZOJ 3461. 【NOIP2013模拟联考5】小麦亩产一千八(kela)

    3461. [NOIP2013模拟联考5]小麦亩产一千八(kela) (Standard IO) Time Limits: 1000 ms  Memory Limits: 262144 KB  Det ...

  2. jzoj 3461. 【NOIP2013模拟联考5】小麦亩产一千八(math)

    3461. [NOIP2013模拟联考5]小麦亩产一千八 Description "有了金坷垃,肥料一袋能顶两袋撒,小麦亩产一千八,吸收两米下的氮磷钾--",话说HYSBZ(Hen ...

  3. JZOJ 3455. 【NOIP2013模拟联考3】库特的向量(code)

    题目 Description 从前在一个美好的校园里,有一只(棵)可爱的弯枝理树.她内敛而羞涩,一副弱气的样子让人一看就想好好疼爱她.仅仅在她身边,就有许多女孩子想和她BH,比如铃,库特,等等.不过, ...

  4. JZOJ 3468. 【NOIP2013模拟联考7】OSU!(osu)

    Description osu 是一款群众喜闻乐见的休闲软件. 我们可以把osu的规则简化与改编成以下的样子: 一共有n次操作,每次操作只有成功与失败之分,成功对应1,失败对应0,n次操作对应为一个长 ...

  5. JZOJ 3457. 【NOIP2013模拟联考3】沙耶的玩偶(doll)

    题目 Description 在美鱼和理树后援团拯救世界的同时,外表柔弱的理树也开始坚强起来,思考着离开这个世界的办法.误打误撞地,她遇上了正在教室破坏课桌打开迷宫入口的沙耶.沙耶告诉理树,这个世界的 ...

  6. [jzoj 3461]【NOIP2013模拟联考5】小麦亩产一千八 {Fibonacci数列}

    题目 Description "有了金坷垃,肥料一袋能顶两袋撒,小麦亩产一千八,吸收两米下的氮磷钾--",话说HYSBZ(Hengyang School for Boys & ...

  7. 【NOIP2013模拟联考5】军训(training)

    Description HYSBZ 开学了!今年HYSBZ 有n 个男生来上学,学号为1-n,每个学生都必须参加军训.在这种比较堕落的学校里,每个男生都会有Gi 个女朋友,而且每个人都会有一个欠扁值H ...

  8. 【NOIP2013模拟联考6】选课(select)

    Description 你真的认为选课是那么容易的事吗?HYSBZ的ZY同志告诉你,原来选课也会让人产生一种想要回到火星的感觉.假设你的一周有n天,那么ZY编写的选课系统就会给你n堂课.但是该系统不允 ...

  9. 【NOIP2013模拟联考5】休息(rest)

    Description 休息的时候,可以放松放松浑身的肌肉,打扫打扫卫生,感觉很舒服.在某一天,某LMZ 开始整理他那书架.已知他的书有n 本,从左到右按顺序排列.他想把书从矮到高排好序,而每一本书都 ...

最新文章

  1. 关于Socket通信客户端是否需要绑定端口号
  2. linuxmysql乱码
  3. ASP调用存储过程详解。
  4. python教程:sum(1-2+3-4+...+99)
  5. Linux进程间通信分类 以及 pipe的原理实现
  6. svn的安装和服务器搭建, 汉化以及多项目管理
  7. linux阻塞原语,Linux中的内存原语是什么?
  8. 洛谷 4178 Tree——点分治
  9. 使用postman调试web API
  10. labimage 怎样旋转图片_隔断墙见多了,头次见能180旋转任意移动,还多出一面墙来储物...
  11. 多目标优化算法:多目标樽海鞘算法MOSSA(提供MATLAB源码)
  12. 彩扩机项目--两个状态指示灯,巧妙的使用二极管的特性,以及物理按键的滤波...
  13. 手把手教你用jQuery Mobile做相册
  14. LightOJ - 1395
  15. html四舍五入函数,Javascript四舍五入(Math.round()与Math.pow())
  16. C 生化危机 SDUT
  17. 调制与解调(一)——模拟调制系统调制解调方法
  18. 1.1 海思3518 H264编码
  19. 公网IP 申请 与 配置
  20. aix curl 安装_aix6.1 openssh安装

热门文章

  1. apache配置防盗链
  2. IEEE ISO/IEC简介
  3. aehyok.com的成长之路一——开篇
  4. JEE_Ajax技术
  5. jQuery插件开发全解析(转)
  6. 使用Ant制作发布版本
  7. 面试题40. 最小的k个数
  8. Varnish——CDN推送平台(web页面批量清除缓存)
  9. ubuntu 安装搜狗输入法_超详细的ubuntu14.04安装搜狗输入法教程
  10. 美国纽约法院定于明年3月对BitMEX前高管进行审判