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.


Title: A deterministic algorithm for finding r-power divisors
Abstract Building on work of Boneh, Durfee and Howgrave-Graham, we present a deterministic algorithm that provably finds all integers p such that $$p^r \mathrel {|}N$$ p r | N in time $$O(N^{1/4r+\epsilon })$$ O ( N 1 / 4 r + ϵ ) for any $$\epsilon > 0$$ ϵ > 0 . For example, the algorithm can be used to test squarefreeness of N in time $$O(N^{1/8+\epsilon })$$ O ( N 1 / 8 + ϵ ) ; previously, the best rigorous bound for this problem was $$O(N^{1/6+\epsilon })$$ O ( N 1 / 6 + ϵ ) , achieved via the Pollard–Strassen method.  more » « less
Award ID(s):
1946311
PAR ID:
10420421
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Research in Number Theory
Volume:
8
Issue:
4
ISSN:
2522-0160
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present a new elementary algorithm that takes $$ \textrm{time} \ \ O_\epsilon \left( x^{\frac{3}{5}} (\log x)^{\frac{8}{5}+\epsilon } \right) \ \ \textrm{and} \ \textrm{space} \ \ O\left( x^{\frac{3}{10}} (\log x)^{\frac{13}{10}} \right) $$ time O ϵ x 3 5 ( log x ) 8 5 + ϵ and space O x 3 10 ( log x ) 13 10 (measured bitwise) for computing $$M(x) = \sum _{n \le x} \mu (n),$$ M ( x ) = ∑ n ≤ x μ ( n ) , where $$\mu (n)$$ μ ( n ) is the Möbius function. This is the first improvement in the exponent of x for an elementary algorithm since 1985. We also show that it is possible to reduce space consumption to $$O(x^{1/5} (\log x)^{5/3})$$ O ( x 1 / 5 ( log x ) 5 / 3 ) by the use of (Helfgott in: Math Comput 89:333–350, 2020), at the cost of letting time rise to the order of $$x^{3/5} (\log x)^2 \log \log x$$ x 3 / 5 ( log x ) 2 log log x . 
    more » « less
  2. We present a weighted approach to compute a maximum cardinality matching in an arbitrary bipartite graph. Our main result is a new algorithm that takes as input a weighted bipartite graph G(A cup B,E) with edge weights of 0 or 1. Let w <= n be an upper bound on the weight of any matching in G. Consider the subgraph induced by all the edges of G with a weight 0. Suppose every connected component in this subgraph has O(r) vertices and O(mr/n) edges. We present an algorithm to compute a maximum cardinality matching in G in O~(m(sqrt{w} + sqrt{r} + wr/n)) time. When all the edge weights are 1 (symmetrically when all weights are 0), our algorithm will be identical to the well-known Hopcroft-Karp (HK) algorithm, which runs in O(m sqrt{n}) time. However, if we can carefully assign weights of 0 and 1 on its edges such that both w and r are sub-linear in n and wr=O(n^{gamma}) for gamma < 3/2, then we can compute maximum cardinality matching in G in o(m sqrt{n}) time. Using our algorithm, we obtain a new O~(n^{4/3}/epsilon^4) time algorithm to compute an epsilon-approximate bottleneck matching of A,B subsetR^2 and an 1/(epsilon^{O(d)}}n^{1+(d-1)/(2d-1)}) poly log n time algorithm for computing epsilon-approximate bottleneck matching in d-dimensions. All previous algorithms take Omega(n^{3/2}) time. Given any graph G(A cup B,E) that has an easily computable balanced vertex separator for every subgraph G'(V',E') of size |V'|^{delta}, for delta in [1/2,1), we can apply our algorithm to compute a maximum matching in O~(mn^{delta/1+delta}) time improving upon the O(m sqrt{n}) time taken by the HK-Algorithm. 
    more » « less
  3. We consider the problem of designing sublinear time algorithms for estimating the cost of minimum] metric traveling salesman (TSP) tour. Specifically, given access to a n × n distance matrix D that specifies pairwise distances between n points, the goal is to estimate the TSP cost by performing only sublinear (in the size of D) queries. For the closely related problem of estimating the weight of a metric minimum spanning tree (MST), it is known that for any epsilon > 0, there exists an O^~(n/epsilon^O(1))-time algorithm that returns a (1+epsilon)-approximate estimate of the MST cost. This result immediately implies an O^~(n/epsilon^O(1)) time algorithm to estimate the TSP cost to within a (2 + epsilon) factor for any epsilon > 0. However, no o(n^2)-time algorithms are known to approximate metric TSP to a factor that is strictly better than 2. On the other hand, there were also no known barriers that rule out existence of (1 + epsilon)-approximate estimation algorithms for metric TSP with O^~ (n) time for any fixed epsilon > 0. In this paper, we make progress on both algorithms and lower bounds for estimating metric TSP cost. On the algorithmic side, we first consider the graphic TSP problem where the metric D corresponds to shortest path distances in a connected unweighted undirected graph. We show that there exists an O^~(n) time algorithm that estimates the cost of graphic TSP to within a factor of (2 − epsilon_0) for some epsilon_0 > 0. This is the first sublinear cost estimation algorithm for graphic TSP that achieves an approximation factor less than 2. We also consider another well-studied special case of metric TSP, namely, (1, 2)-TSP where all distances are either 1 or 2, and give an O^~(n ^ 1.5) time algorithm to estimate optimal cost to within a factor of 1.625. Our estimation algorithms for graphic TSP as well as for (1, 2)-TSP naturally lend themselves to O^~(n) space streaming algorithms that give an 11/6-approximation for graphic TSP and a 1.625-approximation for (1, 2)-TSP. These results motivate the natural question if analogously to metric MST, for any epsilon > 0, (1 + epsilon)-approximate estimates can be obtained for graphic TSP and (1, 2)-TSP using O^~ (n) queries. We answer this question in the negative – there exists an epsilon_0 > 0, such that any algorithm that estimates the cost of graphic TSP ((1, 2)-TSP) to within a (1 + epsilon_0)-factor, necessarily requires (n^2) queries. This lower bound result highlights a sharp separation between the metric MST and metric TSP problems. Similarly to many classical approximation algorithms for TSP, our sublinear time estimation algorithms utilize subroutines for estimating the size of a maximum matching in the underlying graph. We show that this is not merely an artifact of our approach, and that for any epsilon > 0, any algorithm that estimates the cost of graphic TSP or (1, 2)-TSP to within a (1 + epsilon)-factor, can also be used to estimate the size of a maximum matching in a bipartite graph to within an epsilon n additive error. This connection allows us to translate known lower bounds for matching size estimation in various models to similar lower bounds for metric TSP cost estimation. 
    more » « less
  4. This is the second in a pair of works which study small disturbances to the plane, periodic 3D Couette flow in the incompressible Navier-Stokes equations at high Reynolds number Re . In this work, we show that there is constant 0 > c 0 ≪ 1 0 > c_0 \ll 1 , independent of R e \mathbf {Re} , such that sufficiently regular disturbances of size ϵ ≲ R e − 2 / 3 − δ \epsilon \lesssim \mathbf {Re}^{-2/3-\delta } for any δ > 0 \delta > 0 exist at least until t = c 0 ϵ − 1 t = c_0\epsilon ^{-1} and in general evolve to be O ( c 0 ) O(c_0) due to the lift-up effect. Further, after times t ≳ R e 1 / 3 t \gtrsim \mathbf {Re}^{1/3} , the streamwise dependence of the solution is rapidly diminished by a mixing-enhanced dissipation effect and the solution is attracted back to the class of “2.5 dimensional” streamwise-independent solutions (sometimes referred to as “streaks”). The largest of these streaks are expected to eventually undergo a secondary instability at t ≈ ϵ − 1 t \approx \epsilon ^{-1} . Hence, our work strongly suggests, for all (sufficiently regular) initial data, the genericity of the “lift-up effect ⇒ \Rightarrow streak growth ⇒ \Rightarrow streak breakdown” scenario for turbulent transition of the 3D Couette flow near the threshold of stability forwarded in the applied mathematics and physics literature. 
    more » « less
  5. Abstract The Markoff graphs modulopwere proven by Chen (Ann Math 199(1), 2024) to be connected for all but finitely many primes, and Baragar (The Markoff equation and equations of Hurwitz. Brown University, 1991) conjectured that they are connected for all primes, equivalently that every solution to the Markoff equation moduloplifts to a solution over$$\mathbb {Z}$$ Z . In this paper, we provide an algorithmic realization of the process introduced by Bourgain et al. [arXiv:1607.01530] to test whether the Markoff graph modulopis connected for arbitrary primes. Our algorithm runs in$$o(p^{1 + \epsilon })$$ o ( p 1 + ϵ ) time for every$$\epsilon > 0$$ ϵ > 0 . We demonstrate this algorithm by confirming that the Markoff graph modulopis connected for all primes less than one million. 
    more » « less