Print all safe primes below N
Last Updated :
23 Jun, 2022
Given an integer N, the task is to print all safe primes below N safe primes. A safe prime is a prime number of the form (2 * p) + 1 where p is also a prime.
The first few safe primes are 5, 7, 11, 23, 47, …
Examples:
Input: N = 13
Output: 5 7 11
5 = 2 * 2 + 1
7 = 2 * 3 + 1
11 = 2 * 5 + 1
Input: N = 6
Output: 5 7
Approach: First pre-compute all the primes till N using Sieve of Eratosthenes and then starting from 2 check whether the current prime is also a safe prime. If yes then print it else skip to the next prime.
Below is the implementation of above approach:
C++
#include <bits/stdc++.h>
using namespace std;
void printSafePrimes( int n)
{
int prime[n + 1];
for ( int i = 2; i <= n; i++)
prime[i] = 1;
prime[0] = prime[1] = 0;
for ( int p = 2; p * p <= n; p++) {
if (prime[p] == 1) {
for ( int i = p * 2; i <= n; i += p)
prime[i] = 0;
}
}
for ( int i = 2; i <= n; i++) {
if (prime[i] != 0) {
int temp = (2 * i) + 1;
if (temp <= n && prime[temp] != 0)
prime[temp] = 2;
}
}
for ( int i = 5; i <= n; i++)
if (prime[i] == 2)
cout << i << " " ;
}
int main()
{
int n = 20;
printSafePrimes(n);
return 0;
}
|
Java
class GFG{
static void printSafePrimes( int n)
{
int prime[] = new int [n + 1 ];
for ( int i = 2 ; i <= n; i++)
prime[i] = 1 ;
prime[ 0 ] = prime[ 1 ] = 0 ;
for ( int p = 2 ; p * p <= n; p++)
{
if (prime[p] == 1 )
{
for ( int i = p * 2 ; i <= n; i += p)
prime[i] = 0 ;
}
}
for ( int i = 2 ; i <= n; i++)
{
if (prime[i] != 0 )
{
int temp = ( 2 * i) + 1 ;
if (temp <= n && prime[temp] != 0 )
prime[temp] = 2 ;
}
}
for ( int i = 5 ; i <= n; i++)
if (prime[i] == 2 )
System.out.print(i + " " );
}
public static void main(String []args)
{
int n = 20 ;
printSafePrimes(n);
}
}
|
Python3
from math import sqrt
def printSafePrimes(n):
prime = [ 0 for i in range (n + 1 )]
for i in range ( 2 , n + 1 ):
prime[i] = 1
prime[ 0 ] = prime[ 1 ] = 0
for p in range ( 2 , int (sqrt(n)) + 1 , 1 ):
if (prime[p] = = 1 ):
for i in range (p * 2 , n + 1 , p):
prime[i] = 0
for i in range ( 2 , n + 1 , 1 ):
if (prime[i] ! = 0 ):
temp = ( 2 * i) + 1
if (temp < = n and prime[temp] ! = 0 ):
prime[temp] = 2
for i in range ( 5 , n + 1 ):
if (prime[i] = = 2 ):
print (i, end = " " )
if __name__ = = '__main__' :
n = 20
printSafePrimes(n)
|
C#
using System;
class GFG{
static void printSafePrimes( int n)
{
int [] prime = new int [n + 1];
for ( int i = 2; i <= n; i++)
prime[i] = 1;
prime[0] = prime[1] = 0;
for ( int p = 2; p * p <= n; p++)
{
if (prime[p] == 1)
{
for ( int i = p * 2; i <= n; i += p)
prime[i] = 0;
}
}
for ( int i = 2; i <= n; i++)
{
if (prime[i] != 0)
{
int temp = (2 * i) + 1;
if (temp <= n && prime[temp] != 0)
prime[temp] = 2;
}
}
for ( int i = 5; i <= n; i++)
if (prime[i] == 2)
Console.Write(i + " " );
}
public static void Main()
{
int n = 20;
printSafePrimes(n);
}
}
|
PHP
<?php
function printSafePrimes( $n )
{
$prime = array ();
for ( $i = 2; $i <= $n ; $i ++)
$prime [ $i ] = 1;
$prime [0] = $prime [1] = 0;
for ( $p = 2; $p * $p <= $n ; $p ++)
{
if ( $prime [ $p ] == 1)
{
for ( $i = $p * 2;
$i <= $n ; $i += $p )
$prime [ $i ] = 0;
}
}
for ( $i = 2; $i <= $n ; $i ++)
{
if ( $prime [ $i ] != 0)
{
$temp = (2 * $i ) + 1;
if ( $temp <= $n &&
$prime [ $temp ] != 0)
$prime [ $temp ] = 2;
}
}
for ( $i = 5; $i <= $n ; $i ++)
if ( $prime [ $i ] == 2)
echo $i , " " ;
}
$n = 20;
printSafePrimes( $n );
?>
|
Javascript
<script>
function printSafePrimes(n)
{
let prime = new Array(n + 1);
for (let i = 2; i <= n; i++)
prime[i] = 1;
prime[0] = prime[1] = 0;
for (let p = 2; p * p <= n; p++)
{
if (prime[p] == 1)
{
for (let i = p * 2; i <= n; i += p)
prime[i] = 0;
}
}
for (let i = 2; i <= n; i++)
{
if (prime[i] != 0)
{
let temp = (2 * i) + 1;
if (temp <= n && prime[temp] != 0)
prime[temp] = 2;
}
}
for (let i = 5; i <= n; i++)
if (prime[i] == 2)
document.write(i + " " );
}
let n = 20;
printSafePrimes(n);
</script>
|
Time complexity: O(nlog(logn)) same as Sieve of Eratosthenes
Auxiliary Space: O(n), since n extra space has been taken.
Share your thoughts in the comments
Please Login to comment...