C语言求解布尔矩阵的布尔积运算

让我们先回顾下布尔矩阵的布尔积运算的公式

如果A=[a_ij]是m×p矩阵,B=[b_ij]是p×n矩阵,则A和B的布尔积(Boolean Product)是m×n矩阵C=A·B=[c_ij],其中

       c_ij=(a_i1∧b_1j) ∨(a_i2∧b_2j) ∨……∨(a_ip∧b_pj)  

#include<stdio.h>
//输入n×m的矩阵
void input(int a[][200],int n,int m);
//矩阵A(大小m×l) 和矩阵B(大小l×n)的布尔积运算
void solve(int a[][200],int b[][200],int m,int l,int n);
//输出n×m的矩阵
void output(int a[][200],int n,int m);
int main()
{
    int m,l,n;
    int a[200][200],b[200][200];

    while(scanf("%d %d",&m,&l)!=EOF)
    {
        input(a,m,l);
        scanf("%d %d",&l,&n);
        input(b,l,n);
        solve(a,b,m,l,n);
    }
}
void output(int a[][200],int n,int m)
{
    int i,j;
    for(i=0;i<n;i++)
    {
        for(j=0;j<m;j++)
        {
            if(j!=0)
                printf(" ");
            printf("%d",a[i][j]);
        }
        printf("\n");
    }
}
void input(int a[][200],int n,int m)
{
    int i,j;
    for(i=0;i<n;i++)
        for(j=0;j<m;j++)
            scanf("%d",&a[i][j]);

}
void solve(int a[][200],int b[][200],int m,int l,int n)
{
    int c[200][200];
    int i,j,k;
    int flag;
    for(i=0;i<m;i++)
    {
        for(j=0;j<n;j++)
        {
            flag=0;
            for(k=0;k<l;k++)
            {
                if(a[i][k]==1&&b[k][j]==1)
                {
                    flag=1;
                    break;
                }
            }
            c[i][j]=flag;
        }
    }
    output(c,m,n);
}

测试数据如下

3 4
1 0 0 1
1 1 0 1
1 1 0 1
4 4
0 0 0 0
0 1 1 1
1 1 0 1
0 1 0 1

结果如下

0 1 0 1
0 1 1 1
0 1 1 1