skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on March 1, 2026

Title: Advances in tabulating Carmichael numbers
Abstract We report that there are 49679870 Carmichael numbers less than$$10^{22}$$ 10 22 which is an order of magnitude improvement on Richard Pinch’s prior work. We find Carmichael numbers of the form$$n = Pqr$$ n = P q r using an algorithm bifurcated by the size ofPwith respect to the tabulation boundB. For$$P < 7 \times 10^7$$ P < 7 × 10 7 , we found 35985331 Carmichael numbers and 1202914 of them were less than$$10^{22}$$ 10 22 . When$$P > 7 \times 10^7$$ P > 7 × 10 7 , we found 48476956 Carmichael numbers less than$$10^{22}$$ 10 22 . We provide a comprehensive overview of both cases of the algorithm. For the large case, we show and implement asymptotically faster ways to tabulate compared to the prior tabulation. We also provide an asymptotic estimate of the cost of this algorithm. It is interesting that Carmichael numbers are worst case inputs to this algorithm. So, providing a more robust asymptotic analysis of the cost of the algorithm would likely require resolution of long-standing open questions regarding the asymptotic density of Carmichael numbers.  more » « less
Award ID(s):
2401305
PAR ID:
10626981
Author(s) / Creator(s):
;
Publisher / Repository:
SpringerNature
Date Published:
Journal Name:
Research in Number Theory
Volume:
11
Issue:
1
ISSN:
2522-0160
Page Range / eLocation ID:
8
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We state a general purpose algorithm for quickly finding primes in evenly divided sub-intervals. Legendre’s conjecture claims that for every positive integern, there exists a prime between$$n^2$$ n 2 and$$(n+1)^2$$ ( n + 1 ) 2 . Oppermann’s conjecture subsumes Legendre’s conjecture by claiming there are primes between$$n^2$$ n 2 and$$n(n+1)$$ n ( n + 1 ) and also between$$n(n+1)$$ n ( n + 1 ) and$$(n+1)^2$$ ( n + 1 ) 2 . Using Cramér’s conjecture as the basis for a heuristic run-time analysis, we show that our algorithm can verify Oppermann’s conjecture, and hence also Legendre’s conjecture, for all$$n\le N$$ n N in time$$O( N \log N \log \log N)$$ O ( N log N log log N ) and space$$N^{O(1/\log \log N)}$$ N O ( 1 / log log N ) . We implemented a parallel version of our algorithm and improved the empirical verification of Oppermann’s conjecture from the previous$$N = 2\cdot 10^{9}$$ N = 2 · 10 9 up to$$N = 7.05\cdot 10^{13} > 2^{46}$$ N = 7.05 · 10 13 > 2 46 , so we were finding 27 digit primes. The computation ran for about half a year on each of two platforms: four Intel Xeon Phi 7210 processors using a total of 256 cores, and a 192-core cluster of Intel Xeon E5-2630 2.3GHz processors. 
    more » « less
  2. Abstract The electricE1 and magneticM1 dipole responses of the$$N=Z$$ N = Z nucleus$$^{24}$$ 24 Mg were investigated in an inelastic photon scattering experiment. The 13.0 MeV electrons, which were used to produce the unpolarised bremsstrahlung in the entrance channel of the$$^{24}$$ 24 Mg($$\gamma ,\gamma ^{\prime }$$ γ , γ ) reaction, were delivered by the ELBE accelerator of the Helmholtz-Zentrum Dresden-Rossendorf. The collimated bremsstrahlung photons excited one$$J^{\pi }=1^-$$ J π = 1 - , four$$J^{\pi }=1^+$$ J π = 1 + , and six$$J^{\pi }=2^+$$ J π = 2 + states in$$^{24}$$ 24 Mg. De-excitation$$\gamma $$ γ rays were detected using the four high-purity germanium detectors of the$$\gamma $$ γ ELBE setup, which is dedicated to nuclear resonance fluorescence experiments. In the energy region up to 13.0 MeV a total$$B(M1)\uparrow = 2.7(3)~\mu _N^2$$ B ( M 1 ) = 2.7 ( 3 ) μ N 2 is observed, but this$$N=Z$$ N = Z nucleus exhibits only marginalE1 strength of less than$$\sum B(E1)\uparrow \le 0.61 \times 10^{-3}$$ B ( E 1 ) 0.61 × 10 - 3  e$$^2 \, $$ 2 fm$$^2$$ 2 . The$$B(\varPi 1, 1^{\pi }_i \rightarrow 2^+_1)/B(\varPi 1, 1^{\pi }_i \rightarrow 0^+_{gs})$$ B ( Π 1 , 1 i π 2 1 + ) / B ( Π 1 , 1 i π 0 gs + ) branching ratios in combination with the expected results from the Alaga rules demonstrate thatKis a good approximative quantum number for$$^{24}$$ 24 Mg. The use of the known$$\rho ^2(E0, 0^+_2 \rightarrow 0^+_{gs})$$ ρ 2 ( E 0 , 0 2 + 0 gs + ) strength and the measured$$B(M1, 1^+ \rightarrow 0^+_2)/B(M1, 1^+ \rightarrow 0^+_{gs})$$ B ( M 1 , 1 + 0 2 + ) / B ( M 1 , 1 + 0 gs + ) branching ratio of the 10.712 MeV$$1^+$$ 1 + level allows, in a two-state mixing model, an extraction of the difference$$\varDelta \beta _2^2$$ Δ β 2 2 between the prolate ground-state structure and shape-coexisting superdeformed structure built upon the 6432-keV$$0^+_2$$ 0 2 + level. 
    more » « less
  3. Abstract LetXbe a compact normal complex space of dimensionnandLbe a holomorphic line bundle onX. Suppose that$$\Sigma =(\Sigma _1,\ldots ,\Sigma _\ell )$$ Σ = ( Σ 1 , , Σ ) is an$$\ell $$ -tuple of distinct irreducible proper analytic subsets ofX,$$\tau =(\tau _1,\ldots ,\tau _\ell )$$ τ = ( τ 1 , , τ ) is an$$\ell $$ -tuple of positive real numbers, and let$$H^0_0(X,L^p)$$ H 0 0 ( X , L p ) be the space of holomorphic sections of$$L^p:=L^{\otimes p}$$ L p : = L p that vanish to order at least$$\tau _jp$$ τ j p along$$\Sigma _j$$ Σ j ,$$1\le j\le \ell $$ 1 j . If$$Y\subset X$$ Y X is an irreducible analytic subset of dimensionm, we consider the space$$H^0_0 (X|Y, L^p)$$ H 0 0 ( X | Y , L p ) of holomorphic sections of$$L^p|_Y$$ L p | Y that extend to global holomorphic sections in$$H^0_0(X,L^p)$$ H 0 0 ( X , L p ) . Assuming that the triplet$$(L,\Sigma ,\tau )$$ ( L , Σ , τ ) is big in the sense that$$\dim H^0_0(X,L^p)\sim p^n$$ dim H 0 0 ( X , L p ) p n , we give a general condition onYto ensure that$$\dim H^0_0(X|Y,L^p)\sim p^m$$ dim H 0 0 ( X | Y , L p ) p m . WhenLis endowed with a continuous Hermitian metric, we show that the Fubini-Study currents of the spaces$$H^0_0(X|Y,L^p)$$ H 0 0 ( X | Y , L p ) converge to a certain equilibrium current onY. We apply this to the study of the equidistribution of zeros inYof random holomorphic sections in$$H^0_0(X|Y,L^p)$$ H 0 0 ( X | Y , L p ) as$$p\rightarrow \infty $$ p
    more » « less
  4. Abstract The quantum simulation of quantum chemistry is a promising application of quantum computers. However, forNmolecular orbitals, the$${\mathcal{O}}({N}^{4})$$ O ( N 4 ) gate complexity of performing Hamiltonian and unitary Coupled Cluster Trotter steps makes simulation based on such primitives challenging. We substantially reduce the gate complexity of such primitives through a two-step low-rank factorization of the Hamiltonian and cluster operator, accompanied by truncation of small terms. Using truncations that incur errors below chemical accuracy allow one to perform Trotter steps of the arbitrary basis electronic structure Hamiltonian with$${\mathcal{O}}({N}^{3})$$ O ( N 3 ) gate complexity in small simulations, which reduces to$${\mathcal{O}}({N}^{2})$$ O ( N 2 ) gate complexity in the asymptotic regime; and unitary Coupled Cluster Trotter steps with$${\mathcal{O}}({N}^{3})$$ O ( N 3 ) gate complexity as a function of increasing basis size for a given molecule. In the case of the Hamiltonian Trotter step, these circuits have$${\mathcal{O}}({N}^{2})$$ O ( N 2 ) depth on a linearly connected array, an improvement over the$${\mathcal{O}}({N}^{3})$$ O ( N 3 ) scaling assuming no truncation. As a practical example, we show that a chemically accurate Hamiltonian Trotter step for a 50 qubit molecular simulation can be carried out in the molecular orbital basis with as few as 4000 layers of parallel nearest-neighbor two-qubit gates, consisting of fewer than 105non-Clifford rotations. We also apply our algorithm to iron–sulfur clusters relevant for elucidating the mode of action of metalloenzymes. 
    more » « less
  5. Abstract In a Merlin–Arthur proof system, the proof verifier (Arthur) accepts valid proofs (from Merlin) with probability 1, and rejects invalid proofs with probability arbitrarily close to 1. The running time of such a system is defined to be the length of Merlin’s proof plus the running time of Arthur. We provide new Merlin–Arthur proof systems for some key problems in fine-grained complexity. In several cases our proof systems have optimal running time. Our main results include:Certifying that a list ofnintegers has no 3-SUM solution can be done in Merlin–Arthur time$$\tilde{O}(n)$$ O ~ ( n ) . Previously, Carmosino et al. [ITCS 2016] showed that the problem has a nondeterministic algorithm running in$$\tilde{O}(n^{1.5})$$ O ~ ( n 1.5 ) time (that is, there is a proof system with proofs of length$$\tilde{O}(n^{1.5})$$ O ~ ( n 1.5 ) and a deterministic verifier running in$$\tilde{O}(n^{1.5})$$ O ~ ( n 1.5 ) time).Counting the number ofk-cliques with total edge weight equal to zero in ann-node graph can be done in Merlin–Arthur time$${\tilde{O}}(n^{\lceil k/2\rceil })$$ O ~ ( n k / 2 ) (where$$k\ge 3$$ k 3 ). For oddk, this bound can be further improved for sparse graphs: for example, counting the number of zero-weight triangles in anm-edge graph can be done in Merlin–Arthur time$${\tilde{O}}(m)$$ O ~ ( m ) . Previous Merlin–Arthur protocols by Williams [CCC’16] and Björklund and Kaski [PODC’16] could only countk-cliques in unweighted graphs, and had worse running times for smallk.Computing the All-Pairs Shortest Distances matrix for ann-node graph can be done in Merlin–Arthur time$$\tilde{O}(n^2)$$ O ~ ( n 2 ) . Note this is optimal, as the matrix can have$$\Omega (n^2)$$ Ω ( n 2 ) nonzero entries in general. Previously, Carmosino et al. [ITCS 2016] showed that this problem has an$$\tilde{O}(n^{2.94})$$ O ~ ( n 2.94 ) nondeterministic time algorithm.Certifying that ann-variablek-CNF is unsatisfiable can be done in Merlin–Arthur time$$2^{n/2 - n/O(k)}$$ 2 n / 2 - n / O ( k ) . We also observe an algebrization barrier for the previous$$2^{n/2}\cdot \textrm{poly}(n)$$ 2 n / 2 · poly ( n ) -time Merlin–Arthur protocol of R. Williams [CCC’16] for$$\#$$ # SAT: in particular, his protocol algebrizes, and we observe there is no algebrizing protocol fork-UNSAT running in$$2^{n/2}/n^{\omega (1)}$$ 2 n / 2 / n ω ( 1 ) time. Therefore we have to exploit non-algebrizing properties to obtain our new protocol.Certifying a Quantified Boolean Formula is true can be done in Merlin–Arthur time$$2^{4n/5}\cdot \textrm{poly}(n)$$ 2 4 n / 5 · poly ( n ) . Previously, the only nontrivial result known along these lines was an Arthur–Merlin–Arthur protocol (where Merlin’s proof depends on some of Arthur’s coins) running in$$2^{2n/3}\cdot \textrm{poly}(n)$$ 2 2 n / 3 · poly ( n ) time.Due to the centrality of these problems in fine-grained complexity, our results have consequences for many other problems of interest. For example, our work implies that certifying there is no Subset Sum solution tonintegers can be done in Merlin–Arthur time$$2^{n/3}\cdot \textrm{poly}(n)$$ 2 n / 3 · poly ( n ) , improving on the previous best protocol by Nederlof [IPL 2017] which took$$2^{0.49991n}\cdot \textrm{poly}(n)$$ 2 0.49991 n · poly ( n ) time. 
    more » « less