uniform distribution theory

Similar documents
uniform distribution theory

THE EULER FUNCTION OF FIBONACCI AND LUCAS NUMBERS AND FACTORIALS

When do the Fibonacci invertible classes modulo M form a subgroup?

Diagonal Representation of Certain Matrices

Normic continued fractions in totally and tamely ramified extensions of local fields

Report Documentation Page

Sums of the Thue Morse sequence over arithmetic progressions

Use of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications

Crowd Behavior Modeling in COMBAT XXI

Real-Time Environmental Information Network and Analysis System (REINAS)

THE EULER FUNCTION OF FIBONACCI AND LUCAS NUMBERS AND FACTORIALS

P. Kestener and A. Arneodo. Laboratoire de Physique Ecole Normale Supérieure de Lyon 46, allée d Italie Lyon cedex 07, FRANCE

REGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT)

Predicting Tropical Cyclone Formation and Structure Change

Closed-form and Numerical Reverberation and Propagation: Inclusion of Convergence Effects

A report (dated September 20, 2011) on. scientific research carried out under Grant: FA

System Reliability Simulation and Optimization by Component Reliability Allocation

Attribution Concepts for Sub-meter Resolution Ground Physics Models

Report Documentation Page

CRS Report for Congress

Estimation of Vertical Distributions of Water Vapor and Aerosols from Spaceborne Observations of Scattered Sunlight

Analysis Comparison between CFD and FEA of an Idealized Concept V- Hull Floor Configuration in Two Dimensions. Dr. Bijan Khatib-Shahidi & Rob E.

Metrology Experiment for Engineering Students: Platinum Resistance Temperature Detector

Development and Application of Acoustic Metamaterials with Locally Resonant Microstructures

SW06 Shallow Water Acoustics Experiment Data Analysis

PIPS 3.0. Pamela G. Posey NRL Code 7322 Stennis Space Center, MS Phone: Fax:

Broadband matched-field source localization in the East China Sea*

Report Documentation Page

Quantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS)

Babylonian resistor networks

Understanding Near-Surface and In-cloud Turbulent Fluxes in the Coastal Stratocumulus-topped Boundary Layers

Sensitivity of West Florida Shelf Simulations to Initial and Boundary Conditions Provided by HYCOM Data-Assimilative Ocean Hindcasts

DIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION

FRACTAL CONCEPTS AND THE ANALYSIS OF ATMOSPHERIC PROCESSES

Exact Solution of a Constrained. Optimization Problem in Thermoelectric Cooling

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION REPORT NUMBER

USMC Enlisted Endstrength Model

The Lehmer matrix and its recursive analogue

REPORT DOCUMENTATION PAGE

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models

Thermo-Kinetic Model of Burning for Polymeric Materials

Ocean Acoustics Turbulence Study

DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited.

Coastal Engineering Technical Note

Marginal Sea - Open Ocean Exchange

High-Fidelity Computational Simulation of Nonlinear Fluid- Structure Interaction Problems

Mine Burial Studies with a Large Oscillating Water-Sediment Tunnel (LOWST)

Understanding Near-Surface and In-Cloud Turbulent Fluxes in the Coastal Stratocumulus-Topped Boundary Layers

Advanced Numerical Methods for NWP Models

An Examination of 3D Environmental Variability on Broadband Acoustic Propagation Near the Mid-Atlantic Bight

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models

Improved Parameterizations Of Nonlinear Four Wave Interactions For Application In Operational Wave Prediction Models

LAGRANGIAN MEASUREMENTS OF EDDY CHARACTERISTICS IN THE CALIFORNIA CURRENT

SMA Bending. Cellular Shape Memory Structures: Experiments & Modeling N. Triantafyllidis (UM), J. Shaw (UM), D. Grummon (MSU)

Determining the Stratification of Exchange Flows in Sea Straits

A 1/10th Degree Global Ocean Simulation Using the Parallel Ocean Program

On Applying Point-Interval Logic to Criminal Forensics

REPORT DOCUMENTATION PAGE. Theoretical Study on Nano-Catalyst Burn Rate. Yoshiyuki Kawazoe (Tohoku Univ) N/A AOARD UNIT APO AP

Mixture Distributions for Modeling Lead Time Demand in Coordinated Supply Chains. Barry Cobb. Alan Johnson

Molecular Characterization and Proposed Taxonomic Placement of the Biosimulant 'BG'

