poj 3696 The Luckiest number 欧拉函数在解a^x=1modm的应用
2024-09-04 04:34:15
题意:
给一个L,求长度最小的全8数满足该数是L的倍数。
分析:
转化为求方程a^x==1modm。
之后就是各种数学论证了。
代码:
//poj 3696
//sep9
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
ll L;
ll factor[65536]; ll mul(ll x,ll y,ll p)
{
ll ret=0;
while(y){
if(y&1)
ret=(ret+x)%p;
x=(x+x)%p;
y>>=1;
}
return ret;
} ll pow(ll x,ll n,ll m)
{
ll ret=1;
x%=m;
while(n){
if(n&1)
ret=mul(ret,x,m);
x=mul(x,x,m);
n>>=1;
}
return ret;
} ll euler(ll n)
{
ll ret=n;
for(ll i=2;i*i<=n;++i)
if(n%i==0){
ret=ret/i*(i-1);
while(n%i==0)
n/=i;
}
if(n>1)
ret=ret/n*(n-1);
return ret;
} ll cal()
{
ll m=L*9;
for(int i=0;i<3;++i)
if(m%2==0)
m/=2;
else
break;
if(m%2==0||m%5==0)
return 0;
ll phi=euler(m);
int idx=0;
for(ll i=1;i*i<=phi;++i)
if(phi%i==0)
factor[idx++]=i,factor[idx++]=phi/i;
sort(factor,factor+idx);
for(int i=0;i<idx;++i)
if(pow(10,factor[i],m)==1)
return factor[i];
} int main()
{
int cases=0;
while(scanf("%I64d",&L)==1&&L){
printf("Case %d: ",++cases);
printf("%I64d\n",cal());
}
return 0;
}
最新文章
- [转]java.lang.OutOfMemoryError: PermGen space及其解决方法
- Java多线程与并发库高级应用-传统线程同步通信技术
- phpstorm8.0汉化版下载
- hlg 2130 状压dp
- nyoj-----127星际之门(一)
- jquery 事件绑定以及解绑定
- NYOJ-235 zb的生日 AC 分类: NYOJ 2013-12-30 23:10 183人阅读 评论(0) 收藏
- iOS点击cell查看大图,点击大图还原小图-b
- 计算几何(凸包模板):HDU 1392 Surround the Trees
- Ruby设计模式透析之 —— 适配器(Adapter)
- HTML5 CSS3 精美案例 : 达到VCD盒个性幻灯片
- (@WhiteTaken)设计模式学习——组合模式
- 通过Azure Powershell获取asm及arm虚拟机的配置信息
- ad 线束和网络
- 分库工具mysql
- solr如何让全词匹配结果在最前面
- 2018.4.26 Mac安装Redis5.0.3版本服务器
- 求二叉树中第K层结点的个数
- PhoneGap安装配置
- C# Socket编程 笔记,Socket 详解,入门简单
热门文章
- Ajaxterm + nginx 实现一个简单的堡垒机
- IE8的 JS 引擎如此不堪(二) - 解决方案
- [php] 解析JSON字符串
- TDD开发案例
- [UOJ198]时空旅行
- 【枚举】bzoj1643 [Usaco2007 Oct]Bessie&#39;s Secret Pasture 贝茜的秘密草坪
- 【欧拉函数表】POJ2478-Farey Sequence
- Problem A: 零起点学算法16——鸡兔同笼
- WPF附加属性的Set函数不调用的问题
- 用Qemu模拟vexpress-a9 (五) --- u-boot引导kernel,device tree的使用