跳转至

title: 英文全文 · FrontierMath: A Benchmark for Evaluating Advanced Mathematical Reasoning in AI paper_id: 'arXiv `2411.04872v7' year: '2025'


FrontierMath: A Benchmark for Evaluating Advanced Mathematical Reasoning in AI

Source: converted from arXiv LaTeX source source/main.tex.

-0.5cm

Abstract

We introduce FrontierMath, a benchmark of hundreds of original, exceptionally challenging mathematics problems crafted and vetted by expert mathematicians. The questions cover most major branches of modern mathematics—from computationally intensive problems in number theory and real analysis to abstract questions in algebraic geometry and category theory. Solving a typical problem requires multiple hours of effort from a researcher in the relevant branch of mathematics, and for the upper end questions, multiple days. FrontierMath uses new, unpublished problems and automated verification to reliably evaluate models while minimizing risk of data contamination. Current state-of-the-art AI models solve under 2% of problems, revealing a vast gap between AI capabilities and the prowess of the mathematical community. As AI systems advance toward expert-level mathematical abilities, FrontierMath offers a rigorous testbed that quantifies their progress.

Teams interested in evaluating models on the benchmark should reach out to mailto:math_evals@epochai.orgmath_evals@epochai.org We thank Terence Tao, Timothy Gowers, and Richard Borcherds for their insightful interviews and thank Terence Tao for contributing several problems to this benchmark. We also thank Nat Friedman, Fedor Petrov, Feiyang Lin, Eli Lifland, Marc Pickett, and the Constellation team. We gratefully acknowledge OpenAI for their support in creating the benchmark.

Introduction

Recent AI systems have demonstrated remarkable proficiency in tackling challenging mathematical tasks, from achieving olympiad-level performance in geometry trinh2024solving to improving upon existing research results in combinatorics romera2024mathematical. However, existing benchmarks face some limitations:

Saturation of existing benchmarks

Current standard mathematics benchmarks such as the MATH dataset hendrycks2021measuring and GSM8K cobbe2021training primarily assess competency at the high-school and early undergraduate level. As state-of-the-art models achieve near-perfect performance on these benchmarks, we lack rigorous ways to evaluate their capabilities in advanced mathematical domains that require deeper theoretical understanding, creative insight, and specialized expertise. Furthermore, to assess AI's potential contributions to mathematics research, we require benchmarks that better reflect the challenges faced by working mathematicians.

Benchmark contamination in training data

A significant challenge in evaluating large language models (LLMs) is data contamination—the inadvertent inclusion of benchmark problems in training data. This issue leads to artificially inflated performance metrics that mask models' true reasoning capabilities deng2023investigating, xu2024benchmark. While competitions like the International Mathematics Olympiad (IMO) or American Invitational Mathematics Examination (AIME) offer fresh test problems created after model training, making them immune to contamination, these sources provide only a small number of problems and often require substantial manual effort to grade.

Figure fig:sample-problems: FrontierMath features hundreds of advanced mathematics problems that require hours for expert mathematicians to solve. We release representative samples from the benchmark. Statements and solutions for released problems are provided in sec:sample_problems.

See Appendix A.3 for precise coefficients. To address these limitations, we introduce FrontierMath, a benchmark of original, exceptionally challenging mathematical problems created in collaboration with over 60 mathematicians from leading institutions. FrontierMath addresses data contamination concerns by featuring exclusively new, previously unpublished problems. It enables scalable evaluation through automated verification that provides rapid, reproducible assessment. The benchmark spans the full spectrum of modern mathematics, from challenging competition-style problems to problems drawn directly from contemporary research, covering most branches of mathematics in the 2020 Mathematics Subject Classification (MSC2020).In particular, FrontierMath contains problems across 70% of the top-level subjects in the MSC2020 classification excluding 00 General and overarching topics; collections",01 History and biography", and ``97 Mathematics education".

Comparison of unsolved problems across five mathematics benchmarks. While existing benchmarks are approaching saturation (see Appendix sub:percentage-unsolved-app), FrontierMath maintains a >98% unsolved rate.

Figure fig:benchmark_comparison: Comparison of unsolved problems across five mathematics benchmarks. While existing benchmarks are approaching saturation (see Appendix sub:percentage-unsolved-app), FrontierMath maintains a >98% unsolved rate.

The problems in FrontierMath demand deep theoretical understanding, creative insight, and specialized expertise, often requiring multiple hours of effort from expert mathematicians to solve. We provide a set of representative example problems and their solutions (see Appendix sub:percentage-unsolved-app and fig:sample-problems). Notably, current state-of-the-art AI models are unable to solve more than 2% of the problems in FrontierMath, even with multiple attempts, highlighting a significant gap between human and AI capabilities in advanced mathematics (see Figure fig:benchmark_comparison).

To understand expert perspectives on FrontierMath's difficulty and relevance, we interviewed several prominent mathematicians, including Fields Medalists Terence Tao (2006), Timothy Gowers (1998), and Richard Borcherds (1998), and IMO coach Evan Chen. They unanimously characterized the problems as exceptionally challenging, requiring deep domain expertise and significant time investment to solve. While they noted that the requirement of numerical answers differs from typical mathematical research practice, these problems nonetheless demand substantial expertise and time even from experienced mathematicians.

Data collection

The collection pipeline

We developed the FrontierMath benchmark through collaboration with over 60 mathematicians from universities across more than a dozen countries. Our contributors formed a diverse group spanning academic ranks from graduate students to faculty members. Many are distinguished participants in prestigious mathematical competitions, collectively holding 14 IMO gold medals. One contributing mathematician also holds a Fields Medal.

Their areas of expertise collectively span a vast expanse of modern mathematics—including but not limited to algebraic geometry, number theory, point set and algebraic topology, combinatorics, category theory, representation theory, partial differential equations, probability theory, differential geometry, logic, and theoretical computer science. Since many of our contributors are active researchers, the problems naturally incorporate sophisticated techniques and insights found in contemporary mathematical research.

Each mathematician created new problems following guidelines designed to ensure clarity, verifiability, and definitive answers across various difficulty levels and mathematical domains. Four requirements guided problem creation: itemize Originality: While problems could build on existing mathematical ideas, they needed to do so in novel and non-obvious ways—either through clever adaptations that substantially transform the original concepts, or through innovative combinations of multiple ideas that obscure their origins. This ensures that solving them requires genuine mathematical insight rather than pattern matching against known problems. Automated verifiability: Problems had to possess definitive, computable answers that could be automatically verified. To facilitate automated verification, we often structured problems to have integer solutions, which are straightforward to check programmatically. We also included solutions that could be represented as arbitrary SymPy objects—including symbolic expressions, matrices, and more. By utilizing SymPy, we enabled a wide range of mathematical outputs to be represented and verified programmatically, ensuring that even complex symbolic solutions could be efficiently and accurately checked. This approach allowed us to support diverse answer types while maintaining verification standards. Guessproofness: Problems were designed to avoid being susceptible to guessing; solving the problem had to be necessary to find the correct answer. As can be seen from fig:sample-problems, problems often have numerical answers that are large and nonobvious, which reduces the vulnerability to guesswork. As a rule of thumb, we require that there should not be a greater than 1% chance of guessing the correct answer without doing most of the work that one would need to do to ``correctly" find the solution. This is important to ensure that the challenges assessed genuine mathematical understanding and reasoning.

Computational tractability: The solution of a computationally intensive problem must include scripts demonstrating how to find the answer, starting only from standard knowledge of the field.For example, if a problem requires factoring a large integer, this must be doable from a standard algorithm. These scripts must cumulatively run less than a minute on standard hardware. This requirement ensures efficient and manageable evaluation times. itemize

For each contributed problem, mathematicians provided a detailed solution and multiple forms of metadata. Each submission included both a verification script for automated answer checking (see sec:automated_verification) and comprehensive metadata tags. These tags included subject and technique classifications, indicating both the mathematical domains and specific methodologies required for solution. Each submission then underwent peer review by at least one mathematician with relevant domain expertise, who evaluated all components: the problem statement, solution, verification code, difficulty ratings, and classification tags (see sec:validation). Additionally, the problems were securely produced and handled to prevent data contamination, and reviewed for originality (see sec:originality). Finally, the submissions included structured difficulty ratings that assessed three key dimensions: the required level of background knowledge, the estimated time needed to identify key insights, and the time required to work through technical details (see sec:problem_difficulty).

Automated verification

FrontierMath focuses exclusively on problems with automatically verifiable solutions—primarily numerical answers or mathematical objects that can be expressed as SymPy objects (including symbolic expressions, matrices, sets, and other mathematical structures). For each problem, the evaluated model submits code that computes and saves its answer as a Python object. A script then automatically verifies this answer. If the problem has a unique integer solution, it checks for an exact match. If the problem has a unique symbolic real answer, it uses SymPy evaluation to check that the difference between the submitted answer and the actual answer simplifies to 0. In all other cases, a custom verification script is necessary to check validity of a submitted answer (see Figure fig:verification_script).

Example verification script for an instance of Pell's equation: Find a tuple of integers $(x, y)$ such that $x^2 - 7y^2 = 1$. The left boxes show the code that produces and saves answers, the right box shows the verification script that evaluates them. A custom verification script is necessary here since the problem has infinitely many solutions.

Figure fig:verification_script: Example verification script for an instance of Pell's equation: Find a tuple of integers $(x, y)$ such that $x^2 - 7y^2 = 1$ . The left boxes show the code that produces and saves answers, the right box shows the verification script that evaluates them. A custom verification script is necessary here since the problem has infinitely many solutions.

This design choice enables rapid, objective evaluation of AI models without requiring them to output solutions in formal proof languages like Lean or Coq. The automated verification framework substantially reduces the cost and complexity of evaluation, as verification can be performed programmatically without human intervention. Furthermore, it eliminates potential human bias and inconsistency in the evaluation process, as the verification criteria are explicitly defined and uniformly applied.

While this approach enables efficient and scalable evaluation, it does impose certain limitations. Most notably, we cannot include problems that require mathematical proofs or formal reasoning steps, as these would demand human evaluation to assess correctness and clarity.In theory, we could have accepted formalized proofs in Lean or other formal languages as a problem solution format. However, this presents two challenges. Firstly, current AI models are often not trained to write formalized proofs in specialized languages. We wanted to ensure that the problems were measuring genuine reasoning skills, rather than skill at writing formalized proofs. Secondly, Lean's mathematical library, mathlib, doesn't fully cover the undergraduate math curriculum mathlib, let alone the scope of research math, which limits the fields such a benchmark could measure. However, even without formal proofs, these problems span a substantial breadth of modern mathematics and include problems that eminent mathematicians have characterized as exceptionally difficult (see sec:interviews).

Validation and quality assurance

To maintain high quality standards and minimize errors in the benchmark, we established a rigorous multi-stage review process. Each submitted problem underwent blind peer review by at least one expert mathematician from our team, who evaluated the following criteria:

