Sudoku en C ou C++

Pour chaque case (i,j) de la grille
  Si le domaine de la case(i,j) est réduit à un élément : le placer
  Sinon : tirer une valeur v au hasard parmi les valeurs encore dans le domaine
  Pour chaque case de la ligne i : supprimer la valeur v du domaine
  Pour chaque case de la colonne j : supprimer la valeur v du domaine
  Pour chaque case du carré n x n contenant (i,j) : supprimer la valeur v du domaine
Fin pour

#define n 3 // sudoku 3x3 cadrans de 3x3 cases

// Associer chaque case (ligne,colonne) à un index qui identifie un cadran
// (le cadran en haut à geuche a pour 0, celui à sa droite 1 etc...) ce qui donne
// ce mapping pour n = 3 :
// 000111222
// 000111222
// 000111222
// 333444555
// 333444555
// 333444555
// 666777888
// 666777888
// 666777888

// Retourne l'identifiant du cadran contenant la case (ligne, colonne)
unsigned coordonnees2cadran(unsigned ligne,unsigned colonne){
  if (ligne >= n*n) throw;
  if (colonne >= n*n) throw;
  return (ligne/n)*n + colonne/n; // '/' : division euclidienne
}

Vider les domaines associés à chaque ligne
Vider les domaines associés à chaque colonne
Vider les domaines associés à chaque cadran
Pour chaque ligne i
  Pour chaque colonne j
    Ajouter la valeur de la case(i,j) dans le domaine de la ligne(i)
    Ajouter la valeur de la case(i,j) dans le domaine de la colonne(j)
    Ajouter la valeur de la case(i,j) dans le domaine de du cadran contenant (i,j)
  Fin pour
Fin pour

Pour chaque ligne i
  Si le domaine de la ligne(i) ne contient pas les 9 valeurs : retourner faux
Fin pour

Pour chaque colonne j
  Si le domaine de la colonne(j) ne contient pas les 9 valeurs : retourner faux
Fin pour

Pour chaque cadran
  Si le domaine du cadran ne contient pas les 9 valeurs : retourner faux
Fin pour

Retourner vrai

#include <vector>

typedef std::vector<bool> domaine_t;

struct case_t{
  domaine_t domaine; // pour générer la grille ou pour le solveur
  unsigned valeur;
};

typedef std::vector<std::vector<case_t> > grille_t;

#include <iostream>
#include <vector>
#include <set>
#include <cstdlib>
#include <cassert>

#define N     3
#define N2    N*N

typedef std::set<unsigned> domaine_t;

unsigned randint(unsigned max){
 return (unsigned) (((double) rand()*(max+1))/RAND_MAX);
}

struct case_t{
	domaine_t domaine;
	unsigned  valeur;

	case_t():
		valeur(0)
	{
		for(unsigned i=1;i<N2+1;++i) domaine.insert(i);
	}

	void supprimer(unsigned val){
		domaine_t::iterator fit = domaine.find(val);
		if(fit != domaine.end()) domaine.erase(fit);
	}
};

struct grille_t{
	std::vector<std::vector<case_t> > data;

	grille_t(){
		for(unsigned i=0;i<N2;++i){
			data.push_back(std::vector<case_t>(N2,case_t()));
		}
	}

	void corriger(unsigned i,unsigned j,unsigned val){
		assert(i < N2);
		assert(j < N2);

		// Virer val de la ligne i et de la colonne j
		for(unsigned k=0;k<N2;++k){
			data[i][k].supprimer(val);
			data[k][j].supprimer(val);
		}

		// Virer val du cadran dont la case supérieure gauche
		// est (i0,j0)
		unsigned i0 = i / N2;
		unsigned j0 = j / N2;
		for(unsigned i1=0; i1<N; ++i1){
			for(unsigned j1=0; j1<N; ++j1){
				data[i0+i1][j0+j1].supprimer(val);
			}
		}
	}

	void generer(){
		unsigned val, k, rnd;
		for(unsigned i=0;i<N2;++i){
			for(unsigned j=0;j<N2;++j){
				const domaine_t & d = data[i][j].domaine;
				rnd = randint(d.size()-1);
				
				// Tirer la valeur de la case (i,j)
				// (la rnd-ième valeur du domaine)
				domaine_t::const_iterator
					sit (d.begin()),
					send(d.end());
				for(k=0;sit!=send && k < rnd;++sit, ++k);
				val = *sit;
				data[i][j].valeur = val;
				
				// Corriger les domaines
				corriger(i,j,val);
			}
		}
	}
};

std::ostream & operator << (std::ostream & out,const grille_t & g){
	for(unsigned i = 0; i < N2; ++i){
		if(i%N == 0){
			for(unsigned j = 0; j < N2 + N; ++j) out << '-';
			out << std::endl;
		}
		for(unsigned j = 0; j < N2; ++j){
			if (j%N == 0) out << '|';
			out << g.data[i][j].valeur;
		}
		out << std::endl;
	}
	return out;
}

int main(){
//	for(unsigned k=0;k<20;++k) std::cout << randint(9) << std::endl;
	grille_t grille;
	grille.generer();
	std::cout << grille << std::endl;
	return 0;
}

------------
|847|691|352
|536|478|910
|192|350|476
------------
|283|765|194
|674|932|580
|951|840|237
------------
|725|183|649
|419|207|803
|368|509|721