1

我正在尝试用 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
4

2 回答 2

1

以下是一些可以处理较大值的改进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并为该情况输出最紧凑的三角形，则更复杂...

于 2014-04-04T02:19:01.650 回答
0

你可以试试这个。请注意，如果您输入一个数字 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

在您的程序中，几乎所有内容都是正确的，除非您 print ：如果您不添加可选的 .Fortran，则始终添加新行advance="no"。

关于样式，在这里将格式放在单独的行中并不是很有用，并且.lt.可以<在 Fortran 90 中替换为。此外，您的缩进在内部循环中不正确。

于 2014-04-03T18:00:09.950 回答