题意:

给一个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;
}

最新文章

  1. [转]java.lang.OutOfMemoryError: PermGen space及其解决方法
  2. Java多线程与并发库高级应用-传统线程同步通信技术
  3. phpstorm8.0汉化版下载
  4. hlg 2130 状压dp
  5. nyoj-----127星际之门(一)
  6. jquery 事件绑定以及解绑定
  7. NYOJ-235 zb的生日 AC 分类: NYOJ 2013-12-30 23:10 183人阅读 评论(0) 收藏
  8. iOS点击cell查看大图,点击大图还原小图-b
  9. 计算几何(凸包模板):HDU 1392 Surround the Trees
  10. Ruby设计模式透析之 —— 适配器(Adapter)
  11. HTML5 CSS3 精美案例 : 达到VCD盒个性幻灯片
  12. (@WhiteTaken)设计模式学习——组合模式
  13. 通过Azure Powershell获取asm及arm虚拟机的配置信息
  14. ad 线束和网络
  15. 分库工具mysql
  16. solr如何让全词匹配结果在最前面
  17. 2018.4.26 Mac安装Redis5.0.3版本服务器
  18. 求二叉树中第K层结点的个数
  19. PhoneGap安装配置
  20. C# Socket编程 笔记,Socket 详解,入门简单

热门文章

  1. Ajaxterm + nginx 实现一个简单的堡垒机
  2. IE8的 JS 引擎如此不堪(二) - 解决方案
  3. [php] 解析JSON字符串
  4. TDD开发案例
  5. [UOJ198]时空旅行
  6. 【枚举】bzoj1643 [Usaco2007 Oct]Bessie&#39;s Secret Pasture 贝茜的秘密草坪
  7. 【欧拉函数表】POJ2478-Farey Sequence
  8. Problem A: 零起点学算法16——鸡兔同笼
  9. WPF附加属性的Set函数不调用的问题
  10. 用Qemu模拟vexpress-a9 (五) --- u-boot引导kernel,device tree的使用