HYCOM Caspian Sea Modeling. Part I: An Overview of the Model and Coastal Upwelling. Naval Research Laboratory, Stennis Space Center, USA

Optimizing Robotic Team Performance with Probabilistic Model Checking

Convection and Shear Flow in TC Development and Intensification

NAVGEM Platform Support

Scattering of Internal Gravity Waves at Finite Topography

Wavelet Spectral Finite Elements for Wave Propagation in Composite Plates

Mass Transport by Second Mode Internal Solitary Waves

Estimation of Vertical Distributions of Water Vapor from Spaceborne Observations of Scattered Sunlight

VLBA IMAGING OF SOURCES AT 24 AND 43 GHZ

The Extratropical Transition of Tropical Cyclones

High Resolution Surface Characterization from Marine Radar Measurements

Predictive Model for Archaeological Resources. Marine Corps Base Quantico, Virginia John Haynes Jesse Bellavance

Grant Number: N IP To compare obtained theoretical results with NPAL experimental data.

Analysis of South China Sea Shelf and Basin Acoustic Transmission Data

Volume 6 Water Surface Profiles

DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited.

Generation and Propagation of Internal Solitary Waves on the Continental Shelf and Slope

Abyssal Current Steering of Upper Ocean Current Pathways in an Ocean Model with High Vertical Resolution

INFRARED SPECTRAL MEASUREMENTS OF SHUTTLE ENGINE FIRINGS

Improvements in Modeling Radiant Emission from the Interaction Between Spacecraft Emanations and the Residual Atmosphere in LEO

Extension of the BLT Equation to Incorporate Electromagnetic Field Propagation

Comparative Analysis of Flood Routing Methods

AIR FORCE RESEARCH LABORATORY Directed Energy Directorate 3550 Aberdeen Ave SE AIR FORCE MATERIEL COMMAND KIRTLAND AIR FORCE BASE, NM

Surface Fluxes and Wind-Wave Interactions in Weak Wind Conditions

Design for Lifecycle Cost using Time-Dependent Reliability

Evolution of Tropical Cyclone Characteristics

Computer Simulation of Sand Ripple Growth and Migration.

A GENERAL PROCEDURE TO SET UP THE DYADIC GREEN S FUNCTION OF MULTILAYER CONFORMAL STRUCTURES AND ITS APPLICATION TO MICROSTRIP ANTENNAS

Curating and Maintenance of the Sediment Library and Dredge Collection

Improvement of Mesoscale Numerical Weather Prediction For Coastal Regions of Complex Terrain FY2003

Rogue Wave Statistics and Dynamics Using Large-Scale Direct Simulations

ANALYSIS AND MODELING OF STRATOSPHERIC GRAVITY WAVE ACTIVITY ALONG ER-2 FLIGHT TRACKS

THIRD-ORDER FAMILY OF METHODS IN BANACH SPACES

Assimilation of Synthetic-Aperture Radar Data into Navy Wave Prediction Models

Experimental and Theoretical Studies of Ice-Albedo Feedback Processes in the Arctic Basin

Dynamics of Droplet-Droplet and Droplet-Film Collision. C. K. Law Princeton University

Effects of Constrictions on Blast-Hazard Areas

Experimental and Theoretical Studies of Ice-Albedo Feedback Processes in the Arctic Basin

Vortex Rossby Waves and Hurricane Evolution in the Presence of Convection and Potential Vorticity and Hurricane Motion

Study of Electromagnetic Scattering From Material Object Doped Randomely WithThin Metallic Wires Using Finite Element Method

Transcription:

