poj 1597 Uniform Generator【生成指定范围内所有随机数】

本文参考资料:http://hi.baidu.com/bnjyjncwbdbjnzr/item/1f997cfdd225d5d143c36a58

题意:一个生成随机数的函数,

Seed[x+1] = ( seed[x] + STEP ) % MOD

输入step和mod,问能否生成0~MOD-1之间所有的数,是Good Choice,否则Bad Choice

题意其实就是:给出S和M,求0*S%M,1*S%M,2*S%M......(M-1)*S%M能否组成一个集合包含0.1.。。。M-1;(这个是原题意改造而来);

算法:判断两个数是否互质;or 暴力解决

其实暴力完全可以解决这个问题(⊙﹏⊙b),只是其中用数学方法更加高效,巧妙;

下面证明“如果SM互质则满足题意”:

设G=gcd(S,M);则S=A*G,M=B*G;

另设X=K*S%M=K*S-T*M(T为整数,满足X属于0到M-1。其实就是说T是K*S除以M的商。);

X=K*A*G-T*B*G;因此取余后的整数一定是G的倍数。

所以,上述要求证明的命题其实就是要证明:G只能取1才能满足条件;

(个人感觉上述论证就是废话呵呵——S和M互质不就是gcd(S,M)==1吗?)

(下面继续刚才的证明)

充分性的证明:(即当S与M互质,则0到M-1的S倍对M取余一定能遍历0到M-1)

只需证明的是,该余数中两两之间互不相等;

假设k*S和b*S对M取余相等(k和b∈[0,M),并且k和b不等);

则k*S=q1*M+r≡q2*M+r=b*S     <==>      (k-b)*S=M*(q1-q2);

因为题设S与M互质,再由上式子可得M|(k-b),这与k和b∈[0,M),并且k和b不等矛盾;

因此题目命题得证。

另外,偶然看到一个很牛叉的辗转相除法;

int gcd(int a,int b)

{

while(b) b^=a^=b^=a%=b;

return a;

}

此代码,很好很强大;把涉及位运算的交换的程序加入,便到得这段简洁高效的代码;

注:A和B;经过A^=B^=A^=B,结果就得到A和B的交换;

 #include<stdio.h>
long gcd(long a,long b);
int main()
{
long step,mod;
while(scanf("%ld%ld",&step,&mod)!=EOF)
{
if(gcd(step,mod)==) printf("%10d%10d Good Choice\n\n",step,mod);
else printf("%10d%10d Bad Choice\n\n",step,mod);
}
return ;
}
long gcd(long a,long b)
{
long c;
if(b==) return a;
c=a%b;
while(c)
{
a=b;
b=c;
c=a%b;
}
return b;
}

下面是其他网友对该问题的论述:(感觉没有上面的讲的清楚,可以直接忽略)

以数论的观点解释 poj1597 && 更相减损术

http://blog.163.com/jiazheng2222@126/blog/static/1696323832009513104450808/

可以想象成有一个人从起点0处绕着一个圈长为mod个单位“跳步”,每跳一次必须跳step个单位(也就是可不能出现在i*step<pos<(i+1)*step的位置),而至少跳mod次恰好回到起点的话,就是Good Choice,如果在第mod次前就恰好到达起点则为Bad Choice

而我们知道,回到起点把所需的最短时间即是step和mod的最小公倍数,而若为Good Choice,显然最短时间又等于step*mod,由数论的知识gcd(step,mod)=step*mod/(step和mod的最小公倍数)=1

bool check(long step,long mod)
{
    long m=step>mod?step:mod;
    long n=step>mod?mod:step;
    long q=1;

while(q!=0)
    {
        q=m-n;
        m=n>q?n:q;
        n=n>q?q:n;

}
    if(m==1)
        return true;
    else
        return false;
}

http://hi.baidu.com/sjzezoi/item/6a3a71ccab9fee32b67a24fa

题意:一个生成随机数的函数,

Seed[x+1] = ( seed[x] + STEP ) % MOD

问能不能生成0~MOD-1之间所有的数,是Good Choice,否则Bad Choice

