Поиск наиболее часто встречающейся подстроки в строке

#include <iostream>
#include <string>
using namespace std;
 
int main()
{
    string str, podstr, maxpodstr;
    // строка, текущая подстрока, подстр., которая встречается максимальное количество раз
    int maxk = 0, k, pos; // количество maxpodstr, количество текущей podstr, последняя найденная позиция podstr
    cin >> str;
    for (int i = 0; i < str.length(); i++) // с первого по последний символ
        for (int j = 1; j < str.length()-i; j++) // берем с одного до конца строки символов
        { // и смотрим сколько раз эта подстрока встретится
            podstr = str.substr(i, j);
            k = (str.find(podstr, 0) == 0? 1 : 0), pos = 0; // если на нулевой позиции есть эта подстрока,
                                                            // то одно вхождение уже есть
            do
            {
                pos = str.find(podstr, pos+1); // находим следующее вхождение (при pos==0 ищем со второго символа)
                if (pos != string::npos) // если нашли, увеличиваем количество вхождений, иначе выходим из цикла
                    k++;
                else
                    break;
            }
            while (true);
            if (k > maxk) // если нашли чаще встречающуюся подстроку
            {
                maxk = k; // запоминаем ее
                maxpodstr = podstr;
            }
        };
    cout << maxpodstr;
    return 0;
}

/////////////////////////////////////////////////////////////////////////////////////////
//дано натуральное число k <= 10^9. найдите все такие упорядоченные пары (a, b), 
//(a > 0, b > 0, b <= 10^18, a <= 10^18), такие что сумма чисел в этой паре ровно в к раз 
//меньше чем произведение чисел в паре.
//
//Эта задачка была на соревновании, надеюсь никто не будет ругаться(надеюсь, тут никого нет, 
//кто писал это соревнование xD).
//
//Формат ввода: одно число k(k > 0 && k <= 10^9);
//Формат вывода: выведите все пары, удовлетворяющие условию.
//
//Примеры:
//input: 1
//output: 2 2
//
//Еще один
//input: 2
//output:3 6
//6 3
//4 4
/////////////////////////////////////////////////////////////////////////////////////////
#include <iostream>
/////////////////////////////////////////////////////////////////////////////////////////
typedef long long   T_int;
/////////////////////////////////////////////////////////////////////////////////////////
void    print_pairs( T_int    k )
{
    for( T_int  p = 1; p <= k; ++p )
    {
        if  (
                k * k % p   == 0
            )
        {
            std::cout   <<  k + p
                        <<  '\t'
                        <<  k + k * k / p
                        <<  std::endl;
        }
    }//for
}
/////////////////////////////////////////////////////////////////////////////////////////
int     main()
{
    for(;;)
    {
        std::cout   <<  std::endl
                    <<  std::endl
                    <<  std::endl
                    <<  std::endl
                    <<  "k = ";
        T_int   k   =   0;
        std::cin    >>  k;
        print_pairs(k);
    }//for
}

typedef unsigned long long int uintll;
 
int _tmain(int argc, _TCHAR* argv[]) 
{
    uintll k, f[100], p[500][2]; // размеры массивов от балды, можно увеличить если надо
    cout << "Input k: "; cin >> k;
 
    uintll k_=k; int i_f=0, i=1;
    while (k_>1 && i++) if(k_%i==0) {f[i_f++]=i; k_/=i; i=1;} // факторизация
    int nf=i_f;
 
    p[0][0]=1; p[0][1]=1; int np=1;
    uintll fx = 1;
    for (int i_f=0; i_f<nf; i_f++) {
        fx = (i_f==0 || f[i_f]!=f[i_f-1]) ? 1 : fx*=f[i_f];
        int iip = np;
        for (int i_p=0; i_p<np; i_p++)
            if (fx==1 || p[i_p][1]%fx==0) {p[iip][0]=p[i_p][0]; p[iip][1]=p[i_p][1]*f[i_f]; iip++;}
        for (int i_p=1; i_p<np; i_p++)
            if (fx==1 || p[i_p][0]%fx==0) {p[iip][0]=p[i_p][0]*f[i_f]; p[iip][1]=p[i_p][1]; iip++;}
        np = iip;
    }
 
    int count_pairs=0;
    for (int i_p=np-1; i_p>0; i_p--) {
        uintll u = k/p[i_p][0] + k/p[i_p][1], a=p[i_p][0]*u, b=p[i_p][1]*u;
        cout << a << ' ' << b << '\n' << b << ' ' << a << '\n'; count_pairs+=2;
    }
    uintll u = k/p[0][0] + k/p[0][1], a=p[0][0]*u, b=p[0][1]*u;
    cout << a << ' ' << b << '\n'; count_pairs++;
    cout << "Total: "<<count_pairs<<"\n";
 
    system("pause");
    return 0;
}

    uintll k_=k; int i_f=0, i=1;
    while (k_>1 || i*i<=k) {i++; if(k_%i==0) {f[i_f++]=i; k_/=i; i=1;}} // факторизация
    if(k_>1) f[i_f++]=k_;
    int nf=i_f;

