DigitCycleLength[ r_Rational, b_Integer?Positive ] := MultiplicativeOrder[ b, FixedPoint[ Quotient[ #, GCD[ #, b ] ] &, Denominator[ r ] ] ]; a = 0; Do[ If[ Prime[ n ] - DigitCycleLength[ 1/Prime[ n ], 10 ] == 1, a++ ], {n, 2, PrimePi[ 10^7 ]} ] Print[ a ]