Introduce the nonnegative integer surplus variables and . Solving these two equalities gives the optimal relaxed simplex dictionary
with objective
The printed tableau has a sign error in the first row: in equality-row form the coefficient of must be , rather than . The defining constraints and the printed objective row both confirm this correction. The second row, from which the required cut is obtained, is unaffected.
The cutting-plane method solves the linear programming relaxation and, if the optimum is fractional, adds an inequality valid for every feasible integer point but violated by that optimum. For an all-integer tableau row
write each coefficient as its integer part plus its fractional part. The quantity is an integer, while is nonnegative. Thus the integer is at most , and subtraction gives the Gomory fractional cut
It retains every feasible integer point and cuts off the fractional basic point with all nonbasic .
Apply this to the row . Since , and , the resulting valid inequality is
Multiply by five and add an integer slack . Its row is . At the old relaxed optimum, is negative, while the minimization objective coefficients remain nonnegative. This is the starting configuration for the dual simplex algorithm: retain dual feasibility and pivot to repair the negative basic value.
The cut row permits either or to enter. The dual ratio test gives
Thus enters and leaves. Solving the cut row and substituting gives
Setting the nonbasic variables gives nonnegative basic values . The new objective coefficients are nonnegative, so this is optimal for the cut relaxation. Its decision variables are already integer and satisfy the original inequalities. Therefore
Indeed the positive coefficients of both nonbasic variables make this relaxed optimizer unique. No additional cut is needed.
Let be an extreme point and consider its active constraints. Their row vectors span . Otherwise a nonzero vector is orthogonal to every active row. All inactive inequalities have strictly positive slack; since there are finitely many, a sufficiently small makes both and feasible. Their midpoint is , contradicting extremality.
Choose independent active rows and write the corresponding equalities as . The matrix has integer entries, so is a nonzero integer and . By Cramer's rule,
where replaces column by . Every entry of has absolute value at most , even when is real. The determinant expansion has terms, each with absolute value at most . Therefore
This proves the integer-matrix vertex coordinate bound; only the denominator determinant requires integrality.
For an exact finite-input ellipsoid method, the numerical data must have a specified encoding. In the standard model the input is the dimensions, the integer matrix and a rational right-hand side in binary, together with bounds derived from their encoding lengths. The vertex estimate above remains valid for arbitrary real , but arbitrary unencoded reals do not by themselves specify an exact computational decision problem. Clear the rational row denominators to obtain an equivalent integer system, and let bound the resulting integer coefficients and right-hand sides.
Two precautions make a feasibility account valid even for unbounded or lower-dimensional polyhedra. First, nonempty polyhedra need not have extreme points. Nevertheless they have a feasible point in the box with . To prove this feasible-point box bound for a polyhedron with lines, intersect the polyhedron with an orthant containing one feasible point. Minimize the sum of the signed coordinates in that orthant. Its nonempty sublevel sets are closed and bounded, so the minimum is attained. The compact minimum face has an extreme point, which is also extreme in the orthant intersection. Its active rows come from the integer system and signed coordinate inequalities, still bounded by . Applying the proved vertex estimate gives the box bound. This is an existence argument; the algorithm need not enumerate all orthants.
Second, a nonempty feasible set can have zero ambient volume. Use a small uniform relaxation with a certified rational infeasibility gap, rather than stopping on volume alone. For the integer system , put
The normalized integer Farkas infeasibility gap says that if the system is empty there is a multiplier with , and . Its bound follows by maximizing the certificate's value in the normalized multiplier polytope: a vertex has at most positive entries, and a determinant denominator bounded by . Any point satisfying the relaxed inequalities would then give the contradiction
Thus relaxation preserves emptiness. If the original system is nonempty, take its point . The bounded relaxed set
contains a Euclidean ball around of radius : each row changes by at most , and the expanded box leaves room for the ball. Hence is either empty or has an explicit positive inner-volume bound. This is the full-dimensional relaxation of integer inequalities principle.
Start with a ball, viewed as an ellipsoid, centered at zero with radius , which contains the search box. At each iteration test its center against all inequalities defining . If the center is feasible, the preserved-emptiness result answers that the original problem is nonempty. Otherwise a violated row supplies a separation oracle cut. Move that cut parallel to itself through the current center if necessary, and compute a smaller containing ellipsoid for the surviving half. Every feasible point remains in the new ellipsoid.
The central-cut ellipsoid volume bound gives a shrinkage factor at most for . After cuts, the volume is below that of a radius- ball. If no feasible center has been found, and therefore the original system are empty. Dimension one uses interval bisection. The logarithms of , and are polynomial in the encoded input size. Standard outward-controlled finite-precision ellipsoid updates preserve enclosure with polynomial bit complexity.
The coordinate bound therefore supplies the initial bounded search region; the rational gap supplies a valid stopping threshold. Neither unboundedness nor zero volume of the original polyhedron can be ignored when deciding emptiness.
For the metric travelling salesman problem, a tour is a polynomial-size complexity certificate: list the vertices in cyclic visiting order. In polynomial time verify that every vertex occurs exactly once, add the encoded integer edge costs and compare the sum with the threshold. The certificate has bits, and addition uses a number of bits polynomial in the input length. If validity of the metric input is not treated as a promise, nonnegativity, symmetry and all triangle inequalities can also be checked in polynomial time. Therefore the decision problem belongs to NP.
An evaluation algorithm immediately decides a threshold instance by computing the optimum and comparing it with . Conversely, suppose there is a polynomial-time decision algorithm. The optimum is an integer in , where is the largest edge cost; any fixed cyclic ordering gives the upper bound. Binary search on this interval using threshold queries finds the exact minimum. It makes queries, and every queried threshold has polynomial encoding length. Thus the time is polynomial in the bit length, even when is numerically large. Zero-cost and trivial small instances can be handled directly. Consequently
This concerns the optimum value; constructing an optimal tour is the separate optimization task.
For the hardness reduction, start with an undirected graph on vertices and form the complete graph with costs
These costs satisfy the triangle inequality, since every direct cost is at most two and every two-edge path costs at least two. Set the threshold to . Every tour has exactly edges, each costing at least one. Hence it has cost at most precisely when all its edges came from , which is precisely the existence of a Hamiltonian cycle in . The construction and threshold encoding are polynomial-time. Thus the assumed NP-completeness of the Hamiltonian circuit problem gives NP-hardness of metric TSP; together with its NP membership this proves NP-completeness of metric TSP decision.
For the approximation, compute a minimum spanning tree . Deleting any one edge of an optimal tour leaves a spanning tree, so nonnegative costs imply
Double every tree edge. The resulting connected multigraph has even degrees and hence an Euler circuit, of cost . Traverse it and skip previously visited vertices, except for the final return to the start. Each skip replaces a path by its direct edge; repeated use of the triangle inequality shows that no replacement increases cost. The result is a tour visiting every vertex exactly once and satisfying
The Euler traversal and shortcutting take polynomial time, so the assumed polynomial-time MST routine supplies the requested polynomial-time algorithm. This is the double-tree approximation for metric TSP. In the relative-error convention, when the optimum is positive, explaining the printed term “1-approximation”; in the multiplicative convention it is a factor-two approximation. If the optimum is zero, the same bound forces the returned tour to have zero cost as well.
Write for a balanced complete multipartite blow-up with classes of size . In the nonvacuous range , the answer is
We prove both bounds, including the quantitative lower bound.
First, suppose a bipartite graph with parts , of sizes , has at least edges. For any integer , counting pairs consisting of an -subset of and a common neighbour in gives
The inequality follows from discrete convexity: is nondecreasing, so balancing integer degrees of vertices minimizes their sum. Consequently some has
The product estimate uses for . This is a common neighbourhood from bipartite density bound.
We next prove by induction on that, for every fixed , a graph of order with minimum degree at least contains , for some and all sufficiently large . Impossible minimum-vertex degree parameters need no argument. For , take two copies of with adjacency given by . The preceding count with and , where , gives at least common neighbours. The selected vertices and their common neighbours are disjoint in , since a vertex is not its own neighbour. Selecting common neighbours gives .
For , the minimum degree hypothesis implies the induction hypothesis for , for example with excess . Thus contains with for a positive constant . Pair the vertices in its classes into disjoint transversal cliques , each of order . By the union bound, every has at least
common neighbours. Form an auxiliary bipartite graph between these cliques and the vertices, using common-neighbour incidence, and put . Choose small enough that and . The count above, with , gives cliques with at least simultaneous common neighbours. Their vertices form classes of size , and of those common neighbours form the last class. None of the common neighbours belongs to a selected clique, because that would require a loop. This proves the logarithmic clique blow-up from minimum degree assertion.
To pass from an edge count to minimum degree, put and . Repeatedly remove a vertex whose current degree of a vertex is below times the current order. If the process reaches order , its deleted edges and remaining edges give
Thus . In particular, for a fixed sufficiently small , the process cannot reach : the displayed inequality would fail. It stops earlier with a subgraph of order and minimum degree at least . The induction just proved supplies , hence a lower bound on . In particular, the guaranteed part size tends to infinity.
For the upper bound, choose a fixed with and take a binomial random graph . Its edge count has expectation and variance . The Chebyshev inequality shows that with probability tending to one. The expected number of labeled copies of is at most
For and , this tends to zero. The Markov inequality therefore shows that with positive probability both the required edge density holds and there is no such copy. This proves and completes the random obstruction to larger clique blow-ups argument.
There is a necessary qualification to the density parameter. If , the only eligible graph is the complete graph, and . If and , there are no eligible simple graphs; the universal assertion is vacuous for every , so the printed maximum has no finite value in that range.
We use the following precise form of the Szemerédi regularity lemma. For every and positive integer , there are integers such that every graph on vertices has a partition
with at most irregular unordered pairs among the nonexceptional classes. Here for disjoint nonempty sets, and is an -regular pair if
whenever , , , .
For the proof, let be the adjacency indicator, including , and define by averaging on even for diagonal cells. The regularity energy of a partition is
Every exceptional vertex is represented by its own singleton cell in this energy. If a rectangle is refined into smaller rectangles with densities and weights , expanding the square gives the refinement variance identity for regularity energy
Thus refinement never decreases energy.
Suppose an equitable partition with classes of size fails the desired conclusion. For every irregular pair choose witness subsets , , each of size at least , whose density differs from by more than . Split each class by all witness subsets involving it. Each class has at most atoms. On a witness rectangle, the weighted mean of is . The Cauchy-Schwarz inequality therefore bounds the energy increment on this rectangle below by
The different original cell pairs have disjoint rectangles. More than irregular unordered pairs, with both orientations counted in the energy, thus increase it by more than . As long as , this is greater than , and in particular greater than .
Restore equitability by dividing each atom into pieces of common size , moving each leftover vertex into the exceptional set as a singleton energy cell. This remains a refinement, so the energy gain persists. There are at most atoms, so the additional exceptional vertices number at most
If , then , and the new number of equal cells satisfies . For the lower bound, within each old class the leftover count is at most , so at least one full new cell survives. For the upper bound, , using . This proves equalization with a controlled exceptional set without any energy loss.
Put . Choose with . Start with equal classes and fewer than exceptional vertices. Iterate the bounded recurrence times to obtain a bound . Choose so large that the initial remainder is at most and, for every , . The total exceptional set remains at most , since at most rounds each add at most . Each failing round raises an energy in by more than , so fewer than rounds can fail. This proves the lemma with all constants independent of the input graph.
For the application, write and . We may replace the given tolerance by a smaller one, so assume . Set , take , and choose
Apply the lemma and delete all edges incident with , all edges inside a class, all edges across irregular pairs, and all edges across regular pairs of density below . The resulting spanning subgraph loses at most
Let be the reduced graph of a regularity partition, whose edges represent the surviving regular pairs of density at least . We show that a copy of in forces a copy of in .
Choose a proper -coloring of and assign its color classes to the corresponding clusters. Embed its vertices one at a time. For each unembedded vertex maintain the candidate set in its cluster consisting of unused vertices adjacent to all already embedded neighbours. After at most restrictions its size is at least . This last quantity is at least when is large, so every such candidate set is large enough for regularity. In an -regular pair of density , fewer than vertices on one side have fewer than neighbours in any fixed set on the other side of size at least : otherwise these vertices and violate regularity. For the vertex currently being embedded, there are at most relevant future candidate sets. At most vertices are bad for one of them. Its own candidate set has at least vertices, so an unused good vertex can be chosen. All future candidate sets shrink by at most the factor , apart from the removal of a used vertex; the stated size bound follows inductively, since the accumulated removals are at most . This completes the graph embedding lemma for regular pairs in the form needed here.
Since is -free, is -free. Any in would use distinct clusters, because intraclass edges were deleted, and project to a in . Therefore
The required threshold on depends only on and , since makes every nonexceptional cluster sufficiently large.
An integer square is divisible by four exactly when its base is even: an even base has square , while an odd base has square . Apply this observation to the triangular number obtained above. Since is even,
Checking the four residue classes of modulo four, the products have residues , respectively. Therefore
The original PDF has two alternatives modulo four; the TeX transcription incorrectly turns the second alternative into a modulus-three expression.
This is not an equivalence relation. For , , so reflexivity fails. It is symmetric, because is the conjugate of , and an integer is real. It also fails transitivity: take and . Then are integers but is not. A single failure would suffice, but these examples identify exactly which properties break.
This is an equivalence relation. The difference is an even integer, proving reflexivity. If then , proving that the relation is symmetric. If and then , proving transitivity. Its equivalence classes are the translates ; two real numbers need not be integers themselves to belong to the same class.
For integers not both zero, their greatest common divisor is the greatest positive integer dividing both. For the degenerate pair use the conventional value zero. A number divides both exactly when it divides both : one direction uses closure of divisibility under sums; the reverse subtracts . The positive common-divisor sets are therefore identical, proving
This also holds for under the stated convention and does not depend on the signs of the integers.
Bezout identity states that there exist integers such that . If a prime divides but does not divide , then , so choose with . Multiplying by gives . Both terms on the left are divisible by , so . Thus
which proves the Euclid lemma from the requested identity. The separately numbered entries use these arithmetic facts for the Fibonacci numbers.
Write for the Chebyshev theta function and for the Second Chebyshev function. Throughout, denotes a prime and all logarithms are natural. Each prime divides the central binomial coefficient, so
Apply this at and telescope: . Bounding a general by the next power of two proves . Splitting the prime-counting function at now gives
This is the required Chebyshev estimate from central binomial coefficients, with no use of a prime asymptotic.
The higher prime powers also satisfy
For an integer , the prime factorization of a factorial gives
The Stirling formula, or integral comparison for , says . After division by we obtain . Moreover
Subtract this bounded contribution and replace by . The factorial proof of Mertens first theorem therefore yields
For fixed , apply partial summation to and :
The main contribution is
The bounded error in contributes at the endpoint and in the integral, since . Thus the positive logarithmic moments of reciprocal primes satisfy
Use the Jacobi theta function , . Its transformation law follows directly from the Gaussian Fourier transform. With the Fourier transform convention ,
For example, differentiating the transform and using integration by parts gives ; its value at zero is the Gaussian integral . Periodize the Gaussian. Its Fourier coefficients are these transforms at integer frequencies, and both series converge absolutely, so evaluating its Fourier series at zero proves the Poisson summation formula here:
For , integrating the exponentially convergent series termwise gives the Mellin representation of the completed Riemann zeta function
Split at one and substitute in the lower part, using the transformation just proved. The contribution of is . We obtain the pole-subtracted theta integral for the completed zeta function
The integral is entire: decays exponentially, uniformly dominating the integrand and all its derivatives on compact sets. Both the rational term and the integral are invariant under . This proves the meromorphic identity
Equivalently the Riemann xi function
is entire and satisfies . The reflection and duplication formulas for the Gamma function turn the same identity into
understood by meromorphic continuation at removable exceptional points. This establishes the functional equation of the Riemann zeta function, rather than assuming it.
If counts zeros in the critical strip with , including zero multiplicity, the Riemann–von Mangoldt formula is
At a zero ordinate one may instead use the conventional half-weight boundary count; that only alters the stated error. In particular .
To invert without assuming simple zeros, fix and set . Then , so
For large , , and hence . Let . The asymptotic inversion of the zeta zero count gives
As in the zero-counting formula, the sequence lists ordinates with zero multiplicity; “increasing” means nondecreasing here. The argument does not require the Riemann hypothesis or simplicity of zeros.
The Selberg upper-bound sieve replaces an exact but long inclusion-exclusion expansion by a short nonnegative square, and then minimizes its main term. For an interval , the count of multiples of every integer is
Crucially the remainder bound is independent of .
Let and take real weights , supported on squarefree , with . If , the sum of these weights over divisors of is just one; for all other its square is nonnegative. Thus
Expanding and counting common multiples gives
This separates the local density from a uniform interval error. The weights will be chosen to minimize .
The Euler totient function identity yields the Selberg sieve diagonalization
Only squarefree occur. Inverting this finite triangular sum by Möbius inversion gives
Write . The Cauchy-Schwarz inequality gives , with equality at . The inverse formula then produces the optimal Selberg sieve weights
The constraints are satisfied, since the diagonal transformation is invertible and .
We need the leading constant, not just . Here is an elementary proof of the Selberg sieve denominator asymptotic. For define the multiplicative arithmetic function by
Checking the prime-local factors shows . Also : the local factors for absolute values differ from one by , and their logarithmic first moments sum . The signed local sums are all one, so absolute convergence gives . Using the harmonic number estimate,
For the last equality, is bounded and is bounded too.
The optimal Selberg weights have modulus at most one. To prove this here, for squarefree use . The product is the sum of over , , , with squarefree. All are distinct squarefree integers at most , and so these are a subset of the terms of . The displayed weight formula now gives . Consequently
Every prime in the interval exceeding survives the sieve; at most smaller primes must be restored. Choose . Then
and the uniform interval bound from optimal Selberg weights follows:
All constants are independent of . Given , make the absolute error in the coefficient smaller than by choosing large enough. Uniformly for every and ,
Fix a modulus and a residue class with . The Dirichlet characters modulo are the one-dimensional characters of a representation of , extended by zero on nonunits. Let be the principal Dirichlet character.
For every nonprincipal Dirichlet character the partial sums are bounded: the sum over a full period is zero, and a remaining incomplete period has bounded length. Partial summation therefore continues its Dirichlet L-function holomorphically to by
This shows in particular that it has no pole at one. By contrast
has a simple pole there. We use the stipulated nonvanishing for every real nonprincipal Dirichlet character, and establish the remaining nonvanishing as follows.
For real , take the product over all Dirichlet characters. The Euler product logarithms and Orthogonality of Dirichlet characters give
Thus this product is real and at least one. If a nonreal Dirichlet character had , complex conjugation would also give . These are distinct factors. Two zeros and only one possible simple pole force the whole product to tend to zero as , contradicting the lower bound. Hence all nonprincipal Dirichlet characters have , under the given real-character assumption.
For , absolute convergence of each Euler product gives
with bounded remainder as : the terms of prime-power exponent at least two have total modulus at most . Since a nonprincipal defines a holomorphic function that is nonzero near one, its Euler-product logarithm differs from a local holomorphic logarithm by a fixed integer multiple of on a short real interval. It is therefore bounded there. For the principal Dirichlet character the pole instead gives
Apply the Dirichlet character orthogonality relation at each prime. The prime-character sum near one yields
This tends to infinity. A finite set of primes would have a bounded sum as , so every coprime residue class contains infinitely many primes. This proves the Dirichlet theorem on primes in arithmetic progressions from the specified assumption. In fact it proves divergence of the corresponding reciprocal-prime series and gives Dirichlet density .
Finally, for real each prime satisfies
The integrands are nonnegative, so Tonelli theorem permits summation inside the integral. At a fixed , exactly primes from the residue class have entered. Therefore
Both sides are finite, since they are bounded respectively by and . The values of a counting function at its isolated jump endpoints do not change this integral.
For a continuous map , use the classes from part (a). Since , there is an integer with . Naturality of the cup product and yield
With the product orientation on the source and the complex orientation on the target, the mapping degree is defined by . Evaluation on the fundamental class gives
Thus the induced group homomorphism on is zero. This is a cup-power obstruction to nonzero degree, rather than a consequence of the additive homology groups alone.
A linear polyhedron is full-dimensional if its affine hull is . For a convex set, this is equivalent to containing an open Euclidean ball: an open ball has full affine span, while affinely independent points in the set span a simplex with nonempty interior. In particular, a nonempty linear polyhedron need not be full-dimensional; a singleton or a hyperplane is an example.
Assume , increasing the coefficient bound to one if necessary. If and denotes the th row, then . For and , the Cauchy-Schwarz inequality gives
Thus with , proving that is a full-dimensional linear polyhedron. This part works for every positive ; the particular small value also ensures that relaxation cannot create feasibility when is empty.
Here is that exact-feasibility point. Put . If is empty, a Farkas certificate for linear inequalities gives with and . Normalize , and maximize over the compact linear polyhedron
Choose a maximizing extreme point. Its positive support has size : a dependence among its supported columns would permit small perturbations of both signs in , contradicting extremality. Select independent equations for those coordinates. Cramer's rule expresses them over a common nonzero integer determinant , with . Since is integer and , its numerator is a positive integer after choosing a positive denominator. Hence . This is the normalized integer Farkas infeasibility gap. A point in would imply
a contradiction. Consequently if and only if .
The ellipsoid method takes the integer matrix , vector , dimension , and a bound (or their binary encodings). It computes the prescribed , an inner radius , and an outer radius whose logarithms have polynomial size. For an explicit outer bound, let . A nonempty has a rational feasibility certificate for integer inequalities with . To see this bound, intersect with an orthant containing a feasible point. This intersection is a nonempty linear polyhedron containing no whole nontrivial affine line, and it has an extreme point. Indeed, if its active rows do not yet span the ambient space, move along a nonzero direction annihilating them. The orthant prevents motion in both directions from being feasible forever, so in at least one direction a new independent row becomes active. Repeating at most times produces an extreme point. At that point choose independent active constraints, including coordinate constraints where needed. Cramer's rule and the Hadamard determinant inequality bound the numerators by , while a nonzero integer denominator has absolute value at least one. Thus the point has norm at most . Choosing ensures that, whenever is nonempty, an entire radius- Euclidean ball in lies inside .
Start with the ellipsoid and retain . At each iteration test the centre against every relaxed inequality. If it passes all tests, is nonempty and hence so is . Otherwise a violated row gives the separation oracle: the feasible set lies in . Cut the current ellipsoid by this half-space through its centre and replace it by the standard containing ellipsoid. Every step retains and, for , reduces Lebesgue measure by a factor at most . In dimension one the corresponding step bisects an interval.
If more than cuts occur without finding a feasible centre, the Lebesgue measure is smaller than that of a radius- ball. This contradicts the retained ball whenever is nonempty, so the algorithm can declare empty. Both and are polynomial in and , and each inequality test uses arithmetic operations for rows. The bit implementation uses rational rounding with outward enlargement of the containing ellipsoids; the encoding bounds allow polynomial precision while keeping a constant fraction of the volume decrease. Detailed update and rounding formulae are not needed for this account. The essential role of the full-dimensional relaxation of integer inequalities is to supply a positive inner-volume bound without changing the emptiness decision; applying a Lebesgue measure stopping test directly to a lower-dimensional would fail.
The Euclidean algorithm gives
Back substitution produces the Bézout's identity
Hence a positive integer solution is . More generally, subtracting this identity from any other solution shows . Since and are coprime, the integer solutions are , ; both are positive precisely when .
Divide the linear congruence by the common factor , also dividing its modulus:
The Bézout's identity from part (a) makes the inverse of modulo . Multiplication by therefore gives . Thus all integer solutions are . Equivalently, the three solution classes modulo are ; using modulus must not discard two of these classes.
For the odd-square sum, the initial case gives . Assume the formula at . Adding the next square gives
The initial case and this induction step prove for every positive integer .
For the divisibility assertion, . The change when is replaced by is
One of the consecutive integers is even, so the displayed change is divisible by . If is divisible by , its successor is too. Thus induction proves for every .
Number the natural numbers from one. If the power set could be enumerated as , form . For each , membership of in is opposite to its membership in , so . This contradicts completeness of the enumeration. Thus the Cantor diagonal argument proves that the power set is an uncountable set.
For a finite subset , define . This is a nonnegative integer and the map is injective: if two finite subsets differ, let be their largest differing element. Its contribution exceeds the sum of all smaller possible contributions, so their codes differ. There are infinitely many singleton subsets. Hence the countability of finite subsets of a countable set here is countable infinitude.
To treat all bijections, associate to every subset the permutation that swaps and when , and fixes both when . These disjoint pair swaps define a bijection of . Different subsets give different permutations, so the map from the power set into is injective. Therefore is an uncountable set.
Every finitely supported permutation fixes all elements above some . A bijection fixing that tail must permute : it cannot send a smaller element into the fixed tail without violating injectivity. There are exactly such permutations. Thus
is a countable union of finite sets. Enumerate first by , then by one of the finitely many permutations at that ; duplicates can be omitted. Infinitely many distinct pair swaps belong to , so is countably infinite. The full family is uncountable, while the finitely supported family is countably infinite.
For a prime number and an integer coprime to , multiplication by permutes the nonzero residue classes modulo . Indeed, implies , and none of the products is zero. Taking the product over the classes gives
The factorial is coprime to and can be cancelled, proving Fermat little theorem: .
Now and is coprime to , so it is coprime to . Since , the proved theorem gives
Thus the two specified integer congruences modulo are equivalent. In this particular step the additional hypothesis is not needed; it will hold for the squarefree factorization used next.
If is a squarefree integer, its distinct prime factors are pairwise coprime. Consequently divides exactly when every prime factor does. Applying the preceding equivalence with proves
Finally retain this squarefree hypothesis. If divides , write . Fermat little theorem gives for every prime factor. Their pairwise coprimality then gives . For composite this is the sufficient direction of the Korselt criterion for a Carmichael number; primality is not excluded by the requested implication.