此题可以暴力模拟,但正解应该是gcd(STEP,MOD)==1

设G=gcd(a,b),于是有a=A*G,b=B*G(1<=A,B,gcd(A,B)=1) 对于a的多次加可以看成K*a(1<=k),转化成(K*a)%b的所有结果能否表示成0..b-1中的所有数,假(K*a)%b=M,M=K*a-W*b(W为使M>0的最大整数),M=K*A*G-W*B*G M%G==0,既结果是G的倍数,如果想取得0..b-1中的所有数,那么必须G=1才可能..

(代码仅供参考和学习,请不要直接粘贴刷AC数,期待你写出更好的代码)贡献一次PE

暴力:

#include<cstdio>

#include<cstring>

int a,b;

bool vis[100010];

int main()

{

while(scanf("%d%d",&a,&b)+1)

{

memset(vis,0,sizeof vis);

int cnt=0,now=0;

while(cnt<b)

{

if(vis[now]) break;

vis[now]=1;

now=(now+a)%b;

cnt++;

}

if(cnt==b) printf("%10d%10d    Good Choice\n\n",a,b);

else printf("%10d%10d    Bad Choice\n\n",a,b);

}

return 0;

}

正解:

#include<cstdio>

int a,b;

inline int gcd(int x,int y)

{

if(!y) return x;

return gcd(y,x%y);

}

int main()

{

while(scanf("%d%d",&a,&b)==1)

{

if(gcd(a,b)==1) printf("%10d%10d    Good Choice\n\n",a,b);

else printf("%10d%10d    Bad Choice\n\n",a,b);

}

return 0;

}

http://hi.baidu.com/zhuangxie1013/item/cc5a130b05e77ad573e676b0

这题其实是道数论题,但我把它当作模拟题来做了,虽然做出来了,但所用的空间和时间也大的吓人,此题输出时需要注意:输出格式为printf("%10d%10d    Bad
Choice\n\n",step,mod);即第二个数和Bad Choice之间有4个空格,下面是我模拟出来的程序:

#include<iostream>
#include<set>
#include<algorithm>
#include<cstdlib>
using namespace std;

int main()
{
int step,mod,i,temp;
set<int>s;
while(scanf("%d %d",&step,&mod)!=EOF)
{
   temp=0;
   s.insert(temp);
   for(i=1;i<mod;i++)
   {
    temp=(temp+step)%mod;
    s.insert(temp);
   
if(s.size()!=(i+1))                                                                
//出现重复元素
    {
     printf("%10d%10d    Bad
Choice\n\n",step,mod);
     break;
    }

}
   temp=(temp+step)%mod;
   if(temp==0&&s.size()==mod)
    printf("%10d%10d    Good
Choice\n\n",step,mod);
     s.clear();
}
system("pause");
return 0;
}

看了别人的方法,讲的都是数论,所以将其转载:

题目的意思就是一个生成随机数的函数,

Seed[x+1] = ( seed[x]
+ STEP ) % MOD

其中seed就是我们生成出来的随机数,至于seed[0]是哪个数并不重要,后面会证明。STEP就是我们每次往前一个所加的值,再module上MOD得到下一个随机数。

判断这个随机生成函数的好坏的依据是如果能够产生0~MOD-1内的所有数,就是一个好的,否则坏。

我们知道了同余的特性,便可以假设在k步之后,生成的seed[k] = seed[0],所以由此有

Seed[k] = ( seed[0]
+ STEP*k ) % MOD

那么,

STEP * k = MOD

而我们如果要生产MOD个数,必须使k≥MOD,而且k不可能大于MOD,因为这个条件下生成的数又开始重复,所以k=MOD;在前面的条件下,如果STEP和MOD有大于1的公约数例如d,那么会有STEP*(k/d) = MOD,这就是一个循环了,只会产生k/d<MOD个随机数。

结论:iff gcd(STEP, MOD) == 1, good choice

上一篇:PC/UVa 题号: 110106/10033 Interpreter (解释器)题解 c语言版


下一篇:C语言学习第六章