プログラマ専用SNS ミクプラ

コード:

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
int main(void){
   double a=1103515245;
   double c=12345;
   double Xn=1; //初期値
   double M=1;
   double n=0;
   int i;

  //最大周期 2^31 を求める
   for(i=31;i>0;i--)
   {
      M=M*2;
   }

   while(1){

      Xn=fmod(a*Xn+c,M); //乱数の計算
      n++;

      if(fmod(n,10000000)==0)
         printf("%.0lf\t Xn=%10.0f\n",n,Xn); //10000000ごとに出力
 
      if(Xn==1) break; //1周期終えたら、ループを抜ける

   }

   // 結果出力
   printf("%.0lf\t Xn=%10.0f\n",n,Xn);
   printf("周期　%.0f\n",n);
   return EXIT_SUCCESS;

}

コード:

#include<stdio.h>
#include<stdlib.h>
#include<math.h>

int main(void){
  long long a=1103515245;
  long long c=12345;
  long long Xn=1; //初期値
  long long M=pow(2, 31);
  long long n=0;

  while(1){

    Xn=(a*Xn+c) % M; //乱数の計算
    n++;

    if(n % 10000000 == 0)
      printf("%lld\t Xn=%lld\n",n,Xn); //10000000ごとに出力

    if(Xn==1) break; //1周期終えたら、ループを抜ける

  }

  // 結果出力
  printf("%lld\t Xn=%lld\n",n,Xn);
  printf("周期　%lld\n",n);
  return EXIT_SUCCESS;
}