1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90
| #include <iostream>
using namespace std ;
int main ()
{
const int n=10 ;
// Map, M(x,y)
float map[n+2][n+2] ;
for (int i=0;i<n;i++)
{
for (int j=0;j<n;j++)
{
map[i][j]=0 ;
}
}
// Densité
float d ;
// Densité max
float dmax=10 ;
for (int cadre=1;cadre<n;cadre=cadre+1)
{
for (int x=0;x<10;x=x+1)
{
if (map[x][cadre]<dmax)
{
if (x=0)
{
d=(((map[x+1][cadre])/2)+map[x][cadre+1])/2 ;
if (map[x][cadre]+2*d<=dmax)
{
map[x][cadre]=map[x][cadre]+2*d ;
map[x+1][cadre]=map[x+1][cadre]-d ;
map[x][cadre+1]=map[x][cadre+1]-d ;
}
else
{
map[x][cadre]=dmax ;
d=d-dmax/2 ;
map[x+1][cadre]=map[x+1][cadre]-d ;
map[x][cadre+1]=map[x][cadre+1]-d ;
}
}
else
{
d=(((map[x+1][cadre])/2)+(map[x][cadre+1])/2)/2 ;
if (map[x][cadre]+2*d<=dmax)
{
map[x][cadre]=map[x][cadre]+2*d ;
map[x+1][cadre]=map[x+1][cadre]-d ;
map[x][cadre+1]=map[x][cadre+1]-d ;
}
else
{
map[x][cadre]=dmax ;
d=d-dmax/2 ;
map[x+1][cadre]=map[x+1][cadre]-d ;
map[x][cadre+1]=map[x][cadre+1]-d ;
}
}
}
}
}
return 0 ;
} |