【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 ...
随机推荐
- Vue2 v-bind:href 中如何使用过滤器
<a class="topic_title" v-bind:href="info.id|getTitleHref" v-bind:title=" ...
- 【转】: 塞尔达组在GDC2017演讲的文字翻译:创新的勇气
大家好,我是藤林秀麿,以导演的身份参与<荒野之息>的制作,感谢大家的出席.我曾经作为设计者和导演制作了诸多塞尔达游戏(大地与时空之章.缩小帽.四支剑.幻影沙漏.天空之剑),回首望去,我已经 ...
- 在Windows2008下添加iscsi存储出现磁盘Offine(The disk is offine because of policy set by an adminstrator)的解决方法
打开CMD命令行输入如下命令: DISKPART.EXE DISKPART> san SAN Policy : Offline Shared DISKPART> san policy=On ...
- kubeadm构建k8s之Prometheus-operated监控(0.18.1)
介绍: 大家好,k8s的搭建有许多方式,也有许多快速部署的,为了简化部署的复杂度,官方也提供了开源的kubeadm快速部署,最新1.10.x版本已经可以实现部署集群, 如果你对k8s的原理已经非常了解 ...
- CSS Grid布局指南
简介 CSS Grid布局 (又名"网格"),是一个基于二维网格布局的系统,主要目的是改变我们基于网格设计的用户接口方式.如我们所知,CSS 总是用于网页的样式设置,但它并没有起到 ...
- eclipse启动一闪而退
eclipse启动一闪而退 打开eclipse,启动画面一闪而过退出. 解决方法: 以下每一步结束都重启eclipse一下,看能不能正常启动. 1. 在C:/WINDOWS/system32 系统文件 ...
- TP框架代码学习 学习记录 3.2.3
文件:think.class.php PHP提供register_shutdown_function()这个函数,能够在脚本终止前回调注册的函数,也就是当 PHP 程序执行完成后执行的函数.regis ...
- postman的巨坑 之 cookie
问题描述:一个后端接口,该接口需要校验登录态,登录态通过cookie中的一个传参k判断.在保证登录的前提下,调用机器A上的接口一直报“用户未登录”,调机器B上就没问题,于是开始排查问题. 解决过程: ...
- 网站UI分析
本次网站UI分析我选择的是我们石家庄铁道大学的网站,首先对于网站的分析建立在我经常使用鼠须的基础上,我可以很好的站在用户的角度来进行分析,否则对于你不熟悉的网站你可能是不能很好地体验到他的 结构. U ...
- HDU 5206 Four Inages Strategy 水题
题目链接: hdu:http://acm.hdu.edu.cn/showproblem.php?pid=5206 bc(中文):http://bestcoder.hdu.edu.cn/contests ...