poj 2462 Period of an Infinite Binary Expansion

原创 2013年12月02日 18:40:06

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

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<vector>
#include<cmath>
#include<set>
#include<map>
#define LL long long
#define CLR(a, b) memset(a, b, sizeof(a))
#define REP(i, n) for(int i = 0; i < n; i ++)
using namespace std;
const int N = 400100;
bool isp[N];
vector<int> p;
vector<LL> hav;

void get_P()
{
    CLR(isp, true);p.clear();
    for(int i = 2; i < N; i ++)
    {
        if(isp[i])
        {
            p.push_back(i);
            if(i < 1111) for(int j = i * i; j < N; j += i)
            {
                isp[j] = false;
            }
        }
    }
}

LL Euler_phi(LL n)
{
    LL ret = n;
    for(int i = 0; (LL)p[i] * p[i] <= n; i ++) if(n % (LL)p[i] == 0)
    {
        ret = ret / p[i] * (p[i] - 1);
        while(n % p[i] == 0) n /= p[i];
    }
    if(n > 1) ret = ret / n * (n - 1);
    return ret;
}

LL Mul(LL a, LL b, LL mod)
{
    LL ret = 0;
    while(b)
    {
        if(b & 1)
            ret = (ret + a) % mod;
        a = a * 2 % mod;
        b >>= 1;
    }
    return ret;
}

LL Pow(LL a, LL b, LL mod)
{
    LL ret = 1;
    while(b)
    {
        if(b & 1) ret = Mul(ret, a, mod);
        a = Mul(a, a, mod);
        b >>= 1;
    }
    return ret;
}

LL gcd(LL a, LL b)
{
    return b ? gcd(b, a % b) : a;
}

void get_hav(LL n)
{
    hav.clear();
    for(int i = 0; i < p.size() && n > 1; i ++)
    {
        while(n % (LL)p[i] == 0)
        {
            n /= p[i];
            hav.push_back(p[i]);
        }
    }
    if(n > 1) hav.push_back(n);
}

int main()
{
    int cas = 1;
    LL ans, m, x, a, b, g;get_P();
    while(scanf("%I64d/%I64d", &a, &b) != EOF)
    {
        g = gcd(a, b);
        a /= g;b /= g;ans = 1;
        while(b % 2 == 0)
        {
            ans ++;
            b /= 2;
            a %= b;
            g = gcd(a, b);
            a /= g;b /= g;
        }
        x = Euler_phi(b);
        get_hav(x);
        for(int i = 0; i < hav.size(); i ++)
        {
            if(Pow(2LL, x / hav[i], b) == 1)
                x /= hav[i];
        }
        printf("Case #%d: %I64d,%I64d\n", cas ++, ans, x);
    }
}


版权声明:本文为博主原创文章,未经博主允许不得转载。转载请标注:blog.csdn.net/ok_again

POJ 1852 Ants(弹性碰撞问题)

Ants Time Limit: 1000MS Memory Limit: 30000K Total Submissions: 12607 Accepted:...
  • zwj1452267376
  • zwj1452267376
  • 2015年11月13日 17:24
  • 1084

Find Bug 使用说明

说明: 1.Bug是findbug Eclipse插件原生的bug信息描述,Confidence 是fingbug团队认为该代码导致bug的可能性。 2.以下都是我使用findbug在公司项目中找...
  • cai514667671
  • cai514667671
  • 2015年10月21日 13:43
  • 1723

Objective-C 警告相关说明

原文地址:http://nshipster.com/clang-diagnostics/
  • icmmed
  • icmmed
  • 2014年09月03日 16:36
  • 9920

POJ3358 Period of an Infinite Binary Expansion

/* 《数论及应用》P168 */ #pragma warning(disable:4786) #pragma comment(linker, "/STACK:102400000,102400...
  • jijijix
  • jijijix
  • 2016年08月28日 10:59
  • 173

POJ 3358 Period of an Infinite Binary Expansion 解题报告(欧拉函数+因式分解)

Period of an Infinite Binary Expansion Time Limit: 1000MS   Memory Limit: 65536K Tot...
  • kbdwo
  • kbdwo
  • 2014年04月22日 17:05
  • 667

POJ 3358 Period of an Infinite Binary Expansion 欧拉函数+欧拉定理

题意:输入一个有理数p/q(保证是一个小数),然后将其小数部分用二进制表示。求出在此种表示下的循环起点和循环节长度 {x} = 0.a1a2...ar(ar+1ar+2...ar+s)w 题解:来...
  • Tsaid
  • Tsaid
  • 2012年02月28日 23:36
  • 1191

POJ_3358 Period of an Infinite Binary Expansion

http://poj.org/problem?id=3358 题意: 给你一个有理分数,要你求这个分数写成2进制小数时的最小循环节和最开始的位置 思路: 一开始以为可以直接暴力模拟搞,因为没有...
  • ivan_zjj
  • ivan_zjj
  • 2012年05月23日 16:13
  • 900

poj_3358_Period of an Infinite Binary Expansion

Period of an Infinite Binary Expansion Time Limit: 1000MS   Memory Limit: 65536K Tot...
  • u010770930
  • u010770930
  • 2013年08月20日 15:54
  • 439

POJ-3358-Period of an Infinite Binary Expansion 解题报告

欧拉定理,同余运算性质,好题。题意:给你两个整数p和q,请输出p/q作为小数在二进制表示下的第一个循环节在第几位小数和最小循环节长度。比如1/10的二进制小数表示为0.0001100110011001...
  • JZQT_T
  • JZQT_T
  • 2015年03月09日 20:14
  • 480

HDU 3358 Period of an Infinite Binary Expansion 推荐!!(欧拉函数,费马小定理)

转载请注明出处,谢谢 http://blog.csdn.net/ACM_cxlove?viewmode=contents           by---cxlove 题目:http://poj....
  • ACM_cxlove
  • ACM_cxlove
  • 2012年08月05日 16:47
  • 3438
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:poj 2462 Period of an Infinite Binary Expansion
举报原因:
原因补充:

(最多只允许输入30个字)