itemize Verifying the correctness of the problem statement and the provided solution. Checking for ambiguities or potential misunderstandings in the problem wording. Assessing the problem for guessproofness, ensuring that solutions cannot be easily obtained without proper mathematical reasoning. Confirming the appropriateness of the assigned difficulty ratings and metadata. itemize

Verifying correctness is the most important of these tasks. For problems with unique solutions—the majority of our benchmark—verifying correctness requires checking the mathematical argumentation, computations, and, if applicable, programming scripts the writer provided to justify their answer. For problems with non-unique solutions, e.g. those which ask for a solution to some system of Diophantine equations or for a Hamiltonian path in a given graph, the reviewer only checks that the provided verification script matches the problem requirements and test it on the answer.

All problems in the benchmark have undergone at least one blind peer review. Through this review process, we identify problems with incorrect answers as well as those that fail to meet our guessproofness criterion. Such problems are only included in the benchmark after authors successfully address these issues through revision.

During first reviews, our reviewers commonly identify several types of issues with submitted problems. These include answers that aren't easily verifiable, problems where guessing is easier than proving, and cases where simple brute-force methods circumvent the intended difficulty. Such issues occur predominantly with authors who are new to submitting questions to the benchmark and have not yet fully internalized the guidelines. Beyond these concerns, straightforward errors in question statements or computational mistakes in solutions leading to incorrect answers are also not uncommon.

We commissioned blind reviews from second reviewers on 30 randomly chosen questions so far to get some idea about the noise threshold we might expect due to errors, and we also randomly selected 5 questions to be removed from the dataset and serve as a public sample of problems (see sec:sample_problems). In total, 35 of the accepted problems have received scrutiny beyond the first round of peer review.

Two out of 35 questions had an incorrect answer provided by the author, undetected in the first review. Assuming that a question being incorrect is a Bernoulli trial with probability ( p ) and that the second reviewers catch all errors, using a Jeffreys prior on ( p ) yields a posterior error rate of ( 2 + 0.535 + 1 6.9% ) for this benchmark. On one hand, the Jeffreys prior is conservative since our prior belief about ( p ) is likely lower due to the careful construction of the benchmark. On the other hand, we must account for potential undetected errors that even the second review might have missed. Therefore, estimating a critical error rate of approximately ( 10% ) is reasonable. This estimate aligns with error rates typical in machine learning benchmarks; for instance, [northcutt2021pervasivelabelerrorstest] estimate a ( >6% ) label error rate in the ImageNet validation set, while [gema2024we] find that over (9%) of questions in the MMLU benchmark contain errors based on expert review of 3,000 randomly sampled questions.

The second reviews also identified less critical errors. Six out of 35 questions had missing hypotheses in their statements which technically made them not fully rigorous. While domain experts could reasonably impute this missing context, it might still pose difficulties for models attempting to solve these problems. For two of the 35 questions, reviewers proposed strategies for guessing the solution with substantially less effort or computation than was necessary for a rigorous mathematical justification of the answer, violating the guessproofness criterion. In all cases, they proposed adjustments to the problem which corrected the error while preserving the original mathematical intent of the writer.

Although the error rate remains low, the mistakes we have spotted highlight the issues of a passive review system, in which the reviewer sees the full solution document and simply has to state their approval in order for a problem to be accepted into the benchmark. Going forward, we are adopting a more active review process, and will henceforth require concrete confirmation of a problem's soundness before it is accepted. For problems which seek a solution with some easily verified characteristic, such as a tuple of integers satisfying a Diophantine equation, it is enough to test that criterion directly. If a problem seeks a symbolic real number or some efficiently computed value, then it would be strong evidence if a heuristically determined guess approximates the given solution. For more abstract problems, it may be necessary to have a reviewer solve the problem themselves, given only the key steps and ideas of the solution rather than the full write-up.

Additionally, we observed inconsistent difficulty ratings between first and second reviewers; due to the subjective nature of this task, ratings rarely matched and often showed substantial differences.See sec:problem_difficulty for discussion on our ongoing effort to normalize these ratings. In several cases, second reviewers found no fundamental problems with the questions but suggested ways to increase their difficulty, indicating potential room to improve the difficulty of the remaining questions if reviewers were to spend more time on refining them.

Originality and data contamination prevention

Contributors were explicitly instructed to develop novel mathematical problems. While they could build upon existing mathematical ideas, they were required to do so in non-obvious ways—either through clever adaptations that substantially transformed the original concepts, or through novel and innovative combinations of multiple ideas. This approach ensured that solving the problems required genuine mathematical insight rather than pattern matching against known problems.

To minimize the risk of problems and solutions being disseminated online, we encouraged all submissions to be conducted through secure, encrypted channels. We employed encrypted communication platforms to coordinate with contributors and requested that any written materials stored online be encrypted—for instance, by placing documents in password-protected archives shared only via secure methods.While we endeavored to adhere strictly to these protocols, we acknowledge that, in some cases, standard email clients were used when communicating about a subset of the problems outside our immediate team.

Our primary method for validating originality relied on expert review by our core team of mathematicians. Their extensive familiarity with competition and research-level problems enabled them to identify potential similarities with existing problems that automated systems might miss. The team conducted thorough manual checks against popular mathematics websites, online repositories, and academic publications. This expert review process was supplemented by verification using plagiarism detection software.

To provide further confidence in our problems' originality, we ran them through the plagiarism detection tools Quetext and Copyscape. Across the entire dataset, the verification process revealed no significant matches to existing content, with minimal flags attributed only to standard mathematical terminology or software oversensitivity.We limited our plagiarism checking to problem statements only, excluding solutions to minimize online exposure. We specifically chose Quetext and Copyscape over more widely-used tools like Turnitin because they posed the lowest risk of our problems being stored in databases that might later be used to train AI models.

Problem difficulty

Estimates of mathematical problem difficulty are useful for FrontierMath's goal of evaluating AI mathematical capabilities, as such estimates would provide more fine-grained information about the performance of AI models. To assess problem difficulty, each problem author provided ratings along three key dimensions:

enumerate Background: This rating ranges from 1 to 5 and indicates the level of mathematical background required to approach the problem. A rating of 1 corresponds to high school level, 2 to early undergraduate level, 3 to late undergraduate level, 4 to graduate level, and 5 to research level. Creativity: Estimated as the number of hours an expert would need to identify the key ideas for solving the problem. This measure has no upper limit. Execution: Similarly estimated as the number of hours required to compute the final answer once the key ideas are found, including time writing a script if applicable. This measure also has no upper limit. enumerate

To verify and refine the initial difficulty assessments, each difficulty assessment underwent a peer-review process. Reviewers assessed the accuracy of the initial difficulty ratings. Any discrepancies between the authors' ratings and the reviewers' assessments were discussed, and adjustments were made as needed to ensure that the difficulty ratings accurately reflected the problem's demands in terms of background, creativity, and execution.

Accurately assessing the difficulty of mathematics problems presents significant challenges Chen2024, Chen2024Mistake. Problems that seem impossible may become trivial after exposure to certain techniques, and multiple solution paths of varying complexity often exist. Moreover, while we designed our problems to require substantial mathematical work rather than allow solutions through guessing or pattern matching, the possibility of models finding unexpected shortcuts could undermine such difficulty estimates.

Given these challenges, we view our difficulty ratings as providing rough guidance. More rigorous validation, such as through systematic data collection on human solution times, would be needed to make stronger claims about these difficulty assessments.In unpublished preliminary work, we found that problems rated as more difficult correlate with lower solution rates by GPT-4o, providing some support for our difficulty assessment system. However, more systematic validation would be needed to make strong claims about the reliability of these ratings.

Dataset composition

The FrontierMath benchmark covers a broad spectrum of contemporary mathematics, spanning both foundational areas and specialized research domains. This comprehensive coverage is crucial for effectively evaluating AI systems' mathematical capabilities across the landscape of modern mathematics. Working with over 60 mathematicians across different specializations, we captured most top-level MSC 2020 classification codes, reflecting the breadth of mathematics from foundational fields to emerging research areas.

Table tab:msc-distribution: Percentage distribution of Mathematics Subject Classification (MSC) 2020 codes in the FrontierMath dataset, representing the proportion of each classification relative to the total MSC assignments.

\begin{table}[h]
\centering
\small
\begin{tabularx}{\textwidth}{@{}p{6cm}r@{\hspace{0.75cm}}p{5.75cm}r@{}}
\toprule
\textbf{MSC Classification} & \textbf{Percentage} & \textbf{MSC Classification} & \textbf{Percentage} \\
\midrule
11 Number theory & 17.8\% & 57 Manifolds and cell complexes & 2.1\% \\
05 Combinatorics & 15.8\% & 13 Commutative algebra & 2.1\% \\
20 Group theory and generalizations & 8.9\% & 54 General topology & 1.4\% \\
60 Probability theory and stochastic processes & 5.1\% & 35 Partial differential equations & 1.4\% \\
15 Linear and multilinear algebra; matrix theory & 4.8\% & 53 Differential geometry & 1.4\% \\
14 Algebraic geometry & 4.8\% & 42 Harmonic analysis on Euclidean spaces & 1.4\% \\
33 Special functions & 4.8\% & 41 Approximations and expansions & 1.4\% \\
55 Algebraic topology & 3.1\% & 52 Convex and discrete geometry & 1.4\% \\
12 Field theory and polynomials & 2.4\% & 82 Statistical mechanics, structure of matter & 1.0\% \\
30 Functions of a complex variable & 2.4\% & 44 Integral transforms, operational calculus & 1.0\% \\
68 Computer science & 2.4\% & 17 Nonassociative rings and algebras & 1.0\% \\
18 Category theory; homological algebra & 2.4\% & Other (< 3 problems each) & 9.9\% \\
\bottomrule
\end{tabularx}
\caption{Percentage distribution of Mathematics Subject Classification (MSC) 2020 codes in the FrontierMath dataset, representing the proportion of each classification relative to the total MSC assignments.}
\label{tab:msc-distribution}
\end{table}

Number theory and combinatorics are most prominently represented, collectively accounting for approximately 34% of all MSC2020 tags. This prominence reflects both our contributing mathematicians' expertise and these domains' natural amenability to problems with numerical solutions that can be automatically verified. The next most represented fields are algebraic geometry and group theory (together about 14% of MSC tags), followed by algebraic topology (approximately 3%), linear algebra (5%), and special functions (5%). Problems involving probability theory and stochastic processes constitute about 5% of the MSC tags, with additional significant representation in complex analysis, category theory, and partial differential equations, each comprising between 1-3% of the MSC tags.

Mathematical subject interconnections in FrontierMath. Node sizes indicate the frequency of each subject's appearance in problems, while connections indicate when multiple mathematical subjects are combined within single problems, demonstrating the benchmark's integration of many mathematical domains.

Figure fig:msc-cooccurrence-network: Mathematical subject interconnections in FrontierMath. Node sizes indicate the frequency of each subject's appearance in problems, while connections indicate when multiple mathematical subjects are combined within single problems, demonstrating the benchmark's integration of many mathematical domains.

