我正在尝试用FORTRAN创建帕斯卡三角形。我做了算法。用 C 编译并成功,但由于某种原因,我在 FORTRAN 中没有得到相同的预期结果。谁能帮我解决这个问题?
C语言代码(工作):
#include <stdio.h>
#include <stdlib.h>
int main()
{
unsigned int c, i, j, k,n;
scanf("%d",&n);
for(i=0; i < n; i++) {
c = 1;
for(j=1; j <= (n-1-i); j++) printf(" ");
for(k=0; k <= i; k++) {
printf("%2d", c);
c = c * (i-k)/(k+1);
}
printf("n");
}
return 0;
}
FORTRAN中的代码(不起作用,在这里需要帮助):
program pascal
implicit none
integer i,j,k,p,n
read(*,*)n
i=0
do while(i.lt.n)
p=1
do j=1,n-1-i
write(*,5)
5 format(1x)
enddo
do k = 0,i
write(*,1)p
1 format(i2)
p = p*(i-k)/(k+1)
enddo
i=i+1
write(*,2)
2 format(/)
enddo
endprogram
一些可以处理较大n
值的改进:
program pascal
implicit none
integer i,j,k,p,n
write (*, '("input n: ")', advance="no")
read(*,*) n
do i=0,n-1
p=1
do j=1,n-1-i
write(*,'(3X)', advance="no")
enddo
do k = 0,i
write(*,'(I6)', advance="no") p
p = p*(i-k)/(k+1)
enddo
write(*, '(/)')
enddo
endprogram
该程序的格式化IO对我来说似乎很简单。 如果您想找出将输出的特定值的最大整数n
并输出该情况下最紧凑的三角形,则更复杂......
试试这个。请注意,如果您输入数字 n> 9,它将不起作用(对于 n> 5,它不是那么漂亮)。
program pascal
implicit none
integer :: i, j, n, c
read(*, *) n
do i = 0, n - 1
c = 1
do j = 0, n - i
write(*, "(' ')", advance="no")
end do
do j = 0, i
write(*, "(I2)", advance="no") c
c = (c * (i - j)) / (j + 1)
end do
write(*, *)
end do
end program
在您的程序中,几乎所有内容都是正确的,除了当您打印时:如果您不添加可选的advance="no"
,Fortran 总是添加一个新行。
关于样式,在这里将格式放在单独的行中不是很有用,.lt.
可以在Fortran 90中替换为<
。此外,内部循环中的缩进不正确。