Решение задачи "Разложение на простые множители" с Меньшиков

#include <iostream>
#include <vector>
#include <cmath>
 
using namespace std;
int n;
vector<int> muls;
void FindAllSumMul()
{
    int curMul = 2;
    int sqrt_N = sqrt((double)n);
    while (curMul<=sqrt_N)
    {
        while (n % curMul == 0)
        {
            muls.push_back(curMul);
            n /= curMul;
        }
        curMul++;
        sqrt_N = sqrt((double)n);
    }
    if (n!=1)
        muls.push_back(n);
}
void output()
{
    cout<<muls[0];
    for (int i=1;i<muls.size();i++)
        cout<<'*'<<muls[i];
}
int main()
{
    cin>>n;
    FindAllSumMul();
    output();
    return 0;
}