The network visualization (Figure fig:msc-cooccurrence-network) reveals how mathematical subjects combine within individual problems. Number theory and combinatorics appear together most frequently, with 13% of problems requiring both subjects, followed by combinatorics-group theory (9%) and number theory-group theory (8%). These three fields—number theory (44% of all problems), combinatorics (39%), and group theory (22%)—form the core of the benchmark, each combining with more than a dozen other mathematical domains in novel ways.

There is a wide range of techniques required to solve the problems in our dataset. In particular, there are over 200 different techniques listed as being involved in the solutions of our problems. Although generating functions, recurrence relations, and special functions emerge as common techniques, they each appear in less than 5% of problems, underscoring the benchmark's methodological diversity. Even the most frequently co-occurring techniques appear together in at most 3 problems, highlighting how problems typically require unique combinations of mathematical approaches.

Evaluation

Experiment-enabled evaluation framework

To evaluate AI models on FrontierMath problems, we developed a framework that allows models to explore and verify potential solution approaches through code execution, mirroring how mathematicians might experiment with ideas when tackling challenging problems. This framework enables models to test hypotheses, receive feedback, and refine their approach based on experimental results, as illustrated in Figure fig:problem-solving-framework.

The evaluation process follows an iterative cycle where the model analyzes the mathematical problem, proposes solution strategies, implements these strategies as executable Python code, receives feedback from code execution results, and refines its approach based on the experimental outcomes. For each problem, the model interacts with a Python environment where it can write code blocks that are automatically executed, with outputs and any error messages being fed back into the conversation. This allows the model to verify intermediate results, test conjectures, and catch potential errors in its reasoning.

Experimental framework for mathematical problem solving. The model can either provide a final answer directly or experiment with code execution before reaching a solution.

Figure fig:problem-solving-framework: Experimental framework for mathematical problem solving. The model can either provide a final answer directly or experiment with code execution before reaching a solution.

