<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>The Testing of Pseudorandom Sequences using Multidimensional Statistics</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Taras Shevchenko National University of Kyiv</string-name>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Bohdana Havrylyshyna str.</string-name>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Ukraine spopereshnyak@gmail.com</string-name>
        </contrib>
        <contrib contrib-type="author">
          <string-name>University of Library Studies</string-name>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Information Technologies</string-name>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Tsarigradsko Shose</string-name>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Sofia</string-name>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Bulgaria geo.p.dimitrov@gmail.com</string-name>
        </contrib>
      </contrib-group>
      <abstract>
        <p>The available approaches to testing pseudorandom sequences show low flexibility and versatility in the means of finding hidden patterns in the data. To solve this problem, it is suggested to use algorithms based on multidimensional statistics. The paper proposed a new approach for testing pseudorandom sequences, obtained an explicit form of the joint distribution of numbers of 2-chains and numbers of 3-chains of various options random bit sequence of a given small length. Examples, tables, diagrams that can be used to test for randomness of the location of zeros and ones in the bit section are presented. In future as a result an information system will be created that will allow analyzing the pseudorandom sequence of a small length and choosing a quality pseudorandom sequence for use in a particular subject area.</p>
      </abstract>
      <kwd-group>
        <kwd>Algorithms</kwd>
        <kwd>multidimensional Statistics</kwd>
        <kwd>Random Sequence</kwd>
        <kwd>schains</kwd>
        <kwd>Cryptography</kwd>
        <kwd>Pseudorandom Sequence</kwd>
        <kwd>Statistical Testing</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>Random sequences have found the widest application from the gaming computer
industry to mathematical modeling and cryptology.</p>
      <p>We list some areas of their usage: modeling, cryptography and information
security, decision making in automated expert systems, optimization of functional
dependencies, fun and games.</p>
      <p>There are various approaches to the formal definition of the term “randomness”
based on the concepts of computability and algorithmic complexity [1-2].</p>
      <p>By implementing some algorithm, software generators produce numbers (although
not obvious) depending on the set of previous values, so the received numerical
sequences are not truly random and are called pseudo-random sequences (PRS). At the
moment, more than a thousand software PRS generators are known, which differ in
algorithms and values of parameters. Statistical properties are significantly different
from the number sequences that are generated by them.</p>
      <p>Copyright © 2019 for this paper by its authors. Use permitted under Creative Commons License
Attribution 4.0 International (CC BY 4.0)
2019 DCSMart Workshop.</p>
      <p>The presented and not presented results allow us to characterize the state of modern