Uniform Distribution Theory 9 (2014), no.1, 21 25 uniform distribution theory ON THE FIRST DIGITS OF THE FIBONACCI NUMBERS AND THEIR EULER FUNCTION Florian Luca Pantelimon Stănică ABSTRACT. Here, we show that given any two finite strings of base b digits, say s 1 and s 2, there are infinitely many Fibonacci numbers F n such that the base b representation of F n starts with s 1 and the base b representation of φ(f n) starts with s 2. Communicated by Oto Strauch 1. Introduction Let b 2 be and integer. Let s 1 = c 1 c k(b) be a positive integer s 1 written in base b. Washington [4] proved that there exist infinitely many Fibonacci numbers F n whose base b representation starts with s 1. In fact, the first digits of the Fibonacci sequence obey Benford s law in that the proportion of the positive integers n such that F n starts with s 1 is precisely log((s 1 + 1)/s 1 )/. Here, we take this one step further. Let φ(m) be the Euler function of the positive integer m. We put s 2 = d 1... d l(b) for some other positive integer written in base b and prove the following theorem. Theorem. Given positive integers s 1 = c 1 c k(b) and s 2 = d 1... d l(b) written in base b, there exist infinitely many positive integers n such that the base b representation of F n starts with the digits of s 1 and the base b representation of φ(f n ) starts with the digits of s 2. We use the fact that with (α, β) = ((1+ 5)/2, (1 5)/2) the Binet formula F n = αn β n α β holds for all n 0. For a positive real number x we write log x for the natural logarithm of x, and x, respectively, {x}, for the integer part, respectively, fractional part of x. 2010 M a t h e m a t i c s S u b j e c t C l a s s i f i c a t i o n: 11B39,11K36. K e y w o r d s: Euler function, Fibonacci numbers, digits distribution. 21

Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to a penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number. 1. REPORT DATE 2014 2. REPORT TYPE 3. DATES COVERED 00-00-2014 to 00-00-2014 4. TITLE AND SUBTITLE On the First Digits of the Fibonacci Numbers and Their Euler Function 5a. CONTRACT NUMBER 5b. GRANT NUMBER 5c. PROGRAM ELEMENT NUMBER 6. AUTHOR(S) 5d. PROJECT NUMBER 5e. TASK NUMBER 5f. WORK UNIT NUMBER 7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) Naval Postgraduate School,Department of Applied Mathematics,Monterey,CA,93943 8. PERFORMING ORGANIZATION REPORT NUMBER 9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSOR/MONITOR S ACRONYM(S) 12. DISTRIBUTION/AVAILABILITY STATEMENT Approved for public release; distribution unlimited 13. SUPPLEMENTARY NOTES 11. SPONSOR/MONITOR S REPORT NUMBER(S) 14. ABSTRACT Here, we show that given any two nite strings of base b digits say s1 and s2, there are in nitely many Fibonacci numbers Fn such that the base b representation of Fn starts with s1 and the base b representation of (Fn) starts with s2. 15. SUBJECT TERMS 16. SECURITY CLASSIFICATION OF: 17. LIMITATION OF ABSTRACT a. REPORT unclassified b. ABSTRACT unclassified c. THIS PAGE unclassified Same as Report (SAR) 18. NUMBER OF PAGES 5 19a. NAME OF RESPONSIBLE PERSON Standard Form 298 (Rev. 8-98) Prescribed by ANSI Std Z39-18

