随机数频数_游程检验-安全算法组件生成方法

Class1.cs中的代码

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Threading.Tasks;

namespace ClassLibrary1
{
    public class Myclass
    {
        public double cephes_erfc(double x)
        {
            const double one_sqrtpi = 0.564189583547756287, rel_error = 1e-12;
            double a = 1, b = x, c = x, d = x * x + 0.5;
            double q1, q2 = b / d, n = 1.0, t;

            if (Math.Abs(x) < 2.2)
                return 1.0 - cephes_erf(x);
            if (x < 0)
                return 2.0 - cephes_erfc(-x);
            do
            {
                t = a * n + b * x;
                a = b;
                b = t;
                t = c * n + d * x;
                c = d;
                d = t;
                n += 0.5;
                q1 = q2;
                q2 = b / d;
            } while (Math.Abs(q1 - q2) / q2 > rel_error);
            return one_sqrtpi * Math.Exp(-x * x) * q2;
        }
        public double cephes_erf(double x2)
        {
            const double two_sqrtpi = 1.128379167095512574, rel_error2 = 1E-12;
            double sum = x2, term = x2, xsqr = x2 * x2;
            int j = 1;

            if (Math.Abs(x2) > 2.2)
                return 1.0 - cephes_erf(x2);
            do
            {
                term *= xsqr / j;
                sum -= term / (2 * j + 1);
                j++;
                term *= xsqr / j;
                sum += term / (2 * j + 1);
                j++;
            } while (Math.Abs(term) / sum > rel_error2);

            return two_sqrtpi * sum;
        }
        public double Caculater_frequency(string str)
        {
            int[] E = new int[200];
            int i = 0;
            int n = str.Length;
            do
            {
                E[i] = Convert.ToInt32(str[i]) - 48;
                i++;
            } while (i < n);
            double  S_n = 0;
            for (i = 0; i < n; i++)
                S_n = (2 * E[i] - 1) + S_n;
            if (S_n < 0)
                S_n = -S_n;
            double S_obs = S_n / Math.Sqrt(n);
            return cephes_erfc(S_obs / Math.Sqrt(2));

        }
        public double Caculater_runs(string str)
        {
            int[] E = new int[200];
            int i = 0;
            int n = str.Length;
            double V_obs = 0;
            int r_i;
            do
            {
                E[i] = Convert.ToInt32(str[i]) - 48;
                i++;
            } while (i < n);
            double  S_n = 0;
            for (i = 0; i < n; i++)
                S_n = E[i] + S_n;
            double t = S_n / n;
            if (t == 0 || t == 1)
                return 0 ;
            else if (Math.Abs(t - 0.5) >= (2 / Math.Sqrt(n)))
                return 0 ;
            else
                for (i = 0; i < n - 1; i++)
                {
                    if (E[i] == E[i + 1])
                        r_i = 0;
                    else
                        r_i = 1;
                    V_obs = r_i + V_obs;
                }
            V_obs = V_obs + 1;
            return cephes_erfc(Math.Abs(V_obs - 2 * n * t * (1 - t)) / (Math.Sqrt(2 * n)) * 2 * t * (1 - t));

        }
    }
}