인쇄 낮은 삼각 포장 매트릭스 나는 포장 낮은 삼각 행렬로 배열에 저장된 실제 * 8 개 숫자의 행렬이 77

Question

: 나는 다음과 같은 형식으로 인쇄 할 수 있도록하려면

|1 * * * * *| 
|2 7 * * * *| 
|3 8 12 * * *| 
|4 9 13 16 * *| => [1,2,3,4,5,6,7,8,9...21] 
|5 10 14 17 19 *| 
|6 11 15 18 20 21| 

을 :

가

[row|col] 1 2 3 4 5 
    1  1 * * * * 
    2  2 7 * * * 
    3  3 8 12 * * 
    4  4 9 13 16 * 
    5  5 10 14 17 19 
    6  6 11 15 18 20 

[row|col] 6 
    1  * 
    2  * 
    3  * 
    4  * 
    5  * 
    6  6 

내가 가진 문제는 내가 모르는 얼마나 열별로 요소를 저장하고 인쇄하려면 별도의 배열을 생성 할 필요없이 행렬 요소로 루프. 이것은 지금까지 시도한 것입니다.

 Implicit Real*8 (A-H,O-Z) 


     INTEGER  ARRAY(21) 



     10 Format(5X,'[Row|Col]',5(8X,I6)) 
     Icol=6 
     Num1=1 
C  Test array 
     DO K=1,21 
     Array(K)=K*1.0d0 
     ENDO 

C  Print the elements row-by-row.  
44 Num2=Icol-Num1 
     Num=Num2 
     If ((Num2).gt.5) Num=5 
     Write(*,10) (I,I=Num1,Num1+Num) 
     INum1=INum1+Num 
     if (Inum1.ne.Icol) goto 44  
     STOP 
     END 

Answer 1

키는 낮은 삼각형 행렬의 0이 아닌 요소를 지정하는 수식을 식별하는 것입니다. i와 j가 각각 행과 열에 대한 인덱스를 실행하는 경우 요소 (j < i)는 아래 삼각형 부분에 있습니다. 여기에 코드입니다 :

PROGRAM print_low_tri_matrix 
IMPLICIT NONE 
INTEGER :: i,j,n,m,p 
REAL ,ALLOCATABLE, DIMENSION(:,:) :: a 
REAL, ALLOCATABLE, DIMENSION(:) :: r 

n = 6 
p = n*(n+1)/2 ! Number of non-zero elements in a lower triangular matrix of size n by n 

ALLOCATE(a(n,n),r(p)) 

DO i = 1, p 
     r(i) = i*1.0 ! Array containing the non-zero elements 
END DO 

m = 1   ! Index tracking the non-zero elements array r 
DO j = 1,n 
     DO i = 1,n 
       IF(j .LE. i) THEN  ! Non-zero indices of the matrix 
         a(i,j) = r(m) 
         m = m + 1 
       ELSE 
         a(i,j) = 0. 
       END IF 
     END DO 
END DO 


!  Printing the full matrix a - uncomment for checking 
!   DO i = 1,n 
!       WRITE(*,*) (a(i,j),j=1,n) 
!   END DO 


!  Printing the matrix a as needed 
DO i = 1,n 
     DO j = 1,n 
       IF(j .LE. i) THEN 
         WRITE(*,'(F12.4)',ADVANCE="NO") a(i,j) 
       ELSE 
         WRITE(*,'(A12)',ADVANCE="NO") '*' 
       END IF 

       IF (j .EQ. n) WRITE(*,*) 
     END DO 
END DO 

END PROGRAM print_low_tri_matrix 

그리고 출력은 다음과 같습니다 당신이 원하는 경우 물론

1.0000   *   *   *   *   * 
2.0000  7.0000   *   *   *   * 
3.0000  8.0000  12.0000   *   *   * 
4.0000  9.0000  13.0000  16.0000   *   * 
5.0000  10.0000  14.0000  17.0000  19.0000   * 
6.0000  11.0000  15.0000  18.0000  20.0000  21.0000 

, 당신은 정수 배열로 행렬을 변경할 수 있습니다.

P.S : - f77에서이 작업을 수행하려면 WRITE(*,'(F12.4)',ADVANCE="NO") a(i,j)에서 WRITE(*,'(F12.4,$)') a(i,j)으로 형식 지정자를 변경할 수 있습니다. 배열의 할당은 직접 제거하고 정의해야합니다.