F. LUCA P. STĂNICĂ 2. The proof By replacing s 1 with s 1 b m for some positive integer m, if needed, whose effect is adding m zeros at the end of the base b representation of s 1, we may assume that s 1 > s 2. By replacing s 1, s 2 by s 1 b m, respectively, s 2 b m for an arbitrary positive integer m, we may assume that the length of the base b representation of s 1, that is k, is as large as we wish. In Section 4 of [2], it is shown that φ(f n )/F n is dense in [0, 1]. So, we take ε (0, 1/(15b 2k )) and choose a positive integer a such that φ( ) ( s2 + ε, s ) 2 + 2ε. (1) s 1 s 1 Now we take any prime p > and look at p. Since p >, it follows that ( ) Fap p =, and the two factors and p / on the right above are coprime (indeed, the only common prime factor of these two numbers could be p, which is not the case since p > ). Any prime factor q of p / is a primitive prime factor of F dp for some divisor d of a. Recall that a prime number q is said to be a primitive prime factor of F n if q divides F n, but does not divide any F m for 1 m < n. One of the properties of primitive prime factors q of F n when n > 5 is that q ±1 (mod n). In particular, every prime factor q of p / is congruent to ±1 (mod p). Let q 1,..., q t be all the prime factors of p /. Then (2p 1) t q 1 q t p p α ap. Thus, t = O(p/ log p). Then ( ) t ) φ(p ) φ(fa ) = (1 1qi p 22 = φ() = φ() = φ() i=1 t 1 exp + O 1 q i=1 i q 2 q q1 ( ( )) t exp O q 1 ( ( )) 1 exp O log p

FIRST DIGITS OF THE FIBONACCI NUMBERS AND THEIR EULER FUNCTION = φ() ( ( )) 1 1 + O. log p It implies that, if p > exp(κε 1 ), where κ > 0 is some absolute constant, then ( φ(p ) s2 + 0.5ε, s ) 2 + 1.5ε. (2) p s 1 s 1 We now follow Washington s argument [4] to prove that there exist infinitely many primes p such that the base b representation of p starts with s 1. For this, it is enough to show that p = s 1 b N + ζ ap for some integer 0 ζ ap b N 1. (3) Note that since q 1 2p 1, it follows that if p is sufficiently large (say, p > b k ), then p cannot equal s 1 b N, and in particular, if in the above formula (3) we have ζ ap 0, then in fact ζ ap 1. The above formula (3) yields α ap = 5s 1 b N + 5ζ ap + β ap = 5s 1 b N (1 + x ap ). Since ζ ap 1, it follows that 5ζ ap + β ap > 5 1 > 1, and 5ζap + β ap 0 < x ap = 5s1 b. N So, if x ap (0, 1/b k ), and p > b k is sufficiently large, it then follows that ζ ap < b N, which is what we want. Thus, or ap log α = log( + N + log(1 + x ap ), ap log α N log( { = log( } + log(1 + x ap). (4) Observe that log( / is never an integer. Assume that k is sufficiently large such that { 1 b k < 1 log( }. Then putting δ = log(1 + 1/bk ), we see that a relation like (4) with x ap (0, 1/b k ) holds provided that { ( )} ({ a log α log( } { log( } ) p, + δ. (5) 23

F. LUCA P. STĂNICĂ The number γ = a log α/ is irrational. By a result of Vinogradov [3], the sequence of fractional parts {pγ} p prime is uniformly distributed. In particular, containment (5) holds for a positive proportion of primes p, and therefore certainly for infinitely many of them. So, indeed relation (3) holds. Relation (2) now shows that φ(p ) = s 2 b N + θ, where θ ( ( ) ( ) ) s2 ζ ap + 0.5ε + 0.5εb N s2, ζ ap + 1.5ε + 1.5εs 1 b N. s 1 s 1 Since ε < 1/(15b 2k ), the above upper bound is ( ) ( s2 b ζ ap + 1.5ε + 1.5εs 1 b N < (b N k 1 1) s 1 b k + 0.1 ) b k < b N 1, where the last inequality above is implied by 1 9 < bn 1 b N, + 0.1(bk 1) b 2k b N which holds true for all b 2 and N 1. This completes the proof of the theorem. 3. Comments It was shown in [1] that with σ(m) being the sum of divisors of the positive integer m, the ratio σ(f n )/F n is dense in [1, ). The present method now shows that there are infinitely many positive integers n such that the base b representation of F n starts with the digits of s 1 and the base b representation of σ(f n ) starts with the digits of s 2. Also, one may replace the Fibonacci sequence F n in the above statements with some other sequence u n for which it has been proved that φ(u n )/u n and u n /σ(u n ), respectively, are dense in [0, 1]. For example, one can take u n = 2 n 1 (see [1]) and the main result of this paper still holds provided that b is not a power of 2. We give no further details. 24

FIRST DIGITS OF THE FIBONACCI NUMBERS AND THEIR EULER FUNCTION 4. Acknowledgements We thank the referee for a careful reading of the paper and for comments which improved its quality. This paper was written during a visit of F. L. to the Applied Mathematics Department of Naval Postgraduate School in December, 2012. During the preparation of this paper, F. L. was supported in part by Project PAPIIT IN104512 (UNAM), VSP N62909-12-1-4046 (Department of the US Navy, ONR Global) and a Marcos Moshinsky fellowship. REFERENCES [1] F. Luca, On the sum of divisors of the Mersenne numbers, Math. Slovaca 53 (2003), 457 466. [2] F. Luca, V. J. Mejía Huguet and F. Nicolae, On the Euler function of Fibonacci numbers, J. Integer Sequences 9 (2009), A09.6.6. [3] I. M. Vinogradov, A new estimation of a trigonometric sum containing primes, Bull. Acad. Sci. URSS Ser. Math. 2 (1938), 1 13. [4] L. C. Washington, Benford s law for Fibonacci and Lucas numbers, Fibonacci Quart. 19 (1981), 175 177. Received January 10, 2012 Accepted June 19, 2013 F. Luca Mathematical Institute, UNAM Juriquilla 76230 Santiago de Querétaro, México and School of Mathematics, University of the Witwatersrand P. O. Box Wits 2050, South Africa E-mail: fluca@matmor.unam.mx P. Stănică Naval Postgraduate School Applied Mathematics Department Monterey, CA 93943, USA E-mail: pstanica@nps.edu 25