# Questions tagged [number-theory]

Questions on the number-theoretic functionality of Mathematica.

319 questions
Filter by
Sorted by
Tagged with
94 views

### My version of PowerMod breaks down around 10^308

I have been trying to write a function that duplicates PowerMod[a, b, n], computing a^b mod n...
195 views

### Table of Chebyshev psi function

This should be easy, but for some reason I'm struggling. The Chebyshev psi function is given as the sum from 0 to x of the von ...
95 views

### List of invertible congruence classes

I am attempting to create a list of the invertible congruence classes $\bmod 120$. The code I have is ...
63 views

### Using a Do-loop to find divisors mod 13 [closed]

I want to check sum of divisors of i mod 13 fori = 1 to i = 20. I tried writing a Do-Print ...
197 views

### Mathematica code for computing the $p$-adic expansion of rational numbers

Does anyone know any Mathematica code for computing the $p$-adic expansion of rational numbers? I.e. given a rational number $a/b,~a,b\in \mathbb{Z}$ and a prime number $p$, then compute the $p$-adic ...
139 views

### Computing the seven roots of a polynomial

This question was originally asked by @fsrong70 six months ago. The OP deleted it shortly after posting and has not returned to this site since. I had just figured it out when it was deleted. I ...
129 views

### Not comprehending PowerMod

Bug introduced in 10.3 or earlier and persisting through 11.3.0 or later The bug is not present in 10.0. By definition, PowerMod[a, 1/r, m] finds a modular rth root of a mod m. Here's a pair of ...
842 views

### Determining if a number is divisible by 1000 [closed]

I have a number such as: a = 875952; And I want to find if it is divisible by 1000. Is there a concise way of doing that?
461 views

### Generate numbers relatively prime with a given number

I am interested in a function such that f[m, i] = n where m, n are positive integers and n ...
77 views

### Find the maximum of a function involved with Floor function

My function is $$f(H,p) = \left\lfloor \dfrac{\lfloor H/p\rfloor + 3 - \sqrt{(\lfloor H/p\rfloor + 1)^2 - 4H}}{2} \right\rfloor$$ The constraints are $H \geq p(4p-1)$, $p$ is prime although ...
74 views

### What is the best build-in function to pick or select given element in a huge list

im looking for a best and efficient function works like a search engine it takes for example m=5+6I ,then it goes searching in the list V={1,1+I,2+3I,...} until catch it . My Dr said to use "Select[]"...
349 views

### Efficient code for minimum integer with given number of factors

I'm seeking an efficient implementation of the number-theoretic function giving the smallest integer $n$ that has exactly $k$ factors (not necessarily prime): ...
152 views

### Finding $x$, the exponent of a $2^x$ when we need an specific output

Say we need a program Findx[n_Integer, m_Integer] where n is an integer from 1 to 9 and m is an integer from 1 to 1000. The output of the program is a number x, which is the exponent of 2 that ...
157 views

### Solving Diophantine equations 5 [closed]

Given a positive integers x,y,m would like to be able to find integer solutions z from Diophantine equation x^2-y^2 = m*z in Z.
83 views

### Problem in Counting the Number of Divisors with a Function [closed]

I want to create a function in mathematica which returns the number of divisors of a parameter x. So, I created the following function: d[x_] = Length[Divisors[x]]; And when I evaluate: d d ...
28 views

### Finding smallest domain within which variables can satisfy inequality

Given an arbitrary number of variables $\epsilon_i$ that can be picked from a domain $[0, W]$ and some inequality relation between all the variables $G(\epsilon_1, \epsilon_2, ...)$, is there some way ...
520 views

### Does Mathematica have a twin prime equivalent of PrimePi?

Well, it's all there in the title! I'd like to be able to plot the number of twin primes =<x. Is there an inbuilt function that can do this?
109 views

### How to find, load and use a legacy Mathematica package

Question: How to find, load and use a legacy Mathematica package from a previous Mathematica version? What is the preferred approach and process? Context: I am working through a set of examples in ...
41 views

### FunctionDomain in the Reals numbers [closed]

Why don't I obtain that $x$ belongs to the Reals? ...
93 views

### Simultaneous equation [closed]

How to solve Simultaneous equation in Mathematica: x^2 - (a^2 - 1)y^2 = 1 and y^2 - pz^2 = 1 where ...
96 views

### Number-theoretic function using Table

I'm trying to make a function to calculate the Log sum of primes over a limited range $1/2n$ to $n$ or Chebychev theta function over limited range $1/2n$ to $n$. This will be used only for even ...
351 views

### Find the order $m$ of a matrix ${\bf A}$ such that ${\bf A}^m= {\bf 1}$

I have a square matrix ${\bf A}$ defined over the field $\mathbb Z_2$ and I want to find its order such that ${\bf A}^m=1$. I tried using ...
61 views

