Cyclic Number 资料

最新推荐文章于 2022-04-28 09:46:07 发布

weixin_30757793

最新推荐文章于 2022-04-28 09:46:07 发布

阅读量172

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/yefeng1627/archive/2013/05/14/3077399.html

本文深入探讨了Cyclic Number的概念及其与单位分数之间的联系，详细阐述了如何通过伪代码构造Cyclic Number，并指出在不同进制下的特殊性质。重点介绍了通过计算b进制下的1/p (p为素数)来生成Cyclic Number的方法，同时讨论了完全平方数进制下不存在长度超过1的Cyclic Number的情况。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

全文资料来自wiki:http://en.wikipedia.org/wiki/Cyclic_number

Cyclic number

Cyclic Number 是一个整数是，且其连续的倍数 (1,2,..,L-1) ，是其本身数字的一个循环。常见的如下

142857 × 1 = 142857
142857 × 2 = 285714
142857 × 3 = 428571
142857 × 4 = 571428
142857 × 5 = 714285
142857 × 6 = 857142

而形如下非连续的倍数不是

076923 × 1 = 076923
076923 × 3 = 230769
076923 × 4 = 307692
076923 × 9 = 692307
076923 × 10 = 769230
076923 × 12 = 923076

若不允许前导0，则十进制下，仅有一个Cyclic Number, 142857

Cyclic Number的形式

　　从上面单位分数形式，得到如下形式： $\frac{b^{p-1}-1}{p}$ 其中b是进制，p是素数，且 b%p != 0 .

　　不是所有的素数p都能通过以上形式产生 Cyclic Number, 例如 p=13 gives 076923076923，会出现多个循环。

构造Cyclic Number

伪代码如下：

Let b be the number base (10 for decimal)
Let p be a prime that does not divide b.
Let t = 0.
Let r = 1.
Let n = 0.
loop:
Let t = t + 1
Let x = r · b
Let d = int(x / p)
Let r = x mod p
Let n = n · b + d
If r ≠ 1 then repeat the loop.
if t = p − 1 then n is a cyclic number.

此程序是通过计算 b进制下的Cyclic Number 1/p . 其中r是每一步的剩余系，d是当前步的产生数，当前值为 n = n*b + d.

注意，若 t > p/2, 则这个数必定为 Cyclic Number, 无需继续计算后面部分了。

当 b 为完全平方数时，不存在长度超过1的Cyclic Number.

用以上方法写的代码：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;

typedef long long LL;

bool IsPrime(int x){
    for(int i = 2; i*i <= x; i++)
        if(x%i == 0) return false;
    return true;
}
int main(){
    int L, x, p;
    scanf("%d%d", &L, &x);
    p = L+1;
    if( (x == 2) || !IsPrime(p) ){
        puts("-1"); return 0;        
    }
    for(int b = x-1; b > 1; b--){
        if( b%p == 0 ) continue;
        LL t = 0, r = 1, n = 0, x, d;
        while(1){
            t = t+1; x = r*b; d = x/p; r = x%p; 
            n = n*b + d;
            if( r == 1 ) break;
            if( t > p/2 ){ printf("%d\n",b); return 0; }    
        }
        if( t == L ){ printf("%d\n",b);return 0; }    
    }        
    puts("-1");    
    return 0;
}