technologies of designing the PRS (focusing on the most progressive of them by the
following basic provisions [3-6].
2</p>
    </sec>
    <sec id="sec-2">
      <title>Problem Statement</title>
      <p>Before responsible using in mathematical modeling and cryptology, PRS should be
tested. Unfortunately, for many PRS tests, there are some limitations:
• checked out only one of the probable ones properties that are characterize</p>
      <p>PRS;
• not fix family alternatives;
• do not have theoretical ones ratings power.
• do not give a correct an estimate of chance sequences provided a little
sample.</p>
      <p>Problems small and large samples refer to the main problems that arise in practical
application methods analysis data. Let's be use the next classification samples by
number [2], based on requirements presented in the program criteria:
• very small sample - from 5 to 12,
• small sample - from 13 to 40,
• medium sample - from 41 to 100,
• large sample - from 101 and more.</p>
      <p>The minimum size of the sample limits not so much the algorithm of calculating
the criterion, but the distribution of its statistics. For a row algorithms with too much
small ones numbers sample normal approximation distribution of statistics criterion
will be under question.</p>
      <p>During the research, the localization of the local sections of the bit sequence was
conducted to detect the dependencies in the location of its elements by using the exact
distributions of the corresponding statistics. In the work an explicit form of the joint
distribution of the numbers of 2-chains and numbers of 3-chains of various variants in
a random sequence was obtained. This joint distribution allows more accurate
comparison of the use of one-dimensional statistics, to analyze the bit sequence small
length by chance.
3</p>
      <p>Joint Distribution of number of 2-chains and number of
3chains of a provided type in binary sequence
Consider a sequence of random variables
where   = {0, 1}, i= 1, 2, . . . ,  ,  &gt; 0.</p>
      <p>Subsequences   ,   +1, . . . ,   + −1,
1, 2, . . . ,  −  + 1,  = 1, 2, . . . ,  .</p>
      <p>
        1,  2, . . . ,   ,
(
        <xref ref-type="bibr" rid="ref1">1</xref>
        )
sequences (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) are
called
s-chains,  =
with  1 ,  2, .. . ,   , where   = {0, 1},  = 1, 2, . . . ,  .
      </p>
      <p>
        Denote  ( 1  2 . ..   ) the number of s-chains in the sequence (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) that coincide
3,   {0, 1},  ∗ = 1 −  . Then
      </p>
      <p>
        Theorem. Let sequence (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) consist of n,  &gt; 0 independent identically distributed
0, elsewhere
  ∗ + 2(  −  1 +   ),  = 1,3,  1 = −1,  2 = 1,  3 = 0;
Ζ( ,  ) ≝ {
      </p>
      <p>1, if  =  = 0;
  −−11,</p>
      <p>if  ≥  ≥ 1;
Ρ{ (  ) =  1,  ( ∗ ∗ ∗) =  2} = ∑</p>
      <p>1=0  1  0 ×
{  1

  − 1−1Ζ(  ∗ −   +  1 + 1;   −  1 −  1 − 1)+
 2
  − 1+1Ζ(  ∗ −   +  1 − 1;   −  1 −  2 + 1)+</p>
      <p>2  3− 1Ζ(  ∗ −   +  1;   −  1 −  3)+
 (  −  1 − 1,  2,   ∗)}Ζ(  ;   −  1),
Ρ{ (  ) =  1,  ( ∗  ∗) =  2 } = ∑ 1=0  1  0 ×
{  2
  − 1−2  1+1</p>
      <p>− 1− 2−1Ζ(  ∗;   −  1 − 1) 1(  −  1 − 2)+

  − 1  t∗−− 11Ζ( 1;   −  1 −  2)+ 2   − 1−1   ∗−1
 2  2</p>
      <p>− 1− 2−1   − 1− 2−1 +</p>
      <p>1
 2(  −  1 − 1,  2,   ∗)},
where  1( 1) = {0, elsewhere
1, if  1 ≥ 1,</p>
      <p>,  2( 1,  2,  3) = {0, e1l,siefw h1e=re 2 =  3 = 0,
Ρ{ (  ) =  1,  ( ∗ ∗ ∗)+  ( ∗  ∗) =  2 } = ∑ 1=0  1  0 × {(∑ + ∗= 1 1 ×
 
  − 1−2  ∗− 1−1  1+1</p>
      <p>− 1−δ−1Ζ(  ∗ −   +  1 + 1;   −  1 −  ∗ − 1)+
(∑ + ∗= 2  δ</p>
      <p>∗
  − 1   − 1+1Ζ( 1;   −  1 − δ)Ζ(  ∗ −   +  1 − 1;   −  1 −
 ∗ + 1))+(∑ + ∗= 3  
  − 1−1  ∗− 1  1</p>
      <p>
        − 1−δ−1Ζ(  ∗ −   +  1;   −  1 −  ∗))
+ (  −  1 − 1,  2,   ∗)},
(
        <xref ref-type="bibr" rid="ref2">2</xref>
        )
(
        <xref ref-type="bibr" rid="ref3">3</xref>
        )
(
        <xref ref-type="bibr" rid="ref4">4</xref>
        )
(
        <xref ref-type="bibr" rid="ref5">5</xref>
        )
is the symbol ∑
      </p>
      <p>
        denotes addition over all non-negative integers   and   ∗ such that
As a result of applying this technique for testing pseudo-random sequences for
twodimensional statistics, you can build tables (relations (
        <xref ref-type="bibr" rid="ref2">2</xref>
        ) - (
        <xref ref-type="bibr" rid="ref5">5</xref>
        )) and bubble diagrams
(relations (
        <xref ref-type="bibr" rid="ref3">3</xref>
        ) - (
        <xref ref-type="bibr" rid="ref5">5</xref>
        )) with which you can get the probability of the distribution of zeros
and ones in a given sequences.
      </p>
      <p>As practice shows, the use of ready-made tables for analyzing the sequence of
randomness allows you to get the answer as quickly as possible, in contrast to the
classical testing method.</p>
      <p>
        Consider an example of tables and bubble diagrams for a bit-sequence of small
length. For example, let the length of the bit sequence n, n = 32 for relations (
        <xref ref-type="bibr" rid="ref3">3</xref>
        ) - (
        <xref ref-type="bibr" rid="ref5">5</xref>
        )
and n = 24 for relations (
        <xref ref-type="bibr" rid="ref2">2</xref>
        ).
4.1
      </p>
      <p>
        Illustration of the Use of Equality (
        <xref ref-type="bibr" rid="ref2">2</xref>
        )
In Table 1 and in Fig. 1 shows the use of the relation (
        <xref ref-type="bibr" rid="ref2">2</xref>
        ) for a small sample  ,  =
32, and some values  1 and  2.
      </p>
      <p>In Table 1 the first column contains all possible values  1 and  2, for which
probability is Ρ{ (  ) =  1,  (  ∗ ∗ ∗) =  2} ≥ 0,01. The second column of Table 1 gives
the probabilities (in non-decreasing order)  { (  ) =  1,  (  ∗ ∗ ∗) =  2} for pairs
of numbers ( 1,  2) listed in the first column.</p>
      <p>Each row of the fourth column contains the sum of the accumulated probabilities
before the event is implemented { (  ) =  1,  (  ∗ ∗ ∗) =  2} inclusive where  1
and  2 indicated in the same line in the first column.
4.2</p>
      <p>
        Illustration of the Use of Equality (
        <xref ref-type="bibr" rid="ref4">4</xref>
        )
In Table 2 and in Fig. 2. shows the use of the relation (
        <xref ref-type="bibr" rid="ref4">4</xref>
        ) for a small sample of n, n =
32, and some values of  1 and  2.
32, and some values  1 and  2.
 1
12
11
6
4
6
7
9
4
8
9
5
4
10
11
10
4
11
8
5
 2
11
12
10
6
4
6
4
9
9
7
8
4
4
5
5
8
11
11
10
 1
5
4
5
8
7
4
6
7
4
6
7
5
6
6
0,01018
0,01025
0,01028
0,01165
0,01178
0,01179
0,01209
0,01229
0,01285
0,0129
0,01325
0,014
0,01416
0,01481
0,01521
0,01543
0,01578
0,01706
0,01729
 2
1
4
1
1
1
3
3
1
2
2
2
3
2
1
In Table 4 shows the use of the relation (
        <xref ref-type="bibr" rid="ref2">2</xref>
        ) for a small sample  ,  = 24, and some
values  1,  2 and  3.
      </p>
      <p>In Table 4 in the first, second and third columns are all possible values  1,  2 and
 3, for which probability Ρ{ (  ) =  1,  (  ∗ ∗ ∗) =  2,  (  ∗  ∗) =  3 } ≥ 0,009 ,
and the contents of the fourth and fifth columns are similar to the contents of the third
and fourth columns of the Table 1.
5</p>
      <p>
        Results and Discussion
As a result of applying this technique for testing pseudo-random sequences for
twodimensional statistics (relations (
        <xref ref-type="bibr" rid="ref3">3</xref>
        ) - (
        <xref ref-type="bibr" rid="ref5">5</xref>
        )), you can build a bubble diagram with which
you can get the probability of the distribution of zeros and ones in a given sequence.
      </p>
      <p>Consider examples of bubble diagrams for a bit sequence of small length n, n = 32.
5.1</p>
      <p>
        Graphic Illustration of the Use of Equality (
        <xref ref-type="bibr" rid="ref3">3</xref>
        )
      </p>
      <p>After analyzing Fig. 1 it can be concluded that for the analysis of the sequence of
chains of small and medium length (from 13 to 100 elements), one-dimensional
statistics do not always give the correct result. For example, if we consider the sequence
where the parameter  1 = 8, then we can draw a conclusion with a degree of
probability about 10% of randomness of the sequence with these characteristics, however, if
we pay attention when  1 = 8 and  2 = 5 it can be argued that this sequence is
non6
5
2,71%
2,99%</p>
      <p>2,64%
2,38%
3,31%
3,68%
3,36%</p>
      <p>2,59%
2,25%
3,09%
3,50%
3,32%</p>
      <p>2,71%
2,09%
2,42%
2,40%
2,07%
random, therefore as shown in Fig. 1 we have Ρ{ (  ) =  1,  (  ∗ ∗ ∗) =  2} =
1,45%. What also shows the lack of use of one-dimensional statistics for the analysis
of small and medium bit sequences.</p>
      <p>An approach to testing using n-dimensional statistics allows us to rely on a deeper
justification of the randomness of generated sequences.
5.2</p>
      <p>
        Graphic Illustration of the Use of Equality (
        <xref ref-type="bibr" rid="ref4">4</xref>
        )
In Fig. 2 shows the use of the relation (
        <xref ref-type="bibr" rid="ref4">4</xref>
        ) for a small sample  ,  = 32 , and some
values  1 and  2.
      </p>
      <p>4
5
6
7</p>
      <p>
        Valu8e k1
9
10
11
12
In Fig. 3 shows the use of relation (
        <xref ref-type="bibr" rid="ref4">4</xref>
        ) for a small sample  ,  = 32, and some values
 1 and  2.
      </p>
      <p>Fig. 3 gives a bubble chart in which the first parameter (horizontal axis) is the
value  1, the second parameter (vertical axis) is the value  2, and the third parameter
(bubble size) is the probability of the event occurring { (  ) =  1,  (  ∗ ∗ ∗) +
 (  ∗  ∗) =  2 }, which is represented as a percentage.
6
•
•
•
•
•
2
k
eu 7
l
a
V
6
5
4
2,16%
2,42%
2,61%
2,70%</p>
      <p>2,11%
2,03%
2,78%
2,74%</p>
      <p>2,08%
2,02%
2,62%
2,51%
2,16%</p>
      <p>2,06%
9
10
11
4
5
6
7</p>
      <p>8</p>
      <p>Value k1</p>
      <p>
        In this paper, the exact compatible distributions of some statistics (
        <xref ref-type="bibr" rid="ref1">0, 1</xref>
        ) -sequences
of length 1 &lt;  &lt; ∞ are given. For a bit sequence of small length n, n = 32, the tables
containing the numerical values of the corresponding distribution are given. These
tables, as well as the proposed graphic representations, can be used to test the
hypothesis of the randomness of the arrangement of zeros and units.
      </p>
      <p>The Results of the Comparison the NIST Statistical Test Suite
and Test of PRS of Small Length using Multidimensional
Statistics
Consider the well-known examples that are given in [7, 8]. Let us analyze the
submitted sequences for the corresponding tests, where:</p>
      <p>
        P is the probability of sequence randomness according to the selected criterion
from the first column,
P1 is the probability obtained using relation (
        <xref ref-type="bibr" rid="ref2">2</xref>
        ),
P2 is the probability obtained using relation (
        <xref ref-type="bibr" rid="ref3">3</xref>
        ),
P3 is this is the probability obtained using relation (
        <xref ref-type="bibr" rid="ref4">4</xref>
        ),
P4 is this is the probability obtained using relation (
        <xref ref-type="bibr" rid="ref5">5</xref>
        ).
      </p>
      <p>As can be seen from the table, the use of two-dimensional statics gives a more
accurate result for short sequences. And also, according to [8], the recommended
minimum sequence length n is greater than 100 bits.</p>
    </sec>
    <sec id="sec-3">
      <title>Conclusions</title>
      <p>The available approaches to testing pseudorandom sequences show low flexibility and
versatility in the means of finding hidden patterns in the data. To solve this problem,
it is suggested to use algorithms based on multidimensional statistics.</p>
      <p>The approach to testing using multidimensional statistics allows you to rely on a
deeper justification of the randomness of the generated sequences. This area is
promising for scientific research.</p>
      <p>The paper proposed a methodology for testing a sequence and obtained a correct
view of the joint distribution of the numbers of 2-chains and the numbers of 3-chains
of various variants in a random bit sequence of a given small length.</p>
      <p>These algorithms and scheme of work for verification statistical tests of
randomness sequences (proposed in chapter II) combine all the advantages of statistical
methods and are the only alternative for the analysis of sequences of small and
medium length.</p>
      <p>To implement the proposed approach, a PRS software test package is being
developed, which will include tests using multidimensional statistics, which are well
recommended for testing a small length PRS. As a result of the implementation of this
technique, an information system will be created that will allow analyzing the PRS of
a small length and choosing a quality PRS for use in a particular subject area.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Маsоl</surname>
            <given-names>V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Pоpereshnyаk</surname>
            <given-names>S.</given-names>
          </string-name>
          <article-title>Stаtіstісаl аnаlysіs оf lосаl plоts оf bіts sequenсes</article-title>
          .
          <source>Prоblemy uprаvlenіyа і іnfоrmаtіkі</source>
          ,
          <volume>5</volume>
          ,
          <fpage>92</fpage>
          -
          <lpage>105</lpage>
          (
          <year>2019</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Popereshnyak</surname>
            <given-names>S.</given-names>
          </string-name>
          <article-title>Analysis of pseudorandom small sequences using multidimensional statistics</article-title>
          .
          <source>In: The 3rd IEEE International Conference on Advanced Information and Communication Technologies (AICT)</source>
          , pp.
          <fpage>5</fpage>
          .
          <issue>4</issue>
          .1-
          <issue>5</issue>
          .4.4. IEEE Press, Ukraine (
          <year>2019</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Nejad</surname>
            <given-names>F. H.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Sabah</surname>
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Jam</surname>
            <given-names>A. J.</given-names>
          </string-name>
          <article-title>Analysis of avalanche effect on advance encryption standard by using dynamic S-Box depends on rounds keys</article-title>
          .
          <source>In: The 2014 International Conference on Computational Science and Technology (ICCST)</source>
          , pp.
          <fpage>1</fpage>
          -
          <lpage>5</lpage>
          . IEEE Press, Kota Kinabalu (
          <year>2014</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Bhaskar</surname>
            <given-names>C. U.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Rupa</surname>
            <given-names>C.</given-names>
          </string-name>
          <article-title>An advanced symmetric block cipher based on chaotic systems</article-title>
          . In:
          <article-title>The 2017 Innovations in Power and Advanced Computing Technologies (i-PACT)</article-title>
          , pp.
          <fpage>1</fpage>
          -
          <lpage>4</lpage>
          . IEEE Press, Vellore (
          <year>2017</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <surname>Busireddygari</surname>
            <given-names>P.</given-names>
          </string-name>
          ; Kak S.
          <article-title>Pseudorandom tableau sequences</article-title>
          ,
          <source>In: 51st Asilomar Conference on Signals, Systems, and Computers</source>
          , pp.
          <fpage>1733</fpage>
          -
          <lpage>1736</lpage>
          . IEEE Press (
          <year>2017</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <surname>Gurugopinath</surname>
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Samudhyatha</surname>
            <given-names>B.</given-names>
          </string-name>
          ,
          <article-title>Multi-dimensional Anderson-Darling statistic based goodness-of-fit test for spectrum sensing</article-title>
          .
          <source>In: Seventh International Workshop on Signal Design and its Applications in Communications (IWSDA)</source>
          . pp.
          <fpage>165</fpage>
          -
          <lpage>169</lpage>
          . Bengaluru, India. (
          <year>2015</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7. Moody D.
          <article-title>Post-quantum cryptography: NIST's plan for the future</article-title>
          .
          <source>In: Proceedings of the Seventh International Conference on Post Quantum Cryptography</source>
          . IEEE Press, Japan, (
          <year>2016</year>
          ). https://pqcrypto2016.jp
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <given-names>Special</given-names>
            <surname>Publication</surname>
          </string-name>
          800-
          <fpage>22</fpage>
          .
          <article-title>A Statistical Test Suite for Random and Pseudorandom Number Generators for Cryptographic Applications</article-title>
          . http://csrc. nist.gov
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>