### How to define a function the Dirichlet L-function $L(s,\overline{\chi(5,2)})$ in Mathematica?

In Mathematica: The Dirichlet L-function with character $\chi(5,2)$, $L(s,\chi(5,2))$, is expressed as DirichletL[5,2,s] Let $\overline{\chi(5,2)}$ be the complex ...
1k views

### François Viète's approximation to π

How do I program the approximation to π devised by François Viète, which is given by 2 * 2/Sqrt * 2/Sqrt[2 + Sqrt] * 2/Sqrt[2 + Sqrt[2 + Sqrt]] * ... ...
49 views

### Question about using ParallelMap to speed up computations [duplicate]

First; I define two functions: ...
110 views

### Code for finding $a$ and $b$ such that $a b = 1 \mod 4$ [closed]

I need to find $a$ and $b$ such that $a b = 1 \mod 4$? I do not know how to write the code. Could someone help me?
129 views

### Generating numbers palindromic in two number bases

The purpose of the code below is to generate numbers that are $2d+1$ digit palindromes in number base $b+1$, and are also palindromic in number base $b+3$, where: The ...
186 views

### Easy number theory problem

$p$ is an odd prime number，$S = \{x \mid 1 \leq x \leq 2p, x \in \mathbb{Z} \}$, $A$ is a subset of $S$, satisfying $\operatorname{card}(A) = p$ $\sum\limits_{x\in A} x \equiv 0 \pmod p$...
125 views

### Iterative calculation of a number-theoretical constant [closed]

Recently the decimal expansion of a number theoretic constant was searched for, which is the analog of the Landau-Ramanujan constant in a certain context. The constant starts 0.638909... It can be ...
99 views

### Finding minimum x such that Mod[3^x, m] == 1 for m not multiple of 3 [closed]

I would like to find smallest x value for each m value such that Mod[3^x, m] == 1, where m is not multiple of 3. Here is my ...
65 views

### Listing divisors of a number [closed]

I created a list with all divisors of 18000: list1 = Divisors I did an analysis to identify which of these numbers are divisible by $15$, but the result I ...
246 views

### Question about how to speed up Mathematica code

When looking at the Minimal Goldbach prime partition point {p,q} for each n; where n=10^i and i = 2,3,4,...,10; I notice that these points reside in an interval with center n/2 and radius 250: ...
206 views

561 views

### InverseTotient[ ]?

Maple has a function InverseTotient( c ), which returns all those natural numbers $n$ whose Euler totient function $\phi( n ) = c$. Is there an equivalent inverse ...
321 views

### How to represent integers using Egyptian fractions?

Let $N(n)$ be a set of integers, which can be presented using first $n$ Egyptian fractions: $$N(n):=\{m\in\mathbb{Z}:\ \ m=\sum_{i=1}^n\frac{\epsilon_i}{i},\ \epsilon_i=0\ \text{or}\ 1\}$$ I want ...
111 views

### Solving a modular equation with a large modulus [duplicate]

I am having trouble solving this equation: $17^n\equiv 1 \mod 640826899755722841329632946651705039010285107743541637238400825304943424424715705661705191669759$ for the smallest integer $n > 0$ I ...
210 views

### Generation of Step Numbers

I am working on Project Euler 178 but got stuck in trying to optimize my code. The following text comes from the problem: Consider the number $45656$. It can be seen that each pair of ...
126 views

### Calculate 40 digits of the MRB constant

MRB constant is the upper limit point of the following sequence $$s_n=\sum_{k=1}^{n} (-1)^k k^{\frac{1}{k}}$$ $MRB=\color{blue}{0.1878596}...$ I tried to calculate first few digits: ...
56 views

### Question about how to use NestWhileList

I start with: n = 2228; m = n/2; PreviousPrime[n_] := NextPrime[n, -1] I use NestList to build the following list: ...
142 views

### How can I find the $c$ such $Max[Fibonacci[Range[c]]] = 13$?

How can I find the $c$ such Max[Fibonacci[Range[c]]] = 13 I tried Reduce but there is an error message ...
169 views

### Number of primes between two integers x and y (with x < y and excluding x and y)

There is a formula given at the bottom of the following webpage: https://math.stackexchange.com/questions/288747/how-to-find-number-of-prime-numbers-between-two-integers to calculate the number of ...
143 views

### Listing products of prime powers

Given a positive integer $n$, what is the code to list $2^{a_2}3^{a_3}\cdots p^{a_p}$, where $a_i\ge 0$ are integers, with respect to the lexicographic ordering on $(a_2,a_3,\ldots, a_p)$? The only ...
### Finding all integer solutions of the following inequality $\bigg| \sqrt{2}-\frac{p}{q} \bigg | <\frac{1}{q^{5/2}}$
I want to find integer solutions of the following inequality by using Mathematica $$\bigg| \sqrt{2}-\frac{p}{q} \bigg | <\frac{1}{q^{5/2}}$$ ...