include <stdio.h>

    int main()
    {

    int matrixA[3][3];
    int matrixB[3][3];
    int result[3][3];
    int a = 0;
    int b = 0;

    printf("\n Boolean Power of Matrix\n");

    printf("\n Matrix of order 3x3\n\n");



    for(b=0;b<3;b++)
    {
        for(a=0;a<3;a++)
        {
        printf("Insert  Column for matrix  A %d Row %d : ",a+1,b+1);
        scanf("%d",&matrixA[a][b]);
        matrixB[a][b]=matrixA[a][b];

        }

    }

    printf("\n");


    for(b=0;b<3;b++)
    {
        for(a=0;a<3;a++)
        {
        printf("  %d",matrixB[a][b]);

        }
        printf("\n");

    }

 // for matrix A


    for(b=0;b<3;b++)
    {
        for(a=0;a<3;a++)
        {
        printf("Insert  Column  for matrix B %d Row %d :\n ",a+1,b+1);
        scanf("%d",&matrixB[a][b]);
        matrixB[a][b]=matrixB[a][b];

        }

    }

    printf("\n");


    for(b=0;b<3;b++)
    {
        for(a=0;a<3;a++)
        {
        printf("  %d",matrixB[a][b]);

        }
        printf("\n");

    }

  // for matrix B

    printf("Boolean product Matrix of order 3x3\n");


    result[0][0]=(matrixB[0][0]*matrixA[0][0])+(matrixB[0][1]*matrixA[1][0])+(matrixB[0][2]*matrixA[2][0]);
    if (result[0][0]==2 || result[0][0]==3)
    {
        result[0][0]=1;
    }
    result[0][1]=(matrixB[0][0]*matrixA[0][1])+(matrixB[0][1]*matrixA[1][1])+(matrixB[0][2]*matrixA[2][1]);
    if(result[0][1]==2 || result[0][1]==3)
    {
        result[0][1]=1;
    }
    result[0][2]=(matrixB[0][0]*matrixA[0][2])+(matrixB[0][1]*matrixA[1][2])+(matrixB[0][2]*matrixA[2][2]);
    if(result[0][2]==2 || result[0][2]==3)
    {
        result[0][2]=1;
    }
    result[1][0]=(matrixB[1][0]*matrixA[0][0])+(matrixB[1][1]*matrixA[1][0])+(matrixB[1][2]*matrixA[2][0]);
    if(result[1][0]==2 || result[1][0]==3)
    {
        result[1][0]=1;
    }
    result[1][1]=(matrixB[1][0]*matrixA[0][1])+(matrixB[1][1]*matrixA[1][1])+(matrixB[1][2]*matrixA[2][1]);
    if(result[1][1]==2 || result[1][1]==3)
    {
        result[1][1]=1;
    }
    result[1][2]=(matrixB[1][0]*matrixA[0][2])+(matrixB[1][1]*matrixA[1][2])+(matrixB[1][2]*matrixA[2][2]);
    if(result[1][2]==2 || result[1][2]==3)
    {
        result[1][2]=1;
    }

    result[2][0]=(matrixB[2][0]*matrixA[0][0])+(matrixB[2][1]*matrixA[1][0])+(matrixB[2][2]*matrixA[2][0]);
    if(result[2][0]==2 || result[2][0]==3)
    {
        result[2][0]=1;
    }
    result[2][1]=(matrixB[2][0]*matrixA[0][1])+(matrixB[2][1]*matrixA[1][1])+(matrixB[2][2]*matrixA[2][1]);
    if(result[2][1]==2 || result[2][1]==3)
    {
        result[2][1]=1;
    }
    result[2][2]=(matrixB[2][0]*matrixA[0][2])+(matrixB[2][1]*matrixA[1][2])+(matrixB[2][2]*matrixA[2][2]);
    if(result[2][2]==2 || result[2][2]==3)
    {
        result[2][2]=1;
    }


    for(b=0;b<3;b++)
    {
        for(a=0;a<3;a++)
        {
        printf(" %d ",result[a][b]);

        }
        printf("\n");

    }
    return 0;
    }