【hdu6051】If the starlight never fade
Solution
神仙题qwq好吧我个人感觉是神仙题
这题其实有一个比较野路子的做法。。就是。。打表观察。。反正场上ckw大佬就是这样把这题A穿的%%%
然而实际上正解很神秘或者说很妙。。虽然说是不是用原根是。。套路?反正记录一下:
注意到\(P\)是奇质数,那么我们可以找到一个模\(P\)意义下的原根\(g\)
然后因为原根具有一些十分优秀的性质(存在唯一的\(a\)满足\(g^a=x\),其中\(x\in [1,p),a\in [1,p)\))所以我们可以开始快乐了
我们记\(x=g^a,y=g^b(a,b\in [1,p))\),那么原来求\(f(i)\)的式子可以变成:
(x+y)^i\equiv x^i(mod\ P)\Rightarrow (1+g^{b-a})^i\equiv 1(mod\ P)
\end{aligned}
\]
(就是两边同除以\(x^i\))
然后因为\(1+g^{b-a}>=2\),所以一定存在唯一一个\(k\in [1,p-1)\)满足:
\]
所以我们可以把上面的式子写成:
\]
也就是说\((P-1)|ki\)(因为原根的性质并且\(g^{P-1}\equiv 1 (mod\ P)\)),那么考虑\(k\)的取值,我们可以把\(k\)写成这样:
\]
然后因为\(k\in [1,P-1)\),所以\(s'\in [1,gcd(P-1,i))\),所以\(k\)总共有\(gcd(P-1,i)-1\)种取值
这个时候反过来想,这\(gcd(P-1,i)-1\)种取值,对应的是那么多个不同的\(k\),再反推回去对应的就是那么多个\(b-a\)的取值,也就是说当\(b\)固定的时候,也就是\(y\)固定的时候,有那么多个\(a\)可以取(也就是说有那么多个\(x\)可以取),即对于每一个\(y\),能取的\(x\)都有\(gcd(P-1,i)-1\)个
那所以\(f(i)=m(gcd(P-1,i)-1)\)
(可以这样一路推上去都是因为原根那个一一对应的性质,否则不能进行这么神秘的操作)
那所以我们的问题变成了求:
\]
这个时候我们就要开始大力化式子了qwq:
&\sum\limits_{i=1}^{P-1}i\cdot m(gcd(P-1,i)-1)\\
=&m(\sum\limits_{i=1}^{P-1}i\cdot gcd(P-1,i)-\sum\limits_{i=1}^{P-1}i)\\
=&m\sum\limits_{d|(P-1)}d\sum\limits_{d|i,1<=i<=P-1}i[gcd(P-1,i)=d]-\frac{m(P-1)(P-2)}{2}\\
=&m\sum\limits_{d|(P-1)}d^2\sum\limits_{i=1}^{\frac{P-1}{d}}i[gcd(\frac{P-1}{d},i)=1]\\
\end{aligned}
\]
这个时候我们需要一个黑科技:
&\sum\limits_{i=1}^{n}i[gcd(n,i)=1]\\
=&\frac{1}{2}\sum\limits_{i=1}^n (i+(n-i))[gcd(n,i)=1]\\
=&\frac{n}{2}\sum\limits_{i=1}^n [gcd(n,i)=1]\\
=&\frac{n\cdot \varphi(n)}{2}
\end{aligned}
\]
具体为什么的话就是。。如果说我们知道\(gcd(n,i)=1\),我们可以推出\(gcd(n,n-i)=1\),因为如果\(gcd(n,n-i)=x(x>1)\),那么\(i=n-(n-i)\)应该也是\(x\)的倍数,所以\(gcd(n,i)\neq 1\)(最少为\(x\)),所以矛盾,得证
然后我们把每一对\(i\)和\(n-i\)凑在一起就是上面那个式子的样子了
然后对于\(n=1\)的情况需要特判一下就是如果\(n=1\)那么应该是等于\(1\),也就是\(\frac{n\cdot \varphi(n)+[n=1]}{2}\)
然后有了这个黑科技我们就可以得到最终的式子:
\sum\limits_{i=1}^{P-1}if(i)
=m\sum\limits_{d|(P-1)}d^2\frac{\frac{P-1}{d}\cdot \varphi(\frac{P-1}{d})+[\frac{P-1}{d}=1]}{2}
\end{aligned}
\]
然后\(\varphi\)的话直接根号求,然后枚举一下因数就好了
代码大概长这个样子
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int MOD=1e9+7,inv2=500000004;
int n,m,P,ans,T;
int mul(int x,int y){return 1LL*x*y%MOD;}
int add(int x,int y){return (1LL*x+y)%MOD;}
int Phi(int x){
int ret=1;
for (int i=2;i*i<=x;++i){
if (x%i==0) ret=mul(ret,(i-1)),x/=i;
while (x%i==0) ret=mul(ret,i),x/=i;
}
if (x>1) ret=mul(ret,(x-1));
return ret;
}
int calc(int n){
int ret=0,x,tmp;
for (int i=1;i*i<=n;++i){
if (n%i) continue;
x=n/i;
tmp=add(mul(Phi(x),x),(x==1));
tmp=mul(1LL*i*i%MOD,tmp);
ret=add(ret,tmp);
if (i*i==n) continue;
x=i;
tmp=add(mul(Phi(x),x),(x==1));
tmp=mul(1LL*(n/i)*(n/i)%MOD,tmp);
ret=add(ret,tmp);
}
ret=mul(ret,inv2);
ret=add(ret,MOD-(1LL*(n+1)*n/2)%MOD);
return ret;
}
int main(){
#ifndef ONLINE_JUDGE
freopen("a.in","r",stdin);
#endif
int tmp,x;
scanf("%d",&T);
for (int o=1;o<=T;++o){
scanf("%d%d",&m,&P);
ans=calc(P-1);
printf("Case #%d: %d\n",o,mul(m,ans));
}
}
【hdu6051】If the starlight never fade的更多相关文章
- 关于【bootstrap】中,【tooltip】的不算bug的bug的个人看法
先说下遇到这个问题的背景吧. 就是在页面中有个div,这个div右上角(或者其他位置)有个 × 的图标(这个图标添加tooltip工具提示),光标移到这个图标时,触发tooltip,提示“点击移除”这 ...
- 使用【 ajax 】【 bootstrap 】显示出小窗口 详情内容 一些代码意思可以参考下一个文章
使用[ bootstrap ]显示出小窗口 详情内容 显示页面 <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional// ...
- 【Unity】4.6 灯光
分类:Unity.C#.VS2015 创建日期:2016-04-11 一.简介 灯光(Light,也叫光源)是每一个场景的重要组成部分,用于照亮场景和对象,从而让游戏具有自己的个性和风格,比如利用灯光 ...
- 【Unity】4.3 地形编辑器
分类:Unity.C#.VS2015 创建日期:2016-04-10 一.简介 Unity拥有功能完善的地形编辑器,支持以笔刷绘制的方式实时雕刻出山脉.峡谷.平原.高地等地形.Unity地形编辑器同时 ...
- C#设计模式总结 C#设计模式(22)——访问者模式(Vistor Pattern) C#设计模式总结 .NET Core launch.json 简介 利用Bootstrap Paginator插件和knockout.js完成分页功能 图片在线裁剪和图片上传总结 循序渐进学.Net Core Web Api开发系列【2】:利用Swagger调试WebApi
C#设计模式总结 一. 设计原则 使用设计模式的根本原因是适应变化,提高代码复用率,使软件更具有可维护性和可扩展性.并且,在进行设计的时候,也需要遵循以下几个原则:单一职责原则.开放封闭原则.里氏代替 ...
- 【转载】C# 开源库大全非常好
原文地址:http://m.blog.csdn.net/woddle/article/details/37311877 C#开源大全 商业协作和项目管理平台-TeamLab 网络视频会议软件-VMuk ...
- 【Bootstrap5】精细学习记录
[Bootstrap5]精细学习记录 Bootstrap模板 <!DOCTYPE html> <html> <head> <title>Bootstra ...
- Python高手之路【六】python基础之字符串格式化
Python的字符串格式化有两种方式: 百分号方式.format方式 百分号的方式相对来说比较老,而format方式则是比较先进的方式,企图替换古老的方式,目前两者并存.[PEP-3101] This ...
- 【原】谈谈对Objective-C中代理模式的误解
[原]谈谈对Objective-C中代理模式的误解 本文转载请注明出处 —— polobymulberry-博客园 1. 前言 这篇文章主要是对代理模式和委托模式进行了对比,个人认为Objective ...
随机推荐
- gitlab runner 填坑记
一.Gitlab Runner CI/CD 错误: Couldn't connect to Docker daemon at http+docker://localhost - is it runn ...
- Appium+python的单元测试框架unittest(2)——fixtures(转)
(原文:https://www.cnblogs.com/fancy0158/p/10046333.html) unittest提供的Fixtures用以在测试执行前和执行后进行必要的准备和清理工作,可 ...
- 转 gerrit
开发环境 https://blog.csdn.net/u013207966/article/details/79112740 先记录下我的开发环境以及要正确安装gerrit需要用到的工具: Redha ...
- 面向 Unity* 软件和虚拟现实的优化:运行时生成内容
优化游戏以实现高性能一直是游戏开发过程中的一个重要因素.虽然开发人员一直尝试将硬件推向极致,但当移动游戏成为主流时,优化技术变得尤为突出.Unity* 软件.Unreal* 等常见引擎最初都是面向 P ...
- Python模块random使用详情
python常用模块目录 1.random.random()#用于生成一个0到1的随机浮点数:0<= n < 1.0 import random mcw = random.random() ...
- TCP的三次握手(建立连接)和四次挥手(关闭连接)(转)
转自:(http://www.cnblogs.com/Jessy/p/3535612.html) 参照: http://course.ccniit.com/CSTD/Linux/reference/f ...
- lamp一键配置 --转自秋水
https://teddysun.com/lamp LAMP一键安装脚本 最后修改于:2015年11月08日 / 秋水逸冰 / 54,300 次围观 973 本脚本适用环境: 系统支持:CentOS/ ...
- 期中HTML代码及技术博客
<!DOCTYPE html> <html lang="en"> <head> <meta charset="U ...
- 201621123037 《Java程序设计》第7周学习总结
作业06-接口.内部类 1. 本周学习总结 以你喜欢的方式(思维导图或其他)归纳总结集合相关内容. 答: 思维导图: 其他-笔记: 2. 书面作业 1. ArrayList代码分析 1.1 解释Arr ...
- Objective - C 之协议
一.创建方法: 二.实现过程: 1.遵循协议: @protocol NurseWorkingProtocol <NSObject> //<> 表示遵守协议,创建时就有(Nu ...