bzoj3560: DZY Loves Math V
时间:2019-03-25
本文章向大家介绍bzoj3560: DZY Loves Math V,主要包括bzoj3560: DZY Loves Math V使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
这个题不应该想不到的啊
考虑对于每个质因数分别计算答案,最后再乘起来
那么每个答案就是phi(p^(sigema(1~n)i bi)) bi表示这个位置用了多少当前质因数,显然我们可以算出bi的上界ci,即在ai中pi出现了几次
而函数值非常简单,就是p^(k-1)*(p-1),而根据乘法分配律可以得到答案其实就是product(1~n)i (sigema(1~ci)j p^j)*(p-1)/p
因为一个数的因数个数不会超过30,直接暴力算就可以了
还有一个问题,就是phi(1)的情况,我们可以对柿子进行一点操作令它合法,变成((sigema(1~ci)j p^j)-1)*(p-1)/p+1就好了
#include<cstdio> #include<iostream> #include<cstring> #include<cstdlib> #include<algorithm> #include<cmath> using namespace std; typedef long long LL; const int _=1e2; const int maxn=1e5+_; const int maxp=10+5; const LL mod=1e9+7; LL quick_pow(LL A,int p) { LL ret=1; while(p!=0) { if(p%2==1)ret=ret*A%mod; A=A*A%mod;p/=2; } return ret; } struct node{int p,c;}p[maxp*maxn];int plen; bool cmp(node n1,node n2){return n1.p<n2.p;} void divi(int n,int w) { for(int i=2;i*i<=n;i++) if(n%i==0) { p[++plen].p=i; while(n%i==0)n/=i,p[plen].c++; } if(n!=1)p[++plen].p=n,p[plen].c=1; } int a[maxn]; int main() { int n; scanf("%d",&n); for(int i=1;i<=n;i++) { scanf("%d",&a[i]); divi(a[i],i); } sort(p+1,p+plen+1,cmp); LL ans=1; int l,r=0; while(l<=plen) { l=r+1; while(r<plen&&p[r+1].p==p[l].p)r++; LL num=1; for(int i=l;i<=r;i++) { LL sum=0,d=1; for(int j=0;j<=p[i].c;j++) { sum+=d;if(sum>=mod)sum-=mod; d=d*p[l].p%mod; } num=num*sum%mod; } num=(num+mod-1)*(p[l].p-1)%mod*quick_pow(p[l].p,mod-2)%mod+1; ans=ans*num%mod; } printf("%lld\n",ans); return 0; }
- 初识 Fuzzing 工具 WinAFL
- 如何使用Oozie API接口向Kerberos环境的CDH集群提交Spark2作业
- 如何编译及使用TPC-DS生成测试数据
- ASP.NET MVC编程——缓存
- ASP.NET MVC编程——错误处理与日记
- Jenkins 未授权远程代码执行漏洞(CVE-2017-1000353)
- ASP.NET MVC编程——路由
- ASP.NET MVC编程——模型
- ASP.NET MVC编程——验证、授权与安全
- 如何使用SAML配置CDSW的身份验证
- ASP.NET MVC编程——控制器
- ASP.NET MVC编程——视图
- 设计原则
- Entity Framework——配置文件设置
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- JVM系列之:JIT中的Virtual Call
- 速读原著-UnixLinux基础(五)
- Celery在Django中的简单应用
- 基于数据库Binlog 的业务系统操作日志实现方案(阿里中间件Canal)
- 速读原著-UnixLinux基础(三)
- 手撸实现UDP和TCP通信
- 3分钟短文 | PHP获取函数参数名,和类定义的常量,都要反射!
- 3分钟短文 | Laravel模型获取最后一条插入记录的ID编号
- 3分钟短文 | Laravel 模型的get find first方法你分清咋用了吗?
- 3分钟短文 ! Laravel拼装SQL子查询的最佳实现
- 对java中的泛型的理解
- 解决elasticsearch集群Unassigned Shards无法reroute的问题
- 一次系统扩容引起的elasticsearch故障及恢复
- WPF SharpDx 性能优化方法
- java中的reference(四): WeakReference的应用--ThreadLocal源码分析