hdu 2462 The Luckiest number

原创 2013年12月02日 17:44:05

            欧拉定理。这道数论题好凶残啊!

           首先,由题意可以得出,(10^x - 1)/ 9 * 8 = L * p(p是一个未知数,但必定是整数)。

           然后对上式进行移项处理,得:(10^x - 1) = 9 * L * p / 8。

           设m = 9 * L / gcd(L, 8),则有(10^x - 1) = m * p'。p’是必然存在的一个整数。

           然后问题就转化成为了 10^x = 1(mod m),观察此式,显然,m和10必定互质。

           于是根据欧拉定理,10^(Euler(m)) = 1(mod m) 。由于题目要求最小的解,解必然是Euler(m)的因子。

           需要注意的是,对于10^x,由于m太大,直接快速幂相乘的时候会超long long。。。。好bug,需要乘法转化一下。

#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 L, ans, m, x;get_P();
    while(scanf("%I64d", &L), L)
    {
        m = 9 * L / gcd(L, 8LL);
        if(gcd(m, 10LL) != 1)
        {
            printf("Case %d: 0\n", cas ++);
            continue;
        }
        x = Euler_phi(m);
        get_hav(x);
        for(int i = 0; i < hav.size(); i ++)
        {
            if(Pow(10LL, x / hav[i], m) == 1)
                x /= hav[i];
        }
        printf("Case %d: %I64d\n", cas ++, x);
    }
}


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

hdu3565 Bi-peak Number (有上界和下界的数位dp)

Problem Description A peak number is defined as continuous digits {D0, D1 … Dn-1} (D0 > 0 and n >...
  • Kirito_Acmer
  • Kirito_Acmer
  • 2016年04月06日 21:56
  • 488

HDU 5447 Good Numbers (2015年长春赛区网络赛K题)

1.题目描述:点击打开链接 2.解题思路:本题利用唯一分解定理+gcd解决。这是我第一次用Java写题,忙了一下午终于通过了,真是太不容易了==。本题实际上就是把k1,k2分解后,把所有素因子的指数...
  • u014800748
  • u014800748
  • 2015年09月17日 19:47
  • 1195

【BeiJing2011】【BZOJ2462】矩阵模板

Description给定一个M行N列的01矩阵,以及Q个A行B列的01矩阵,你需要求出这Q个矩阵哪些在 原矩阵中出现过。 所谓01矩阵,就是矩阵中所有元素不是0就是1。 Input输入文件...
  • CreationAugust
  • CreationAugust
  • 2015年10月15日 21:06
  • 1839

HDU - 6231 K-th Number (2017CCPC哈尔滨站 二分+尺取法)

K-th Number Problem Description Alice are given an array A[1..N] with N numbers. Now Al...
  • lzc504603913
  • lzc504603913
  • 2017年11月12日 20:26
  • 299

HDU 1238 Substrings (最长公共子串+DFS)

Substrings Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total...
  • liangzhaoyang1
  • liangzhaoyang1
  • 2016年06月03日 13:04
  • 490

HDU3308:LCIS(线段树区间合并)

Problem Description Given n integers. You have two operations: U A B: replace the Ath number by B...
  • libin56842
  • libin56842
  • 2013年11月10日 15:48
  • 3495

【hdu】1698 just a hook 中文题意&题解

Just a Hook Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Prob...
  • WilliamCode
  • WilliamCode
  • 2016年04月12日 09:09
  • 673

ACM-字符串处理之Substrings——hdu1238

ACM 字符串处理 Substrings hdu1238 字符串函数
  • lx417147512
  • lx417147512
  • 2014年05月01日 13:11
  • 1701

HDU 4279 Number(找规律)

原题:http://acm.hdu.edu.cn/showproblem.php?pid=4279 Number Time Limit: 2000/1000 MS (Java/Others) ...
  • YJX_xx
  • YJX_xx
  • 2014年05月08日 23:09
  • 955

【HDU】2665Kth number(划分树)

只是为了存模板~ #include #include #include using namespace std; const int maxn = 100005; int n,m; int sort...
  • u013451221
  • u013451221
  • 2015年04月10日 17:52
  • 478
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:hdu 2462 The Luckiest number
举报原因:
原因补充:

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