栏目分类:
子分类:
返回
名师互学网用户登录
快速导航关闭
当前搜索
当前分类
子分类
实用工具
热门搜索
名师互学网 > IT > 面试经验 > 面试问答

poj 3358 Period of an Infinit...

面试问答 更新时间: 发布时间: IT归档 最新发布 模块sitemap 名妆网 法律咨询 聚返吧 英语巴士网 伯小乐 网商动力

poj 3358 Period of an Infinit...

#include <iostream>#include <cstdio>#include <cstring>#include <algorithm>#include <cmath>using namespace std;int t, n, m, GCD, phi, ans1, ans2;int temp, num;int fac[1000000];int gcd(int x, int y){    int r;    if(x < y){        swap(x, y);    }    while(x%y){        r = x%y;        x = y;        y = r;    }    return y;}int euler(int x){    int y = (int)sqrt(x+0.5);    int curans = x;    for(int i = 2;i <= y;i++){        if(x%i == 0){ curans = curans/i*(i-1); while(x%i == 0){     x /= i; }        }    }    if(x > 1){        curans = curans/x*(x-1);    }    return curans;}int quickpow(int m, int n, int k){    int b = 1;    while(n > 0){        if(n&1){ b = (long long)b*m%k;        }        n >>= 1;        m = (long long)m*m%k;    }    return b;}int main(){    int T = 0;    while(scanf("%d/%d", &n, &m)!=EOF){        T++;        GCD = gcd(n, m);        n /= GCD;        m /= GCD;        t = 0;        while(m%2 == 0){ t++; m /= 2;        }        ans1 = t+1;        phi = euler(m);        if(phi == 1){ ans2 == 1;        }        else{ num = 0; for(int i = 1;i*i <= phi;i++){     if(phi%i == 0)     {         fac[num++] = i;         fac[num++] = phi/i;     } } sort(fac, fac+num); for(int i = 0;i < num;i++){     temp = quickpow(2, fac[i], m);     if(temp == 1){         ans2 = fac[i];         break;     } }        }        printf("Case #%d: ", T);        printf("%d,%dn", ans1, ans2);    }    return 0;}
转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/372338.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

版权所有 (c)2021-2022 MSHXW.COM

ICP备案号:晋ICP备2021003244-6号