RECHERCHE DES NOMBRES CHANCEUX SELON STANISLAW ULAM

#include<iostream>
#include<iomanip>
using namespace std;

//
//Fonction qui renvoie l'indice du prochain élément 
//du tableau qui vaut 'true', càd le prochain entier
//encore actuellement "chanceux" ou 0 s'il n'y a pas
//de prochain nombre inférieur à la limite.
///////////////////////////////////////////////////
int prochainVrai( bool *tab, int n, int taille )
{
	for ( int i=n+1; i<taille; i++)
		if ( tab[i] )
			return i;
	return 0;
}

//
//Recherche et affichage des nombres 
//chanceux inférieurs à "limite"
////////////////////////////////////
void chanceux( int limite )
{
	bool tab[limite+1];
	int pas=1, tmp, cpt;
	
	//traitement de la première étape, on élimine les nombres impairs
	for (int i =0; i<=limite; i++)
		tab[i]=( i%2 != 0 );
	
	//chaque itération de cette boucle représente une étape,
	//on itère donc jusqu'à ce qu'il n'y ait plus
	//de nombres inférieurs à la limite
	while ( (pas=prochainVrai(tab,pas,limite+1)) != 0 )
	{		
		tmp=0;
		cpt=0;
		
		//suppression des nombres correspondant à l'étape en cours
		while ( (tmp=prochainVrai(tab,tmp,limite+1)) != 0)
		{	
			cpt++;	
			
			if (cpt == pas)
			{	
				tab[tmp]=false;
				cpt=0;
			} 
		}
	}
	
	// affichage des nombres chanceux	
	for(int i=1; i<=limite; i++)
		if ( tab[i] )
			cout << setw(7) << i << endl;
}


int main()
{
	int lim= -1;
	
	while ( lim < 0 )
	{	
		cout << "Entrez la limite de recherche des nombres chanceux : ";	
		cin >> lim;
	}
	
	chanceux(lim);
	
	return 0;
}