使用Divide&amp ;;计算矩阵中的元素。征服

时间:2015-06-19 06:33:26

标签: c++ algorithm matrix divide-and-conquer

我开始学习如何实现分治算法,但是我在练习时遇到了一些麻烦。

我编写了一个算法但很遗憾它返回0值。我需要计算一个数字使用D& C在矩阵中找到的次数。

我的想法是将矩阵划分为4个子矩阵,直到所有坐标都相等,这意味着有一个元素(可能是我的数字)。

我的缩写:

CloseableHttpClient

当我编译时,它运行良好但返回0 (因为我认为我的停止语句)。 所以,我的错误在于solve()函数。

我的代码:

lr - left row   | starting with 
lc- left column | left lower corner
rr - right row    | starting with
rc- right column  |  right upper corner

也许你可以告诉我问题出在哪里或直接解决问题,它会更快,我可以更容易理解。

#include<fstream>
#include<conio.h>
#include<iostream>
using namespace std;
ifstream fin("in.txt");
#define debug cerr<<"OK";
int n,a[100][100];
int solve(int lr,int lc, int rr,int rc,int x,int &contor)
{
    int cnt=0;
    if(lr < rr || lc > rc)
    {
        if(cnt<=0)
            return 0;
        return cnt;
    }
        if(lr == lc && rr == rc)
        if(a[lr][lc] == x)
            cnt++;
        else;
    else
    {
        int l=(lr+rr)/2;
        int c=(lc+rc)/2;
        if(l && c)
        {
        int cnt1=solve(lr-l,lc,rr,rc-c,x,contor); // coordinates for upper left square
        int cnt2=solve(lr-l,lc+c,rr,rc,x,contor); // -||- for upper right square
        int cnt3=solve(lr,lc,rr+l,rc-c,x,contor); // -||- for low left square 
        int cnt4=solve(lr,lc+c,rr,rc-c,x,contor); // -||- for low right square
        contor=cnt1+cnt2+cnt3+cnt4;
        }
    }
}

int main()
{  
    fin>>n;
    int i,j;
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++)
        fin>>a[i][j];
    int contor=0;
solve(n,1,1,n,3,contor);  // i'm searching for 3
   cout<<contor;
    _getch();
   fin.close();
   return 0;
}

2 个答案:

答案 0 :(得分:0)

递归不会更新cnt ...

    if(l && c)
    {
    solve(lr-l,lc,rr,rc-c,x); // coordinates for upper left square 
    solve(lr-l,lc+c,rr,rc,x); // -||- for upper right square
    solve(lr,lc,rr+l,rc-c,x); // -||- for low left square 
    solve(lr,lc+c,rr,rc-c,x); // -||- for low right square
    }

答案 1 :(得分:0)

以下是更正后的源代码。我添加了一些评论,以便您可以关注它。此外,我切换到了正确大小的动态数组。

#include<fstream>
#include<conio.h>
#include<iostream>

using namespace std;

int countOccurences(int** values, int min1, int min2, int max1, int max2, int searchValue)
{
    if (min1 > max1 || min2 > max2)
        return 0; //invalid area

    if (min1 == max1 && min2 == max2)
    {
        //the current area is a single cell
        if (values[min1][min2] == searchValue)
            return 1; //we have found 1 occurence
        else
            return 0; //we have found nothing
    }
    else
    {
        //divide the current range
        int center1 = (min1 + max1) / 2;
        int center2 = (min2 + max2) / 2;

        int occurencesInSubMatrices = 0;

        //accumulate the occurences in the according variable
        occurencesInSubMatrices += countOccurences(values, min1, min2, center1, center2, searchValue); // coordinates for upper left square
        occurencesInSubMatrices += countOccurences(values, center1 + 1, min2, max1, center2, searchValue); // -||- for upper right square
        occurencesInSubMatrices += countOccurences(values, min1, center2 + 1, center1, max2, searchValue); // -||- for low left square 
        occurencesInSubMatrices += countOccurences(values, center1 + 1, center2 + 1, max1, max2, searchValue); // -||- for low right square
        return occurencesInSubMatrices;
    }
}

int main()
{
    ifstream fin("in.txt");
    int n;
    fin >> n;

    //create a 2d array of appropriate size
    int** values = new int*[n];

    for (int i = 0; i < n; i++)
    {
        //allocate memory for the i-th row
        values[i] = new int[n];
        for (int j = 0; j < n; j++)
            fin >> values[i][j];
    }

    int count = countOccurences(values, 0, 0, n - 1, n - 1, 3);  // i'm searching for 3
    cout << count;
    _getch();
    fin.close();
    return 0;
}

请注意,这是纯粹的教育范例。在现实世界中,没有人会像这样计算这种情况。