Найти несколько максимальных чисел и запомнить их индекс

#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <clocale>
#include <string.h>
#include <istream>
#include <ctime>
 
using namespace std;
 
int _tmain(int argc, _TCHAR* argv[])
{
    setlocale(LC_ALL, "rus");
    double **matr, *str;
    int N, M, CHIS;
    double mas[5];
    memset(mas, '0', 5);
    int j = 0, i = 0, b = 0;
    int count = 0, n = 0, max = 0; int index = 0;
    cout << "Введите кол-во строк = ";
    cin >> N;
    cout << "Введите кол-во столбцов = ";
    cin >> M;
    cout << "Введите кол-во наибольшее число матрицы (10 - веществ. число до 1, 100 - до 10 и 1000 - до 100) = ";
    cin >> CHIS;
    str = new double[M];
    matr = new double*[N];
    for (i = 0; i < N; i++)
        matr[i] = new double[M];
    srand(time(NULL));
    for (i = 0; i < N; i++)
    {
        for (j = 0; j < M; j++)
        {
            matr[i][j] = static_cast<double>(rand() % CHIS) / 10;//создание вещественной квадратной матрицы
        }
    }
    for (i = 0; i < N; i++)
    {
        cout << endl;
        for (j = 0; j < M; j++)
            cout << setw(5) << matr[i][j];
        cout << endl;
    }
    for (j = 0; j < M; j++)
    {
        for (i = 0; i < N; i++)
        {
            if (matr[i][j] > N && M)//цикл определения столбцов с числами больше заданного размера
            {
                count++;
                cout << "Столбец " << j << endl;
            }
        }
        str[j] = count;//запоминание в новом массиве кол-ва элементов в каждом столбец 
        count = 0;
    }
    for (j = 0; j < M; j++)
    {
        for (i = 0; i < 5; i++)
        {
            if (str[j] >= max)//и тут у меня находит и запоминает лишь один элемент
            {
                max = str[j];
                mas[i] = j;
            }
        }
    }
    for (i= 0; i < 5; i++)
        cout << mas[i];
    cout << endl;
    for (j = 0; j < M; j++)
        cout << str[j];
    cout << endl;
    cout << max<<endl;
    cout << "Какой вывести столб?";
    cin >> n;
        for (i = 0; i < N; i++)
            cout << matr[i][n]<<" ";
    system("pause");
    return 0;
}