uintll k, f[30], p[3000][2]; // размеры массивов от балды, можно увеличить если надо

while (k_>1 || i*i<=k)

while (k_>1 || i*i<=k_)

    p[0][0]=1; p[0][1]=1; int np=1;
    uintll fx = 1;
    for (int i_f=0; i_f<nf; i_f++) {
        fx = (i_f==0 || f[i_f]!=f[i_f-1]) ? 1 : fx*=f[i_f];
        int iip = np;
        for (int i_p=0; i_p<np; i_p++) {
            if (fx==1 || p[i_p][1]%fx==0) {p[iip][0]=p[i_p][0]; p[iip][1]=p[i_p][1]*f[i_f]; iip++;}
            if (i_p==0) continue; 
            if (fx==1 || p[i_p][0]%fx==0) {p[iip][0]=p[i_p][0]*f[i_f]; p[iip][1]=p[i_p][1]; iip++;}
        }
        np = iip;
    }

int _tmain(int argc, _TCHAR* argv[]) 
{
    uintll k, f[30], p[3000][2]; // размеры массивов от балды, можно увеличить если надо
    //cout << "Input k: "; cin >> k;
 
    freopen("TESTS2.txt", "r", stdin);
    freopen("output.txt", "w", stdout);
 
    while (cin >> k) {
        uintll k_=k; int i_f=0, i=1;
        while (k_>1 || i*i<=k_) {i++; if(k_%i==0) {f[i_f++]=i; k_/=i; i=1;}} // факторизация
        if(k_>1) f[i_f++]=k_;
        int nf=i_f;
        //for (int i_f=0; i_f<nf; i_f++) cout<<f[i_f]<<'\n';
        //cout << "--------------------------------------------------------------------\n";
 
        p[0][0]=1; p[0][1]=1; int np=1;
        uintll fx = 1;
        for (int i_f=0; i_f<nf; i_f++) {
            fx = (i_f==0 || f[i_f]!=f[i_f-1]) ? 1 : fx*=f[i_f];
            int iip = np;
            for (int i_p=0; i_p<np; i_p++) {
                if (fx==1 || p[i_p][1]%fx==0) {p[iip][0]=p[i_p][0]; p[iip][1]=p[i_p][1]*f[i_f]; iip++;}
                if (i_p==0) continue; 
                if (fx==1 || p[i_p][0]%fx==0) {p[iip][0]=p[i_p][0]*f[i_f]; p[iip][1]=p[i_p][1]; iip++;}
            }
            np = iip;
        }
        cout << np*2-1 << "\n";
        for (int i_p=np-1; i_p>0; i_p--) {
            uintll u = k/p[i_p][0] + k/p[i_p][1], a=p[i_p][0]*u, b=p[i_p][1]*u;
            cout << a << ' ' << b << '\n' << b << ' ' << a << '\n';
        }
        uintll u = k/p[0][0] + k/p[0][1], a=p[0][0]*u, b=p[0][1]*u;
        cout << a << ' ' << b << '\n';
    }
 
    system("pause");
    return 0;
}