c++ acm题目1057 DNA培养问题

以下是题目和我的代码，提交之后是wrong answer，还请各位大神帮忙看看哪里错误，万分感谢！！！

Description
A biologist experimenting with DNA modification of bacteria has found a way to make bacterial colonies sensitive to the surrounding population density. By changing the DNA, he is able to “program” the bacteria to respond to the varying densities in their immediate neighborhood.
The culture dish is a square, divided into 400 smaller squares (20x20). Population in each small square is measured on a four point scale (from 0 to 3). The DNA information is represented as an array D, indexed from 0 to 15, of integer values and is interpreted as follows:
In any given culture dish square, let K be the sum of that square’s density and the densities of the four squares immediately to the left, right, above and below that square (squares outside the dish are considered to have density 0). Then, by the next day, that dish square’s density will change by D[K] (which may be a positive, negative, or zero value). The total density cannot, however, exceed 3 nor drop below 0.
Now, clearly, some DNA programs cause all the bacteria to die off (e.g., [-3, -3, …, -3]). Others result in immediate population explosions (e.g., [3,3,3, …, 3]), and others are just plain boring (e.g., [0, 0, … 0]). The biologist is interested in how some of the less obvious DNA programs might behave.
Write a program to simulate the culture growth, reading in the number of days to be simulated, the DNA rules, and the initial population densities of the dish.

Input
Input to this program consists of three parts:

1. The first line will contain a single integer denoting the number of days to be simulated.
2. The second line will contain the DNA rule D as 16 integer values, ordered from D[0] to D[15], separated from one another by one or more blanks. Each integer will be in the range -3…3, inclusive.
3. The remaining twenty lines of input will describe the initial population density in the culture dish. Each line describes one row of squares in the culture dish, and will contain 20 integers in the range 0…3, separated from one another by 1 or more blanks.

Output
The program will produce exactly 20 lines of output, describing the population densities in the culture dish at the end of the simulation. Each line represents a row of squares in the culture dish, and will consist of 20 characters, plus the usual end-of-line terminator.
Each character will represent the population density at a single dish square, as follows:
No other characters may appear in the output.

Sample Input
1
2
0 1 1 1 2 1 0 -1 -1 -1 -2 -2 -3 -3 -3 -3
3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

Sample Output
##!..
#!..
!..
…
…
…
…
…!..
…!#!..
…!#X#!..
…!#!..
…!..
…
…
…
…
…
…
…
…

#include
#include<stdio.h>
using namespace std;
int main(){
int day,d=1;//day是模拟的天数；
int dna[16];//DNA信息
int density1[22][22]={0},density2[22][22]={0};
int i,j;
cin>>day;
for(i=0;i<16;i++)scanf("%d",&dna[i]);
for(i=1;i<21;i++)
for(j=1;j<21;j++)
scanf("%d",&density1[i][j]);
while(day–){
d=1;//d代表现在是第几天
for(i=1;d%2!=0&&i<21;i++)
for(j=1;j<21;j++)
density2[i][j]=density1[i][j]+density1[i+1][j]+density1[i-1][j]+density1[i][j+1]+density1[i][j-1];
for(i=1;d%2==0&&i<21;i++)
for(j=1;j<21;j++)
density1[i][j]=density2[i][j]+density2[i+1][j]+density2[i-1][j]+density2[i][j+1]+density2[i][j-1];
d++;
}
for(i=1;d%2!=0&&i<21;i++){
for(j=1;j<21;j++){
if(density2[i][j]==0){cout<<".";continue;}
else if(density2[i][j]==1){cout<<"!";continue;}
else if(density2[i][j]2){cout<<“X”;continue;}
else cout<<"#";
}
cout<<endl;
}
for(i=1;d%20&&i<21;i++){

	for(j=1;j<21;j++){
	    if(density1[i][j]==0){cout<<".";continue;}
        else if(density1[i][j]==1){cout<<"!";continue;}
        else if(density1[i][j]==2){cout<<"X";continue;}
        else cout<<"#";
	}
	cout<<endl;
}
return 0;

}