这道题目感觉好难,根本就是无从下手的感觉,尝试了以前的所有方法,都没有思路,毫无进展,参考了一下别人的思路,感觉学到了新的知识

接下来开始分析

观察1/10这组数据,按照二进制转化法可以得到: 1/10 2/104/108/1016/1032/10.……

对于每一个分子进行模10处理 可以相应的得到:    1/102/104/108/106/102/10……

出现了重复,这个重复就是要求的最小循环

对于p/q,首先p'=p/gcd(p,q),q'=q-gcd(p,q),然后求p'*2^i ≡ p'*2^j  (mod q'),然后开始变换,p'*2^i*(2^(j-i)-1) ≡ 0 (mod q'),也就是说 q'|p'*2^i*(2^(j-i)-1),因为 gcd(p',q')=1所以q'|2^i*(2^(j-i)-1)

因为2^(j-i)-1肯定为奇数,所以q'有多少个2的幂,i就是多少,而且i就是循环开始前的第一位数字,令q''为q'除去2的幂之后的数,此时q''|2(j-i)-1,实际上就是 求 某个x  使得  2^x ≡ 1(mod q'');因为q''与2是互诉的,所以肯定有解,令 n=q'',   2^φ(n) ≡ 1 (mod n ),由于题目要求的是 最小的x,看似 φ(n) 是最终解,所以不妨 像poj3696那样大胆假设 x其实是 φ(n)的一个因子,推导符合题目要求,再反过来假设x不是φ(n)的因子,  令r=φ(n),mod x,r大于0,同时r<x,注意 2^φ(n)  ≡  1(mod n),且  2^x≡  1(mod n),所以2^r %n=1,那么就存在一个比x更小的正整数 是的 2^r ≡ 1(mod n),所以  第二个假设失败,所以 x为 φ(n)的因子, 不断的寻找φ(n)的因子 然后判断是否符合题目要求即可

#include<iostream>
#include<cstdio>
#include<list>
#include<algorithm>
#include<cstring>
#include<string>
#include<queue>
#include<stack>
#include<map>
#include<vector>
#include<cmath>
#include<memory.h>
#include<set> #define ll long long #define eps 1e-8 #define inf 0xfffffff
const ll INF = 1ll<<61; using namespace std; //vector<pair<int,int> > G;
//typedef pair<int,int > P;
//vector<pair<int,int> > ::iterator iter;
//
//map<ll,int >mp;
//map<ll,int >::iterator p;
// int a,b,r[112]; int Gcd(int a,int b)
{
return b==0?a:Gcd(b,a%b);
} void detal(int x,int m)
{
r[0]=x%m;
for(int i=1;i<32;i++)
r[i]=(ll)r[i-1]*r[i-1]%m;
} int quick(int x,int y,int m)
{
int i=0,ans,j;
int d[50];
while(y)
{
d[i++]=y%2;
y>>=1;
}
for(j=0,ans=1;j<i;++j)
if(d[j])
ans=(ll)ans*r[j]%m;
return ans;
} int main()
{
int Case=0;
char s[102];
while(~scanf("%s",s))
{
a=0,b=0;
bool flag=false; for(int i=0;i<strlen(s);i++)
{
if(s[i] != '/' && !flag)
a=a*10+s[i]-'0';
else
flag=true;
if(flag && s[i] != '/')
b=b*10+s[i]-'0';
} if(a == 0)
{
printf("Case #%d: %d,%d\n",++Case,1,1);
continue;
}
int gcd=Gcd(a,b);
b/=gcd;
int x=0;
while(!(b&1))
{
b>>=1;
x++;
}
x++;
int c[112][2],k=0,a=b,ans=b;
for(int i=2;i*i<=a;i++)
{
if(a%i == 0)
{
ans-=ans/i;
a/=i;
while(a%i == 0)
a/=i;
}
}
if(a > 1)
ans-=ans/a;;
a=ans;
for(int i=2;i*i<=a;i++)
{
if(a%i == 0)
{
c[k][0]=i;
c[k][1]=0;
c[k][1]++;
while(a%i == 0)
{
c[k][1]++;
a/=i;
}
k++;
}
}
if(a > 1)
{
c[k][0]=a;
c[k][1]=1;
k++;
}
detal(2,b);
for(int i=0;i<k;i++)
{
for(int j=0;j<c[i][1];j++)
{
if(quick(2,ans/c[i][0],b) != 1)
break;
ans/=c[i][0];
}
}
printf("Case #%d: %d,%d\n",++Case,x,ans);
}
return EXIT_SUCCESS;
}

poj3358 Period of an Infinite Binary Expansion 数论有难度的更多相关文章

  1. poj3358 Period of an Infinite Binary Expansion

    Period of an Infinite Binary Expansion 题目大意:给你一个分数,求这个分数二进制表示下从第几位开始循环,并求出最小循环节长度. 注释:int范围内. 想法:这题说 ...

  2. poj 3358 Period of an Infinite Binary Expansion

    由乘2取整得到分数的小数位,可以找到规律!!! 例如:1/10,2/10,4/10,8/10,16/10,32/10,64/10…… 取整后:1/10,2/10,4/10,8/10,6/10,2/10 ...

  3. poj 2462 Period of an Infinite Binary Expansion

    欧拉定理.根据分数转换成2进制的过程,分子每次都乘2.对于循环节x,当2^x = 1(mod b)时肯定是循环节.显然当分母不能整除2的时候,即分母和2互质的话,就可以利用欧拉定理,使得2^(Eule ...

  4. Period of an Infinite Binary Expansion 题解

    Solution 简单写一下思考过程,比较水的数论题 第一个答案几乎已经是可以背下来的,在此不再赘述 考虑我们已经知道了\((p,q)\),其中\((p \perp q) \wedge (q \per ...

  5. sgu 137. Funny Strings 线性同余,数论,构造 难度:3

    137. Funny Strings time limit per test: 0.25 sec. memory limit per test: 4096 KB Let's consider a st ...

  6. ACM数学

     1.burnside定理,polya计数法 这个专题我单独写了个小结,大家可以简单参考一下:polya 计数法,burnside定理小结 2.置换,置换的运算 置换的概念还是比较好理解的,< ...

  7. acm数学(转)

    这个东西先放在这吧.做过的以后会用#号标示出来 1.burnside定理,polya计数法    这个大家可以看brudildi的<组合数学>,那本书的这一章写的很详细也很容易理解.最好能 ...

  8. [转] POJ数学问题

    转自:http://blog.sina.com.cn/s/blog_6635898a0100magq.html 1.burnside定理,polya计数法 这个大家可以看brudildi的<组合 ...

  9. ACM数学知识体系

    在盛情收到学弟邀请给他们整理ACM数学方面的知识体系,作为学长非常认真的弄了好久,希望各学弟不辜负学长厚爱!!!非常抱歉因为电脑全盘格式化好多word.PPT都丢失,我尽量具体地给大家找到各知识点学习 ...

随机推荐

  1. Netfilter-packet-flow.svg

    调试网络的方法:(Debugging the kernel using Ftrace)  $ watch -n1 -d sudo cat /proc/net/snmp$ watch -n1 -d su ...

  2. C与C++不同

    常量表示方法不同 C不支持引用,C++支持 注释不同,C89不支持单行注释 (++i)++在C中不合法 (a=3)=4在C中不合法 不能在for循环头部定义变量 C++注重类型,强类型,严格检查类型 ...

  3. 魔兽世界---屠夫(Just a Hook)

    Just a Hook Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Tota ...

  4. JAVA基础 (二)反射 深入解析反射机制

    在谈论到反射这个问题时,你是否有例如以下疑问? 不管是在.NET还是Java中反射的原理和机制是一样的,理解了一种还有一种就能够迎刃而解,想要理解反射首先须要了解底层的一些概念和执行.理解了反射有助于 ...

  5. 最终有SpringMvc与Struts2的对照啦

    眼下企业中使用SpringMvc的比例已经远远超过Struts2,那么两者究竟有什么差别,是非常多刚開始学习的人比較关注的问题,以下我们就来对SpringMvc和Struts2进行各方面的比較: 1. ...

  6. JavaScript之怎样获取元素节点

    JavaScript获取元素节点一共有三种方法,分别是通过元素ID.通过标签名字和通过类名字来获取: 1.通过元素ID属性的ID值来获得元素对象-getElementById() DOM提供了一个名为 ...

  7. CentOS安装常用软件

    下载第三方库rpmforge,找到合适自己版本的rpmforge下载,用以支持NTFS格式硬盘和MP3格式音频或其他 http://pkgs.repoforge.org/rpmforge-releas ...

  8. 获取 web容器中的bean

    public class WebContextBeanFinder { public static Object getBean(String beanId) { ServletContext ser ...

  9. 浅谈Mybatis(二)

    一.resultMap 作用:发现数据库的查询结果与实体之间不匹配时,需要通过ResultMap来进行映射处理.常用于多表查询. 多表查询还是比较复杂的,因为可能的情况很多.这里只说两种情况: 1.1 ...

  10. <jsp:directive.page>标签

    directive 英 [dɪ'rektɪv; daɪ-] 美 [daɪ'rɛktɪv] n. 指示:指令 adj. 指导的:管理的 等效于 <%page import="com.ct ...