When submitting a final answer, models must follow a standardized format by including a specific marker comment (|# This is the final answer|), saving the result using Python's pickle module, and ensuring the submission code is self-contained and independent of previous computations. The interaction continues until either the model submits a correctly formatted final answer or reaches the token limit, which we set to 10,000 tokens in our experiments. If a model reaches this token limit without having submitted a final answer, it receives a final prompt requesting an immediate final answer submission. If the model still fails to provide a properly formatted final answer after this prompt, the attempt is marked as incorrect.

Results

Accuracy on FrontierMath

We evaluated six leading language models on our existing subset of FrontierMath problems: o1-preview ([o1-preview]), o1-mini ([o1-mini]), and GPT-4o (2024-08-06 version) ([GPT-4o]), Claude 3.5 Sonnet (2024-10-22 version) ([sonnet3-5]), Grok 2 Beta ([grokbeta]), and Google DeepMind's Gemini 1.5 Pro 002 ([gemini15pro002]). All models had a very low performance on FrontierMath problems, with no model achieving even a 2% success rate on the full benchmark (see Figure fig:model-performance). This stands in stark contrast to other mathematics evaluations such as GSM8K ([cobbe2021training]), MATH ([hendrycks2021measuring]), AIME 2024 ([o1-learningtoreason]), or Omni-MATH ([gao2024omni]), which are closer to saturation (see Figure fig:benchmark_comparison).

Performance of leading language models on FrontierMath based on mean accuracy across 8 runs. All models show consistently poor performance, with even the best models solving less than 2% of problems in each run on average. When re-evaluating problems that were solved at least once by any model, o1-preview demonstrated the strongest performance across repeated trials (see Section sub:atleastonce).

Figure fig:model-performance: Performance of leading language models on FrontierMath based on mean accuracy across 8 runs. All models show consistently poor performance, with even the best models solving less than 2% of problems in each run on average. When re-evaluating problems that were solved at least once by any model, o1-preview demonstrated the strongest performance across repeated trials (see Section sub:atleastonce).

Based on a single evaluation of the full benchmark, we found that models solved very few problems (less than 2% success rate). Given this low success rate and the fact that we ran only one evaluation, the precise ordering of model performance should be interpreted with significant caution, as individual successes can have an outsized impact on rankings. To better understand model behavior on solved problems, we identified all problems that any model solved at least once—a total of four problems—and conducted repeated trials with five runs per model per problem (see Appendix sub:atleastonce). We observed high variability across runs: only in one case did a model solve a question on all five runs (o1-preview for question 2). When re-evaluating these problems that were solved at least once, o1-preview demonstrated the strongest performance across repeated trials (see Section sub:atleastonce).

Moreover, even when a model obtained the correct answer, this does not mean that its reasoning was correct. For instance, on one of these problems running a few simple simulations was sufficient to make accurate guesses without any deeper mathematical understanding. However, models' low overall accuracy shows that such guessing strategies do not work on the overwhelming majority of FrontierMath problems. We also ran each model five times per problem on our five public sample problems (see Appendix sec:sample_problems), and made the transcripts publicly available.The transcripts can be downloaded at https://epochai.org/files/sample_question_transcripts.zip

Number of responses and token usage

Our evaluation framework allows models to run experiments and reflect on their results multiple times before submitting a final answer. We found significant differences in how models use this opportunity. For example, o1-preview averages just 1.29 responses per problem, whereas Grok 2 Beta averages 3.81 responses. o1-preview and Gemini 1.5 Pro 002 typically submit a final answer before seeing any experimental results, despite explicit instructions that they can take time to experiment and wait until they are confident to make a final guess. In contrast, the other models show more extended interactions.

This pattern is also reflected in token usage. Gemini 1.5 Pro 002 reaches the 10,000 token limit in only 16.8% of questions, while Claude 3.5 Sonnet, GPT-4o and Grok 2 Beta exceed this limit in more than 45% of their attempts. On average, Gemini 1.5 Pro uses the fewest tokens per question (around 6,000), while other models use between 12,000 and 17,000 tokens.These averages can exceed the 10,000 token conversation limit because models that reach the limit still get to see the results of their final experiment and complete their final answer.

Related work

The evaluation of mathematical reasoning in AI systems has progressed through various benchmarks and competitions, each testing different aspects of mathematical capability. Here we review the major existing benchmarks and competitions, and discuss how FrontierMath advances this landscape.

Several established mathematics benchmarks focus on evaluating basic mathematical reasoning capabilities in AI systems. Datasets like GSM8K cobbe2021training and MATH hendrycks2021measuring provide structured collections of elementary to undergraduate-level problems. While these benchmarks provided useful measures of capabilities in mathematical reasoning, state-of-the-art models now achieve near-perfect performance on many of them (see fig:benchmark_comparison), limiting their utility for discriminating between model capabilities.

Recent work has developed more challenging benchmarks for evaluating mathematical reasoning. The Advanced Reasoning Benchmark (ARB) sawada2023arb combines university-level mathematics with contest problems, spanning undergraduate to early graduate-level topics. The GHOSTS dataset frieder2024mathematical provides a comprehensive evaluation of language models' mathematical capabilities through 709 problems spanning graduate-level topics, with a focus on natural language interactions and human-expert evaluation. Several benchmarks have focused specifically on olympiad-level mathematics: OlympiadBench he2024olympiadbench compiled 8,476 problems from mathematics and physics olympiads, PutnamBench tsoukalas2024putnambench draws from the Putnam Competition, and OmniMATH gao2024omni provides 4,428 problems sourced from international competitions (including IMO), Art of Problem Solving forums, and verified competition archives, with detailed domain categorization and difficulty scaling. While these benchmarks present more challenging problems than earlier datasets, their reliance on historical competition materials creates significant data contamination risks, as models may have been trained on the source materials or highly similar problems.

The Putnam-AXIOM benchmark gulati2024putnamaxiom partially circumvents these contamination risks by generating problem variants. Specifically, the authors programmatically generate new problems of equivalent difficulty by renaming variables and changing constants in existing Putnam problems. This results in significant drops in performance - for example, o1-preview's accuracy drops from 50% to 34% using generated instead of unmodified problems. However, this technique only provides adds a limited degree of novelty to generated problems, and does not aid much in designing problems of much higher difficulty.

The AI Mathematical Olympiad (AIMO) xtx2024aimo, an initiative led by XTX Markets, aims to develop AI systems capable of achieving International Mathematical Olympiad gold medal performance. While AIMO creates novel national olympiad-level problems specifically designed to test advanced mathematical reasoning while avoiding data contamination risks, its requirement that runtime models be open-weight means that frontier AI models cannot be officially evaluated on the benchmark.

A parallel line of work explores formal mathematics benchmarks. The MiniF2F benchmark zheng2021minif2f provides 488 formally verified problems from mathematical olympiads and undergraduate courses, manually formalized across multiple theorem proving systems including Metamath, Lean, and Isabelle. While MiniF2F enables rigorous evaluation of neural theorem provers through machine-verified proofs, it faces similar data contamination concerns and focuses on competition-style problems rather than research mathematics.

FrontierMath addresses these fundamental limitations in three key ways. First, it eliminates data contamination concerns by featuring exclusively new, previously unpublished problems, ensuring models cannot leverage training on highly similar materials. Second, it raises the difficulty level beyond existing benchmarks by incorporating research-level mathematics across most branches of modern mathematics, moving beyond beyond both elementary problems and competition-style mathematics. Third, its automated verification system enables efficient, reproducible evaluation of any model—including closed source frontier AI systems—while maintaining high difficulty standards. This positions FrontierMath as a unique tool for evaluating progress toward research-level mathematical capabilities.

Interviews with mathematicians

We conducted interviews with four prominent mathematicians to gather expert perspectives on FrontierMath's difficulty, significance, and prospects: Terence Tao (2006 Fields Medalist), Timothy Gowers (1998 Fields Medalist), Richard Borcherds (1998 Fields Medalist), and Evan Chen, a renowned mathematics educator and IMO coach.Note that Evan Chen is also a co-author of this paper, and Terence Tao contributed several problems to the benchmark.

Assessment of FrontierMath difficulty

All four mathematicians characterized the research problems in the FrontierMath benchmark as exceptionally challenging, noting that the most difficult questions require deep domain expertise and significant time investment. For example, referring to a selection of several questions from the dataset, Tao remarked, "These are extremely challenging. I think that in the near term basically the only way to solve them, short of having a real domain expert in the area, is by a combination of a semi-expert like a graduate student in a related field, maybe paired with some combination of a modern AI and lots of other algebra packages...” However, some mathematicians pointed out that the numerical format of the questions feels somewhat contrived. Borcherds, in particular, mentioned that the benchmark problems “aren’t quite the same as coming up with original proofs.”

Expected timeline for FrontierMath progress

The mathematicians expressed significant uncertainty about the timeline for AI progress on FrontierMath-level problems, while generally agreeing these problems were well beyond current AI capabilities. Tao anticipated that the benchmark would "resist AIs for several years at least," noting that the problems require substantial domain expertise and that we currently lack sufficient relevant training data.

The interviewees anticipated human-AI collaboration would precede full automation. Chen and Tao both suggested that human experts working with AI systems could potentially tackle FrontierMath problems within around three years, much sooner than fully autonomous AI solutions. As Tao explained, on advanced graduate level questions, guiding current AI systems to correct solutions takes "about five times as much effort" as solving them directly—but he expected this ratio to improve below 1 on certain problems within a few years with sufficient tooling and capability improvements.

Key barriers to FrontierMath progress

A major challenge the mathematicians identified was the extremely limited training data available for FrontierMath's specialized research problems. As Tao noted, for many FrontierMath problems, the relevant training data is "almost nonexistent...you're talking like a dozen papers with relevant things." Gowers emphasized that solving these problems requires familiarity with "the tricks of the trade of some particular branch of maths"—domain knowledge that is difficult to acquire without substantial training data. The mathematicians discussed potential approaches to address data limitations, including synthetic data generation and formal verification methods to validate solutions.

Potential uses of FrontierMath-capable systems

The mathematicians interviewed identified several potential contributions from AI systems capable of solving FrontierMath-level problems. As Gowers noted, such systems could serve as powerful assistants for "slightly boring bits of doing research where you, for example, make some conjecture that would be useful, but you're not quite sure if it's true... it could be a very, very nice time saving device." Tao similarly emphasized their value for routine but technically demanding calculations, suggesting that a system performing at FrontierMath's level could help mathematicians tackle problems that currently "take weeks to figure out." However, several interviewees cautioned that the practical impact would depend heavily on computational costs. As Tao observed regarding current systems like AlphaProof, "if your amazing tool takes three days of compute off of all of Google to solve each problem...then that's less of a useful tool." The mathematicians generally agreed that AI systems at FrontierMath's level would be most valuable as supplements to human mathematicians rather than replacements, helping to verify calculations, test conjectures, and handle routine technical work while leaving broader research direction and insight generation to humans.

Discussion

In this paper, we introduced FrontierMath, a new benchmark composed of a challenging set of mathematical problems spanning most branches of modern mathematics. FrontierMath focuses on highly challenging math problems, with an emphasis on research-level problems, whose solutions require multiple hours of concentrated effort by expert mathematicians. We found that leading AI models cannot solve over 2% of the problems in the benchmark. We have further qualitatively validated the difficulty of the problems by sharing a sample of ten problems to four expert mathematicians, including three Fields medal winners, who uniformly assessed the problems as exceptionally difficult.

FrontierMath addresses two key limitations of previous math benchmarks such as the MATH dataset hendrycks2021measuring and GSM8K cobbe2021training. First, it focuses on exceptionally challenging problems that require deep reasoning and creativity, including research-level mathematics, preventing saturation with current leading AI systems. Second, by using exclusively novel, unreleased problems, it helps prevent data contamination, reducing the risk that models succeed through pattern matching against training data.

Despite these strengths, FrontierMath has several important limitations. First, the practical focus on automatically verifiable and numerical answers excludes proof-writing and open-ended exploration, which are significant parts of modern math research. Second, while our problems are substantially more challenging than existing benchmarks, requiring hours of focused work from expert mathematicians, they still fall short of typical mathematical research, which often spans weeks, months or even years of sustained investigation. By limiting problems to hours instead of months, we make the benchmark practical but miss testing crucial long-term research skills.

Another limitation is that current AI models cannot solve even a small fraction of the problems in our benchmark. While this demonstrates the high difficulty level of our problems, it temporarily limits FrontierMath's usefulness in evaluating relative performance of models. However, we expect this limitation to resolve as AI systems improve.

Evaluating AI systems' mathematical capabilities will provide valuable insights into their potential to assist with complex mathematical research. Beyond mathematics, the precise chains of reasoning and creative synthesis required in these problems serve as an important proxy for broader scientific thinking. Through its high difficulty standards and rigorous verification framework, FrontierMath offers a systematic way to measure progress in advanced reasoning capabilities as AI systems evolve.

Future work

We will further develop the FrontierMath benchmark by introducing new, rigorously vetted problems. These new problems will follow our established data collection and review process, with enhanced quality assurance procedures to reduce answer errors and improve the accuracy of difficulty ratings. We will also release another set of representative sample problems to showcase the benchmark's depth and variety.

We will further experiment with different evaluation methodologies to better assess the mathematical proficiency of today's leading models. For example, we will test the effects of increasing the token limit, allowing models to reason for longer and run more experiments per problem. We also plan to conduct multiple runs for each model-problem pair, enabling us to report statistics and confidence intervals across attempts. These evaluations will help us better understand how different evaluation parameters affect model performance on advanced mathematical reasoning tasks.

Through FrontierMath, we are setting a new standard of difficulty for evaluating complex reasoning capabilities in leading AI systems. At Epoch AI (www.epochai.org), we are committed to maintaining FrontierMath as a public resource and regularly conducting comprehensive evaluations of leading models. We will continue working with mathematicians to produce new problems, perform quality assurance, validate difficulty scores, and provide transparent and authoritative assessments of AI systems' mathematical reasoning abilities.

0.8em

Sample problems and solutions

This section presents five problems from the FrontierMath benchmark, with one chosen at random from each quintile of difficulty. The problems are arranged in decreasing order of difficulty, from research-level number theory to undergraduate algebraic geometry. Background levels indicate the prerequisite knowledge required, while creativity and execution ratings reflect hours of innovative thinking and technical rigor needed respectively.

Figure : MSC2020 codes: 11 (Number Theory), 14 (Algebraic Geometry), 20 (Group Theory and Generalizations)

We thank Piotr \'Sniady for pointing out a typo in the problem statement of Sample 4, and Will McLafferty, Timothy Smits, and Frederick Vu for identifying errors in the originally given solution for Sample 5 in an earlier draft of the paper.

Sample problem 1 — high difficulty

Problem

Definitions. For a positive integer $n$ , let $v_p(n)$ denote the largest integer $v$ such that $p^v \mid n$ .

For $p$ a prime and $a \not\equiv 0 \pmod{p}$ , we let $\ord_p(a)$ denote the smallest positive integer $o$ such that $a^o \equiv 1 \pmod{p}$ . For $x > 0$ , we let $$\ord_{p, x}(a) = \prod_{\substack{q \le x \ q \text{ prime}}} q^{v_q(\ord_p(a))} \prod_{{q > x \ q \text{ prime}}} q^{v_q(p-1)}.

$$

\textit{Problem.} Let $S_x$ denote the set of primes $p$ for which $$

_p, x(2) > _p, x(3),

$$ and let $d_x$ denote the density $$

d_x = |S_x||\p x : p is prime|

$$ of $S_x$ in the primes. Let $$

d_ = _x d_x.

$$

Compute\footnote{We have lowered the precision of this and several other problems which were written before the one-minute computation rule was finalized.} $\lfloor 10^6 d_{\infty} \rfloor.$

Answer:

367707

Solution

We will break this into multiple sections.

Background

First of all, the presence of the $x$-parameter is a technicality; informally, we are just asking for the density of primes $p$ such that $$

_p(2) > _p(3).

$$

The motivation of this problem comes from Artin's primitive conjecture which asks: for which integers $a$ are there infinitely many primes $p$ such that $a$ is a primitive root modulo $p$? This is equivalent to requesting $\ord_p(a) = p-1$, i.e. $\ind_p(a) = 1$.

Artin's conjecture is open. For our purposes, the relevant fact is that the conjecture has been solved on the assumption of a generalization of the Riemann hypothesis (for zeta functions of number fields). This was done by C. Hooley in 1967.

The key idea is that $a$ is a primitive root mod $p$ if and only if there is no $n > 1$ such that $a$ is an $n$th power modulo $p$ for some $n \mid p-1$. One can try to compute the number of such primes $p$ not exceeding $x$ with an inclusion-exclusion argument. The Chebotarev density theorem gives us asymptotics for

$$

|\p x : p 1 n, a is a perfect nth power modulo p|

$$

for any fixed $n.$

Bounding the error terms is what we need the generalization of Riemann's hypothesis for (c.f.: the standard RH bounds the error term in the prime number theorem); we don't know how to do it without.

For concreteness, if $a = 2$, then $a$ is a primitive root modulo $p$ for roughly $37\%$ of primes. See \url{https://mathworld.wolfram.com/ArtinsConstant.html}. In general, the density is strictly positive as long as $a$ is not a perfect square nor $-1$ (in which case there are no such primes $p > 3$). The density is not quite $37\%$ in all cases -- see below for more on this.

This argument is very general: it can be adapted to look at primes $p$ such that $\ord_p(2) = \ord_p(3) = p-1$, or $\ord_p(2) = (p-1)/a, \ord_p(3) = (p-1)/b$ for fixed $a, b \in \mathbb{Z}_+$, or $\ord_p(2) = \ord_p(3)$, or -- you guessed it -- $\ord_p(2) > \ord_p(3)$. The "bounding the error terms" part is no more difficult in these cases, and no non-trivial complications arise. But we of course still can't do it without GRH.

We replaced $\ord_p(a)$ with $\ord_{p, x}(2)$ for the reason of deleting the error terms, and for this reason only. As a result, when one does the inclusion-exclusion argument, there are only a finite number of terms needed, and thus we only need asymptotic formulas for \eqref{eq:1}, rather than any sharp bounds for the error term.

So the question then is about computing the main terms. To do so, one needs to be able to compute asymptotics for (1). In the case $a = 2$, it turns out that (1) is asymptotically $1/\phi(n)n$ for any square-free $n$. In some cases, like for $a = 4$ and $a = 5$, this is not the case: For $a = 4$, take $n = 2$. For $a = 5$, this fails for $n = 10$, since $p \equiv 1 \pmod{5}$ implies that $5$ is a square modulo $p$ -- this is related to the fact that $\sqrt{5} \in \mathbb{Q}(\zeta_5)$. But it turns out that there is always only a small, finite number of revisions one needs to make to the argument, allowing for the computation of the main terms.

There are a few papers considering generalizations of Artin's conjecture. It is known that the density of primes $p$ for which $\ord_p(2) > \ord_p(3)$ is positive (as a special case of a more general result), this again under GRH, but the task of computing that density has not been carried out in the literature, justifying the originality of this problem. There is a fair amount of legwork and understanding needed to do this.

Solution outline

Here are the key steps for verifying the correctness of the answer:

\begin{enumerate} \item We have \emph{explicit} (but cumbersome) formulas for the density $d(a, b) = d(a, b, x)$ of primes $p$ such that

$$

_p, x(2) = (p-1)/a, _p, x(3) = (p-1)/b

$$

\begin{enumerate} \item These formulas are "\emph{quasi-multiplicative}": they can be written in the form $$

d(a, b) = f(v_2(a), v_3(a), v_2(b), v_3(b)) _5 q x g_q(v_q(a), v_q(b)).

$$ \item $f$ and $g_q$ are explicit functions that can be evaluated in $O(1)$ time. \end{enumerate} \item For any fixed $a$, the sum $$

_a < b d(a, b, x)

$$ can be evaluated in $O(a)$ time. \begin{itemize} \item Idea: We can write $\sum_{a < b} d(a, b, x) = \sum_{b = 1}^{\infty} d(a, b, x) - \sum_{b = 1}^a d(a, b, x)$. The first sum is just the density of $p : \ord_{p, x}(2) = (p-1)/a$, which can be evaluated similarly to point 1 above. The second sum can be evaluated with brute force. \end{itemize} \item We have, for any $T \ge 100$, roughly\footnote{This omits a lower-order term.} $$

_T < a < b d(a, b, x) < 114T^-2 46/47,

$$ (think of the exponent as $-2 + \epsilon$), and thus we have $$

|_a < b d(a, b, x) - _a < b, a T d(a, b, x)| < 114T^-2 46/47.

$$ \item We have, for any $x_{\infty} > x \ge 10^8$, $$

a < b, a T |d(a, b, x) - d(a, b, x)| 19x.

$$ \end{enumerate}

Then, taking $T = 10^5$ and $x = 10^8$, we get an answer that has error at most $2.5 \cdot 10^{-8}$. This then gives $7$ correct decimal places.

Step 1

In this part, we fix $a, b \in \mathbb{Z}_+$, and we'll compute the density of primes $p$ satisfying \eqref{eq:2}.

For \eqref{eq:2} to hold, we need $2$ to be a perfect $a$th power (and $p \equiv 1 \pmod{a}$), but "no more": $2$ should not be an $aq$th power for any prime $q$ (such that $p \equiv 1 \pmod{aq}$). Similarly for $b$. It's convenient to frame these propositions in terms of algebraic number theory.

Let $K = \mathbb{Q}(\zeta_{\lcm(a, b)}, 2^{1/a}, 3^{1/b})$, and for each prime $q$ let $K_{1, q} = \mathbb{Q}(\zeta_{aq}, 2^{1/aq})$ and $K_{2, q} = \mathbb{Q}(\zeta_{bq}, 3^{1/bq})$. We want to compute the density of primes $p$ such that $p$ splits in $K$, but doesn't split in $K_{i, q}$ for any $i \in {1, 2}, 2 \le q \le x$.

There is an "obvious guess" for what that density should be, by naively assuming that different values of $q$ are independent (e.g. $2$ being a fifth power shouldn't provide evidence on whether $2$ is a seventh power), and that there are no "spurious correlations" (e.g. $2$ being a fifth power shouldn't provide evidence on whether $3$ is a fifth power or not). This intuition is almost correct: it makes a small, finite number of mistakes, namely at the primes $q = 2$ and $q = 3$. There are two reasons for this: \begin{itemize} \item $\zeta_3 \in \mathbb{Q}(\zeta_4, \sqrt{3})$, since $\zeta_3 = (1 + \sqrt{-3})/2$. \item $\sqrt{2} \in \mathbb{Q}(\zeta_8)$, since we have $(1 + i)^2 = 2i$, and thus $1 + i = \zeta_8 \cdot \sqrt{2}$. \end{itemize}

Crucially, the \emph{only} relations between $\zeta_N, 2^{1/N}$ and $3^{1/N}$ are these two. The following two well-known lemmas formalize this:

\textbf{Lemma 1.} Denote $K_n = \mathbb{Q}(\zeta_n, 2^{1/n}, 3^{1/n})$. Let $N$ be an arbitrary positive integer.\footnote{Picture $N$ as large, e.g. $N = 100$.} Let $q_1, q_2, \ldots , q_k$ be pairwise distinct primes greater than three. The fields $$

K_6^N, K_q_1^N, , K_q_k^N

$$ are linearly disjoint\footnote{This is equivalent to the Galois group of the compositum (in this case, $K_{(6q_1 \cdots q_k)^N}$) factorizing as the product of Galois groups of the individual fields: $\Gal(K_{(6q_1 \cdots q_k)^N}/\mathbb{Q}) \cong \Gal(K_{6^N}/\mathbb{Q}) \times \Gal(K_{q_1^N}/\mathbb{Q}) \times \cdots \times \Gal(K_{q_k^N}/\mathbb{Q})$. (Another way of saying this is that the degree of the compositum is the product of the degrees of the individual fields.) For the non-expert, this is just a way of saying "these things are independent from each other".}. Furthermore, the degrees $[K_{q_i^N} : \mathbb{Q}]$ are equal to $\phi(q_i^N)q_i^{2N}$ ("the maximum possible").

\textit{Proof.} We apply several results of \cite{olli}. Recall that $\mathbb{Q}(\zeta_{q_1^N}), \ldots , \mathbb{Q}(\zeta_{q_K^N})$ are linearly disjoint. We can append $K_{6^N}$ to the set while maintaining disjointness, since the largest cyclotomic subfield of $K_{6^N}$ lies in $\mathbb{Q}(\zeta_{2 \cdot 6^N}, \sqrt{2}, \sqrt{3}) \subset \mathbb{Q}(\zeta_{2 \cdot 6^N}, \zeta_8, \zeta_{12})$ (see Proposition 3.14 of the paper). Then, throwing in elements of the form $2^{1/q_i^N}$ and $3^{1/q_i^N}$ maintains linear disjointness, because the degrees of the resulting field extensions are powers of $q_i$, and those are coprime. For the claim about maximality of degrees, apply Proposition 3.9 with $n = q_i^N, K = \mathbb{Q}(\zeta_n), k = n^2$, with $a_i$ the numbers $2^{q_i^x}3^{q_i^y}$ for $0 \le x, y < N$; the quotients of such numbers are not $n$th powers in $K$ by Corollary 3.6. QED.

Of course, the claim about linear disjointness and the degrees being the maximum possible now holds for \emph{subfields} of $K_{6^N}, K_{q_1^N}, \ldots , K_{q_k^N}$, too; we'll be using this fact throughout.

At the primes $q = 2, q = 3$ we lose a factor of 4, as promised earlier:

\textbf{Lemma 2.} For $N \ge 3$ we have $$

[K_6^N : Q] = (6^N)6^2N4.

$$

\textit{Proof.} Apply Proposition 3.9 with $n = 6^N$, $K = \mathbb{Q}(\zeta_n), k = n^2/4$, with $a_i$ the numbers $2^{2^x3^y}3^{2^z3^w}$, where $0 \le x, y, z, w < N$ and $x, z < N/2$. The quotients of such numbers are not $n$th powers, again by Corollary 3.6. QED

With these tools, we'll now attack the problem of computing $d(a, b)$. Assume $x$ is large enough that $q \mid ab \implies q \le x$ (otherwise just replace $a \rightarrow \prod_{q \le x} q^{v_q(a)}$).

By the Chebotarev density theorem, the density of primes $p$ splitting in a field $K$ is $[K : \mathbb{Q}]$.\footnote{Reminder for the non-expert: Splitting of primes in a number field is a somewhat technical concept, but for our purposes you only need to know that: $p$ splitting in $\mathbb{Q}(\zeta_n)$ is equivalent to $p \equiv 1 \pmod{n}$, and $p$ splitting in a field like $\mathbb{Q}(\zeta_n, 2^{1/a}, 3^{1/b})$ (where we always assume $\lcm(a, b) \mid n$) is equivalent to $p \equiv 1 \pmod{n}$, $2$ being an $a$th power and $3$ being a $b$th power mod $p$. So "splitting" is exactly the thing we need.} Thus, by inclusion-exclusion,

$$

d(a, b, x) = n_1, n_2 : p n_i p x (n_1)(n_2)[Q((an_1, bn_2), 2^1/an_1, 3^1/bn_2) : Q].

$$

By the first lemma, this sum is "quasi-multiplicative" in the sense that prime powers $q^n$ with $q > 3$ can be factored out. Explicitly, if we denote $a' = 2^{v_2(a)}3^{v_3(a)}$ and define $b'$ similarly, we have

$$

d(a, b, x) = n_1, n_2 \1, 2, 3, 6\ (n_1)(n_2)[Q((a'n_1, b'n_2), 2^1/a'n_1, 3^1/b'n_2) : Q] 5 q x (_n_1, n_2 \1, q\ (n_1)(n_2)[Q((q^v_q(a)n_1, q^v_q(b)n_2), 2^1/q^v_q(a)n_1, 3^1/q^v_q(b)n_2) : Q]).

$$

The first sum over ${1, 2, 3, 6}^2$ depends only on $v_2(a), v_2(b), v_3(a), v_3(b)$, and the sums over ${1, q}^2$ only depend on $v_q(a), v_q(b)$. This proves 1(a).

For claim 1(b), we first note that for $q \ge 5$, the fields in the sum have "the maximum degree" (remember the remark about the conclusion of Lemma 1 applying to subfields, too): that is, we have $$

[Q(_(q^v_q(a)n_1, q^v_q(b)n_2), 2^1/q^v_q(a)n_1, 3^1/q^v_q(b)n_2)] = ((q^v_q(a)n_1, q^v_q(b)n_2))q^v_q(a) + v_q(b)n_1n_2.

$$ So these sums can be computed. (Writing the formulas down explicitly isn't necessary, since we can now write a program that computes the values.)

For the sum over $n_1, n_2 \in {1, 2, 3, 6}$, we need to be a bit more careful about the degrees to account for the two failures-in-maximality mentioned earlier. Denote the degree of $$

Q(_(a'n_1, b'n_2), 2^1/a'n_1, 3^1/a'n_2)

$$ by $D(a'n_1, b'n_2)$. By Lemma 2 (and by knowing where the factor of $4$ there comes from), we can compute $D(x, y)$. Indeed, we have \begin{itemize} \item $D(x, y) = \phi(\lcm(x, y))xy/4$, if we have ($8 \mid \lcm(x, y)$ and $2 \mid x$) \emph{and} ($12 \mid \lcm(x, y)$ and $2 \mid y$). \item $D(x, y) = \phi(\lcm(x, y))xy/2$, if we have \emph{either} ($8 \mid \lcm(x, y)$ and $2 \mid x$) \emph{or} ($12 \mid \lcm(x, y)$ and $2 \mid y$), but not both \item $D(x, y) = \phi(\lcm(x, y))xy$, if \emph{neither} ($8 \mid \lcm(x, y)$ and $2 \mid x$) \emph{nor} ($12 \mid \lcm(x, y)$ and $2 \mid y$) hold. \end{itemize} That is, we have $D(x, y) = \phi(\lcm(x, y))xy \cdot C(\gcd(x, 24), \gcd(y, 24))$, where $C : {1, 2, \ldots , 24}^2 \to {1, 1/2, 1/4}$ is explicitly given by the list above. Now, our sum can be written as

$$

_n_1, n_2 \1, 2, 3, 6\ (n_1)(n_2)((a'n_1, b'n_2))a'n_1b'n_2 C((a'n_1, 24), (b'n_2, 24)),

$$

which can be computed in $O(1)$ time.

Step 2

Now, we fix $a$, and aim to evaluate the sum $$

_a < b d(a, b, x) = _b = 1^ d(a, b, x) - _b = 1^a d(a, b, x).

$$ The second sum can be evaluated with brute force, and is thus not of interest to us here. Consider then the first sum. We can just forget about $b$ altogether here: because we are summing over all of $b$, we are imposing no constraint on $\ord_{p, x}(3)$; that is, we are only demanding $\ord_{p, x}(2) = (p-1)/a$.\footnote{Technical point: we're assuming here that $\lim_{N \to \infty} \sum_{a, b \le N} d(a, b, x) = 1$; see Step 3 for more on this.}

Now, the density $d(a, x)$ of $p$ such that $\ord_{p, x}(2) = (p-1)/a$ can be computed similarly to the computations in Part 1. But it's much simpler: note that $[\mathbb{Q}(\zeta_n, 2^{1/n}) : \mathbb{Q}]$ is always at least $n\phi(n)/2$ -- we lose only a factor of $2$, not a factor of $4$ like in Lemma 2, since we only have the failure at $8 \mid n$ coming from $\sqrt{2} \in \mathbb{Q}(\zeta_8)$. And so our formula is multiplicative:

$$

d(a, x) = _2 q x _n \1, q\ (n)[Q(_q^v_q(a)n, q^1/q^v_q(a)n) : Q].

$$

Furthermore, the terms are easy to calculate: If $q > 2$, then they are $$

1(q^v_q(a))q^v_q(a)(1 - 1q^2).

$$ For $q = 2$, the values are (for $v_q(a) = 0, 1, 2, \ldots$) $$

12, 38, 116, 364, 3256, 31024,

$$

Step 3

We now want to bound tail sums of $d(a, b, x)$; the basic idea is to give pointwise bounds for $d(a, b, x)$ and then sum those bounds.

\textbf{Lemma 3.} For any $a, b$ and $x$ we have $$

d(a, b, x) 4((a, b))ab.

$$

The proof is simple: the condition $\ord_{p, x}(2) = (p-1)/a, \ord_{p, x}(3) = (p-1)/b$ in particular implies that $p$ splits in $\mathbb{Q}(\zeta_{\lcm(a, b)}, 2^{1/a}, 3^{1/b})$, which has at least degree $\phi(\lcm(a, b))ab/4$ -- see the previous two lemmas.

Here is a very basic tail sum bound.

\textbf{Lemma 4.} We have $$

_T _a, b T d(a, b, x) = 1.

$$

For the proof, note that $$

_a, b : (a, b) > T d(a, b, x) _a, b : (a, b) > T 4((a, b))ab _a > T 1a(a) _1 b T 1b.

$$ The inner sum is $O(\log T)$. For the outer sum, we have well-known bound $\phi(a) \gg \sqrt{a}$, and so $\sum_{a > T} 1/(a\phi(a)) \ll 1/\sqrt{T}$, so the whole sum is $O(\log T / \sqrt{T}) = o(1)$.

Further note that $d(1, 2, x)$ is strictly positive (roughly $12.5\%$). Together with the previous lemma this implies that the density of $p$ with $\ord_{p, x}(2) > \ord_{p, x}(3)$ equals $$

(a, b) Z+^2 : a < b d(a, b, x).

$$ (An intuitive but non-trivial, technical claim -- we in fact haven't proven such results if $x$ is replaced with $\infty$, and this is the reason Artin's primitive conjecture is still open.)

Now, we want quantitative, strong bounds for $$

_T < a < b d(a, b, x).

$$ For this, we'll use the fact that $$

(n) c_1 (nc_2)^,

$$ where $c_1 = 92160, c_2 = 510510, \delta = 46/47$. (See \url{https://math.stackexchange.com/questions/301837/is-the-euler-phi-function-bounded-below}.) We have

$$

_T < a < b d(a, b, x) 4_T < a < b 1((a, b))ab < 2c_2^c_1 _T < a, b 1ab(a,b)^1 + .

$$

Continue:

$$

2c_2^c_1 _T < a, b 1ab(a, b)^ &= 2c_2^c_1_T < a, b (a, b)^(ab)^1 + & 2c_2^c_1_d = 1^ d^ _T/d < a', b' 1(da'db')^1 + &= 2c_2^c_1_d = 1^ d^-2 - (_T/d < n 1n^1 + )^2.

$$

The inner sum is bounded by $\int_{T/d}^{\infty} 1/t^{1 + \delta} dt = \delta^{-1}(T/d)^{-\delta}$ (plus $1$ for $d > T$). And thus, the above is bounded as $$

2c_2^c_1^2_d = 1^ d^-2 - ((dT)^2 + 1_d > T).

$$ The sum over $d$ without the $1_{d > T}$ part can be bounded numerically for $\delta = 46/47$ as $13$, and the $1_{d > T}$ part introduces an additional (small) contribution of $T^{-1 - \delta}/(1 + \delta)$, and so we arrive at $$

2c_2^c_1^2(13T^-2 + T^-1 - /(1 + )).

$$ For $T = 10^5$ this is less than $2 \cdot 10^{-8}$.

\begin{comment}

$$

_T < a < b d(a, b, x) 4_T < a < b 1((a, b))ab < 2c_1c_2^ _T < a, b 1ab(a,b)^1 + .

$$

Continue:

$$

2c_1c_2^ _T < a, b 1ab(a, b)^ &= 2c_1c_2^_T < a, b (a, b)^(ab)^1 + & 2c_1c_2^_d = 1^ d^ _T/d < a', b' 1(da'db')^1 + &= 2c_1c_2^_d = 1^ d^-2 - (_T/d < n 1n^1 + )^2.

$$

The inner sum is bounded by $\int_{T/d}^{\infty} 1/t^{1 + \delta} dt = \delta^{-1}(T/d)^{-\delta}$ (plus $1$ for $d > T$). And thus, the above is bounded as $$

2c_1c_2^^2_d = 1^ d^-2 - ((dT)^2 + 1_d > T).

$$ The sum over $d$ without the $1_{d > T}$ part can be bounded numerically for $\delta = 46/47$ as $13$, and the $1_{d > T}$ part introduces an additional (small) contribution of $T^{-1 + \delta}/(1 + \delta)$, and so we arrive at $$

2c_1c_2^^2(13T^-2 + T^-1+/(1 + )).

$$ For $T = 10^5$ this is less than $1.5 \cdot 10^{-9}$. \end{comment}

Step 4

We want to bound the difference that arises from replacing a particular value of $x$ in $d(a, b, x)$ with an arbitrarily large value of $x = x_{\infty}$. Recall the function $g_q(v_q(a), v_q(b))$ from Part 1. We have $g_q(0, 0) = 1 - 2/q(q-1) + 1/q^2(q-1)$. In particular, $$

_q > x is prime g_q(0, 0) _q > x is prime (1 - 2(q-1)^2).

$$ Now, when $0 \le t \le 1/100$, we have $1 - t \cdot 0.99 \ge e^{-t}$, so the above is bounded from below by $$

_q > x is prime e^-1.98/(q-1)^2 (_q > x is prime -1.98/(q-1)^2).

$$ For $y \ge 10^8$, there are at most $1.1y/\log(y)$ primes in the interval $[y, 2y)$, and thus by dyadic decomposition $$

_q > x is prime 1.98(q-1)^2 _y = 2^kx, k 0 1.98y^2 1.1y(y) = 1.98 1.1x _k 0 12^k(k(2) + (x)),

$$ which for $x \ge 10^8$ is numerically bounded by $1/9x$. Thus, we arrive at the lower bound $$

(-1/9x) 1 - 19x.

$$

Thus, for any $x_{\infty} > x > \max(a, b)$ (the last condition needed so that $v_q(a) = v_q(b) = 0$ for $q \ge x$) with $x \ge 10^8$, we have $$

0 d(a, b, x_) - d(a, b, x) d(a, b, x)(1 - 19x) - d(a, b, x) 19xd(a, b, x).

$$ And, since $\sum_{a < b, a \le T} d(a, b, x) \le 1$, we have $$

a < b, a T |d(a, b, x) - d(a, b, x)| 19x.$$

The code

Finally, we have posted to GitHub a C++ script for efficiently computing the answer. You can find it at https://github.com/epoch-research/artin_code/. The first 120 lines are prime sieves, gcd, fast exponentiation and other standard boilerplate. The next 80 lines are for computing various expressions arising in the explicit formulas for the density of eq:2, after which is a function for actually computing that density. This is done once more but one-dimensionally -- see Step 2. The main function computes the answer to the problem.

Sample problem 2 — high-medium difficulty

Problem

Construct a degree 19 polynomial $p(x) \in \mathbb{C}[x]$ such that $X:={p(x) = p(y)} \subset \mathbb{P}^1 \times \mathbb{P}^1$ has at least 3 (but not all linear) irreducible components over $\mathbb{C}$ . Choose $p(x)$ to be odd, monic, have real coefficients and linear coefficient -19 and calculate $p(19)$ .

Answer

1876572071974094803391179

Solution

Note that $p(x)$ defines a covering $p: \mathbb{P}x^1 \rightarrow \mathbb{P}^1_z$ of degree 19 which at $z=\infty$ has monodromy $\sigma$ , a 19-cycle, which without loss of generality is $(1, 2, 3, \hdots 19).$

Note that $X$ is just the fibre product of $p$ with $p$ (which in general is a singular curve) and therefore $p(x)$ also defines a finite map $\hat{p}: X \rightarrow \mathbb{P}^1_z$ of degree $19^2$ (if we prefer to stick with smooth schemes, we can instead use the normalisation $\hat{X}$ of $X$ , and we will freely move between the irreducible components of $X$ and the corresponding connected components of the smooth $\hat{X}$ with its covering structure).

The idea then is that any such covering $\hat{p}$ can be studied using group theory via its monodromy. Specifically, suppose $q: X \rightarrow \mathbb{P}_z^1$ is of degree $n$ and has branch locus $B \subset \mathbb{P}^1_z$ .Then if you follow the fibres around each branch value and see how this permutes the points in the fibres, we get a homomorphism $\bar{q}: \pi_1(\mathbb{P}^1_z-B,x) \rightarrow \text{Im}(\phi):=G_q \subset S_n$ (defined up to conjugation) and $X$ is irreducible exactly when $\bar{q}$ is transitive. Moreover, if $X$ is reducible, then the irreducible components correspond to the $G_q$ -orbits in $S_n$ .

Suppose $p$ with branch locus $B$ (including $\infty$ ) is encoded as $\bar{p}: G_p \subset S_{19}$ . Then we can check that $\hat{p}$ has branch locus $B$ and is encoded diagonally as $\bar{p}^2: G_p \subset S_{19}^2$ . So we are looking for some $G_p \subset S_{19}$ that has between 3 and 18 orbits in $S_{19}^2$ . In particular $G_p$ is not doubly transitive on $S_{19}$ .

Recall a theorem of Burnside burnside: for $q$ prime and a permutation group $G \subset S_q,$ either $G$ is doubly transitive or $|G|$ divides $q(q-1).$

Applying this and Sylow's theorem to our $G_p$ , we see that the order 19 subgroup $\langle \sigma_{\infty}\rangle \subset G_p$ is normal. Now take any $\sigma_i \in G_p$ and suppose $\sigma_i$ fixes at least 2 points (wlog 1 and $n$ ). Then by normality, $\sigma_{i} \sigma_{\infty} \sigma_i^{-1} = \sigma_{\infty}^s$ for some $s \in \mathbb{Z}$ . Writing this out explicitly, we must have $n = 1 + (n-1)s \text{ mod } p$ and hence $s=1$ . Therefore plugging in 1 to $\sigma_{i} \sigma_{\infty} \sigma_i^{-1} = \sigma_{\infty}$ gives $\sigma_i(2)=2$ . Continuing in this way, we conclude all non-trivial elements $\sigma_i$ must have at most 1 fixed point.

Applying Riemann-Hurwitz to $p$ , we have $19-1=\Sigma_r \text{Index}(\sigma_r)$ where $r$ enumerates the branch points in $B$ , $\sigma_r$ is the corresponding permutation over $r$ and Index geometrically is the ramification index-1. Since each non-trivial $\sigma_i$ fixes at most 1 point, we have that $\text{Index}(\sigma_i)\geq (19-1)/2=9$ , with the lower bound corresponding to precisely $9$ disjoint transpositions. Therefore Riemann-Hurwitz implies $r \leq 2$ . If $r=1$ , we have exactly one 19-cycle, which would give rise to 19 irreducible components. So the only possibility is $r=2$ and $\text{Index}(\sigma_i)=9$ i.e. $\sigma_i$ consist of 9 disjoint transpositions for $i=1,2$ . Explicitly if we assume wlog that $\sigma_1$ contains $(1,2)$ then the relation $\sigma_1 \sigma_2=\sigma_{\infty}$ (coming from the topology of $\mathbb{P}^1-B$ ) implies that $\sigma_2$ contains $(2,19)$ and continuing we find $\sigma_1 = (1,2)(3, 19)(4, 18) \hdots$ and $\sigma_2=(2,19)(3,18) \hdots$ .

By the Riemann existence theorem, a polynomial with specified monodromy is essentially unique (up to affine scaling of the target and source) and is given by the degree 19 Chebyshev polynomial $T_{19}(x)$ (which can be checked has this cycle structure). One checks that it has exactly the diagonal in $S_{19}^2$ and 9 other orbits, giving 10 irreducible components.

This proof therefore shows that the only possibility is a polynomial affine equivalent to $T_{19}(x)$ i.e. of the form $AT_{19}(ax+b)+B$ . Moreover, one can check that, except for the diagonal, all the remaining components are quadratic in both variables, hence of bi-degree (2,2).

Finally we impose the desired properties of $p(x)$ . Since $T_{19}(x)$ is odd, $p(x)$ being odd implies $b=0$ and $B=0$ . $p(x)$ monic implies $A * a^{19} * 2^{18}=1$ and the linear of coefficient of $p(x)$ is $-19Aa=-19$ . Solving gives $A=1/a$ and $a^{18}=1/2^{18}$ . Real coefficients means $a=\pm 1/2$ and hence $p(x)=2 * T_{19}(x/2)$ .

$p(x) = 2 T_{19}(x/2)$ , the degree 19 Chebyshev polynomial, which has 10 irreducible components (including the diagonal). Therefore $p(19)= 1876572071974094803391179.$

Sample problem 3 — medium difficulty

Problem

Let $a_n$ for $n\in\mathbb{Z}$ be the sequence of integers satisfying the recurrence formula [a_n=198130309625 a_n-1+354973292077 a_n-2-427761277677 a_n-3+370639957a_n-4] with initial conditions $a_i=i$ for $0\leq i\leq 3$ . Find the smallest prime $p\equiv 4\bmod7$ for which the function $\mathbb{Z} \rightarrow \mathbb{Z}$ given by $n\mapsto a_n$ can be extended to a continuous function on $\mathbb{Z}_p$ .

Answer

9811

Solution

The idea of the solution comes from the Skolem-Mahler-Lech theorem. Let $\alpha_1,\dots,\alpha_4$ be the roots of the characteristic polynomial of the recurrence, this is, [x^4-198130309625 x^3-354973292077 x^2+427761277677 x-370639957.] Then we have the explicit formula [a_n=c_1_1^n+c_2_2^n+ c_3_3^n+c_4_4^n] for all $n\in\mathbb Z$ , where the coefficients $c_i$ are determined by the initial conditions. We can think of the roots $\alpha_i$ as elements of some extension $L$ of $\mathbb Q_p$ of degree $\leq4!$ . Let $\pi$ be a uniformiser of $L$ , and denote by $v_L$ the $\pi$ -adic valuation normalised so that $v_L(p)=1$ (in particular, $v_L(\pi)=1/e)$ .

Note that the independent coefficient $370639957$ is prime, so all the roots $\alpha_i$ have $v_L(\alpha_i)=0$ for any prime $p<370639957$ . As we shall see, the prime we are looking for is indeed $<370639957$ .

Let us prove that whenever all the roots satisfy $v_L(\alpha_i-1)>0$ , then there exists the required extension of $n\mapsto a(n)$ to $\mathbb Z_p$ . If this is the case, then $\alpha_i\equiv1\bmod\pi$ , so in particular [-370639957=_1_2_3_41 p.] Hence, $p$ must divide $370639957+1=370639958$ . The prime factors of this number are $2$ , $13$ , $1453$ and $9811$ , of which only the last two are $\equiv4\bmod 7$ . Similarly, we have [198130309625_1+_2+_3+_44 p] so $198130309625-4=198130309621$ must be multiple of $p$ . The prime factors are $37$ , $673$ , $811$ , $9811$ , so the only possible prime is $p=9811$ .

For this prime the characteristic polynomial reduces mod $p$ to [x^4-4x^3+6x^2-4x+1=(x-1)^4] so indeed the four roots are $\equiv1\bmod p$ .

Let us show now how to do the extension of $n\mapsto a(n)$ . By the formula above, we need to extend $\alpha^n$ to a function on $\mathbb Z_p$ , and as in real analysis this is accomplished via the exponential and the logarithm. The main problem is that $\log x$ converges only for $v_p(x-1)>0$ in $\mathbb Z_p$ . In our case this is true, so we can simply write [a_n=_ic_i_i^n= _ic_i(n)] and this clearly extends to the continuous function [f(x)=_ic_i(x).] (Note that $f$ is unique since $\mathbb Z$ is dense in $\mathbb Z_p$ ).

For other primes $p$ for which the roots are not $1\bmod \pi$ , the extension to $\mathbb Z_p$ is not possible. This is because of the following. For simplicity let us assume that the roots $\alpha_i$ all lie in $\mathbb Z_p$ (this happens exactly when the characteristic polynomial splits completely modulo $p$ ). The general case is done similarly but with some care.

Since the projection map $\mathbb Z_p^\times\to\mathbb F_p^\times$ induces an isomorphism between the group of $(p-1)$ th roots of unities and $F_p^\times$ , for every $\alpha\in\mathbb Z_p^\times$ we can find some $(p-1)$ th roots of unity $\zeta_\alpha$ such that $\alpha\equiv\zeta_\alpha\bmod p$ . Hence, $v_p(\alpha/\zeta_\alpha-1)>0$ , so $\log(\alpha/\zeta_\alpha)$ exists. Hence, [a_n=ic_i_i^n= _ic_i__i^n(_i/__i)^n= _ic_i__i^n(n(_i/__i)),] and since $\zeta\alpha^{p-1}=1$ , we have [a_n=_ic_i__i^np-1(n(_i/__i)).] Thus, defining for each $0\leq k\leq p-2$ the function [a^(k)(x)=_ic_i__i^k(x(_i/__i))] then we have [a^(k)(n)=a_nif $n\equiv k\bmod{p-1}$ .] Thus, we can extend the function $n\mapsto a_n$ from the set $S_k={n\equiv k\bmod{p-1}}$ to its closure $\mathbb Z_p$ . But each of the functions $a^{(k)}(x)$ is different, so they are not compatible.

Sample problem 4 — medium-low difficulty

Problem

Let $M_{1000}^4$ be the set of $4$ -tuples of invertible $1000\times1000$ matrices with coefficients in $\mathbb{C}$ . Let $S\subset M_{1000}^4$ be the subset of all tuples $(A_1,A_2,A_3,A_4)$ satisfying the conditions: alignat2 A_i^2 &= I, && for all 1 i 4 A_iA_j &= A_jA_i, && if \3j- i, 3i-j\ 5Z_>0 = A_iA_jA_i^-1A_j^-1 &= A_jA_i, && if \3j- i, 3i-j\ 5Z_>0 , alignat where $5\mathbb{Z}{>0}$ refers to the set of positive multiples of $5$ , i.e., $5\mathbb{Z} = {5,10,15,\dots}$ .

The group $G = GL(1000)$ of invertible complex $1000\times 1000$ matrices acts on $S$ by the formula: [B (A_1,A_2,A_3,A_4) = (BA_1B^-1, BA_2B^-1, BA_3B^-1, BA_4B^-1).] Find the number of orbits of this action, i.e., find $|S/G|$ .

Answer

625243878951

Solution

Let $V$ be a $1000$ dimensional complex vector space, with some chosen basis. Then using the fact that the $A_i$ satisfy precisely the Coxeter relations for the symmetric group $S_5$ , we deduce that the data of $A_1,\dots, A_4$ is just a representation of $S_5$ in $V$ . More precisely, the set of $(i,j) \in {1,2,3,4}^2$ such that $3j-i\in 5\mathbb{Z}_{>0}$ is precisely ${(1,2), (2,4), (4,3)}$ . This shows that the "braid relation" only occurs for these three pairs, and all other pairs just give commutation relations.

Declaring that the transpositions $(1\ 2), (2\ 3), (3\ 4), (4\ 5)$ act on $V$ as multiplication respectively by $A_1, A_2, A_4$ and $A_3$ gives a representation of $S_5$ on $V$ . The action of $GL(1000)$ just gives the matrices of the same representation but with respect to different bases. Therefore, the number of orbits of this action is exactly the number of $1000$ dimensional representations of $S_5$ up to isomorphism.

From the character table of $S_5,$ we see that the irreducible representations of $S_5$ have dimensions $1,1,4,4,5,5,6$ . Thus the number we are looking for is the number of ways to write $1000$ as a sum of $1,4,5,6$ (with two different kinds of $1,4,5$ each). This is precisely the coefficient of $t^{1000}$ in the power series [(1+t+t^2+)^2(1+t^4+t^8+)^2(1+t^5+t^10+)^2(1+t^6+t^12+),] which is equal to the power series [1(1-t)^2(1-t^4)^2(1-t^5)^2(1-t^6).]

The above coefficient can be computed to be 625243878951.

Sample problem 5 — low difficulty

Problem

How many nonzero points are there on [x^3y+y^3z+z^3x=0]

over $\mathbb{F}_{5^{18}}$ up to scaling?

Answer

3814708984376

Solution

Using the code presented in fig:proj-solution-counter to count the solutions for $\mathbb{F}5,\mathbb{F}$ gives us $6,26,126$ solutions respectively. You can download the package finite_field_implementation at https://github.com/epoch-research/finite_field_implementation/},\mathbb{F}_{125

Figure fig:proj-solution-counter: Python implementation of a projective solution counter, used to identify the pattern in solution counts over finite fields $F_5$ , $F_25$ , and $F_125.$

We can also verify this curve is smooth at $p=5$ and has genus $3$ .

The Weil conjectures maxhasseweil gives us [|C[F_p^n]|=p^n - _i=1^2g_i^n+1] for some complex $\alpha_i$ of absolute value $\sqrt{p}$ which come in conjugate pairs. calebji

We notice

[_j = e^(2j+1)i/65]

satisfies this equation.

To show that this is unique, let $\alpha_i$ be roots of $z^2 - \sqrt{5}\beta_iz + 5$ , $i=1,2,3$ . Then $\alpha_i^2$ are roots od $z^2 - 5T_2(\beta_i)z + 25$ and $\alpha_i^3$ are roots of $z^2 - \sqrt{125}T_3(\beta_i)z + 215$ , where $T_2$ and $T_3$ are the 2nd and 3rd Chebyshev polynomials. So we know [_1+_2+_3= T_2(_1)+T_2(_2)+T_2(_3)=T_3(_1)+T_3(_2)+T_3(_3)=0.]

So, knowing the values of $\beta_1+\beta_2+\beta_3,\beta_1^2+\beta_2^2+\beta_3^2$ and $\beta_1^3+\beta_2^3+\beta_3^3$ gives you the $b_i$ (and so the ( $\alpha_i$ ) uniquely).

Setting $n=18$ gives us the result:

[|C[F_5^18]|=5^18 + 6 5^9 + 1. ]

Evaluation

Prompts

To guide the models effectively in solving the complex mathematical problems, we wrote prompts that encourage mathematical reasoning, writing code for experiments, and adherence to formatting requirements.

Initial Prompt

The initial prompt introduces the problem and provides detailed instructions. It explains how the model can submit Python code enclosed in triple backticks (|```|) to run experiments. It explains how to submit a final answer: the submitted code should contain the exact comment |# This is the final answer|, and pickle dump the answer to a file named |final_answer.p|. Finally, the prompt reminds the model that it cannot natively perform arithmetic and algebra, and therefore encourages it to use Python code to verify this reasoning.

Continuation Prompts

After each interaction, if the model has not yet arrived at a correct solution, we provide feedback based on the code execution results. The continuation prompt first includes any error messages, standard output, or indications of timeouts from the executed code. It then urges the model to keep working and try new approaches if stuck, and reminds it of the guidelines for submitting the final answer.

Final Prompt

If the model has not submitted a correct answer within token limit, we provide a final prompt that informs the model that the answer is needed now and instructs it to submit the best attempt. This final prompt reiterates the requirements for the final answer code block.

For problems where authors expect a solution to take more than 20 seconds to run, we removed this.

Analysis of problems solved at least once

We selected four problems that could be solved by at least one model, and conducted repeated trials with five runs per model per problem. We report the success rate over the five trials of each model on each problem in Table tab:model-comparison (naming problems from 1 to 4).

Table tab:model-comparison: Model success rates out of five runs on the four test problems.

\begin{table}[h]
\centering
\begin{small}
\begin{tabular*}{\textwidth}{@{\extracolsep{\fill}}l*{6}{c}l@{}}
\toprule
\textbf{Problem} &
\textbf{\makecell{Grok 2\\Beta}} &
\textbf{o1-preview} &
\textbf{o1-mini} &
\textbf{GPT-4} &
\textbf{\makecell{Gemini 1.5\\Pro}} &
\textbf{\makecell{Claude 3.5\\Sonnet}} &
\textbf{MSC2020 Classification} \\
\midrule
1 & \textbf{60\%} & \textbf{40\%} & \textbf{20\%} & 0\% & 0\% & \textbf{20\%} &
\makecell[l]{60 Probability theory\\41 Approximations and expansions} \\[1ex]
\midrule
2 & 0\% & \textbf{100\%} & \textbf{20\%} & \textbf{20\%} & 0\% & \textbf{80\%} &
\makecell[l]{55 Algebraic topology\\57 Manifolds and cell complexes} \\[1ex]
\midrule
3 & 0\% & \textbf{80\%} & 0\% & 0\% & \textbf{60\%} & 0\% &
\makecell[l]{20 Group theory\\12 Field theory} \\[1ex]
\midrule
4 & 0\% & 0\% & \textbf{20\%} & 0\% & 0\% & 0\% &
\makecell[l]{14 Algebraic geometry\\05 Combinatorics} \\
\bottomrule
\end{tabular*}
\vspace{0.5em}
\caption{Model success rates out of five runs on the four test problems.}
\label{tab:model-comparison}
\end{small}
\end{table}

Percentage of unsolved problems across benchmarks

Figure fig:benchmark_comparison shows the fraction of problems not solved by the top-performing model on prominent mathematics benchmarks: the College Mathematics subset of MMLU ([hendrycks2021measuringmassivemultitasklanguage]), GSM8K ([cobbe2021training]), MATH ([hendrycks2021measuring]), AIME ([o1-learningtoreason]), Omni-MATH ([gao2024omni]), MathVista ([lu2024mathvistaevaluatingmathematicalreasoning]), and FrontierMath.

To make this figure, we took one minus the top reported accuracy for each benchmark under the pass@1 setting. For MMLU, we used the o1 accuracy of 98.1% ([o1-learningtoreason]). For GSM8K, we used the Claude 3.5 Sonnet pass@1 accuracy of 96.4% ([sonnet-3-5-addendum]). For MATH, we used o1's reported pass@1 accuracy of 94.8% ([o1-learningtoreason]). For AIME, we used o1's reported pass@1 accuracy of 74% ([o1-learningtoreason]). For Omni-MATH, we took o1-mini's accuracy of 60.54% as reported on the benchmark leaderboard [gao2024omni]. For MathVista, we used o1's pass@1 accuracy of 73.9% ([o1-learningtoreason]). Finally, for FrontierMath, we relied on our own evaluation, which shows that more than 98% of the problems remain unsolved by the best-performing models.

Pass@8 accuracy

Figure fig:model-performance illustrates the mean accuracy across 8 different evaluations, where no model achieves higher than 2%. Figure fig:pass@8 we show the pass@8 accuracy of different models on FrontierMath.

Pass@8 accuracy of different models on the FrontierMath dataset. The best model (o1-preview) achieves around 6% accuracy, compared to around 2% for Grok 2 Beta.

Figure fig:pass@8: Pass@8 accuracy of different models on the FrontierMath dataset. The best model (o1-preview) achieves around 6% accuracy, compared to around 2% for Grok 2 Beta.

Difficulty ratings

Our approach to difficulty rating involves a combination of expert self-assessment and rigorous peer review to ensure the reliability and consistency of the ratings. Each problem author provided an initial difficulty assessment along three key dimensions:

enumerate Background: This rating ranges from 1 to 5 and indicates the level of mathematical background required to approach the problem. A rating of 1 corresponds to high school level, 2 to early undergraduate level, 3 to late undergraduate level, 4 to graduate level, and 5 to research level. Creativity: Estimated as the number of hours an expert would need to identify the key ideas for solving the problem. This measure has no upper limit. Execution: Similarly estimated as the number of hours required to compute the final answer once the key ideas are found, also with no upper limit. enumerate

To verify and refine the initial difficulty assessments, each problem underwent a peer-review process. Reviewers, who are also mathematicians, examined the problems along with the provided solutions and difficulty ratings. They assessed the correctness of the solutions, the clarity of the problem statements, and the accuracy of the initial difficulty ratings. Any discrepancies between the authors' ratings and the reviewers' assessments were discussed, and adjustments were made as needed to ensure that the difficulty ratings accurately reflected the problem's demands in terms of background knowledge, creativity, and execution.

This peer-review process not only validated the difficulty ratings but also enhanced the overall quality of the problem set by ensuring correctness and clarity.

Our initial analyses have shown that substantially easier problems that are not in this benchmark that have received higher difficulty ratings are significantly less likely to be solved by advanced AI models like GPT-4o, supporting the validity of our difficulty assessment system.

Metadata collected

This is the metadata we collect on each problem. Writers put in their own suggestions for this metadata when submitting a problem, but this gets edited by the reviewers to ensure consistency of difficulty ratings and of the relative granularity of subject vs. technique classification.

Table table:metadata: Table of metadata tags and their descriptions

\begin{table}[h!]
    \small
    \centering
    \begin{tabular}{>{\raggedright\arraybackslash}p{3.5cm} p{10cm}}
        \toprule
        \textbf{Metadata tag} & \textbf{Description} \\
        \midrule
        \textbf{Background rating} & A difficulty rating ranging from 1 to 5 quantifying how much background knowledge is required to solve the problem. \\
        & \begin{tabular}{@{}l@{}}
            1: High school level \\
            2: Early undergraduate level \\
            3: Late undergraduate level \\
            4: Graduate level \\
            5: Research level
        \end{tabular} \\
        \midrule
        \textbf{Creativity rating} & This rating is based on how long we think the average human expert with the requisite background knowledge would take to find the key ideas going into the solution of a problem. \\
        \midrule
        \textbf{Execution rating} & Measures the amount of attention to detail and precise reasoning required to solve the problem. \\
        \midrule
        \textbf{Subjects} & A list of broad subjects that the question fits into: “analytic number theory”, “representation theory”, “differential geometry”, etc. \\
        & \begin{tabular}{@{}l@{}}
            Note that a question can fit into more than one subject.
        \end{tabular} \\
        \midrule
        \textbf{Techniques} & A list of techniques, theorems, results, etc. that can be used to solve the problem. The list doesn’t need to be exhaustive and only needs to contain the “most prominent” items that occur to the question author. \\
        & \begin{tabular}{@{}l@{}}
            Example: “generating functions”, “double counting”, “Vieta’s theorem”,\\ etc.
        \end{tabular} \\
        \midrule
        \textbf{Is programming required?} & This is either “yes” or “no”, depending on whether a human expert would require access to a programming environment in order to find the answer to the question. \\
        \bottomrule
    \end{tabular}
    \caption{Table of metadata tags and their descriptions}
    \label{table:metadata}
\end{table}