因数分解C++

Description
求出区间[a,b]中所有整数的质因数分解。
Input
输入两个整数a，b。
Output
每行输出一个数的分解，形如k=a1a2a3…(a1<=a2<=a3…，k也是从小到大的)(具体可看样例)
Sample Input
3 10
Sample Output
3=3
4=22
5=5
6=23
7=7
8=222
9=33
10=25
Hint
提示
先筛出所有素数，然后再分解。
数据规模和约定
2<=a<=b<=10000

#include <iostream>
using namespace std;
int isprime(int n) {
    if(n <= 1)
        return 0;
    else if(n == 2 || n == 3)
        return 1;
    else {
        for(int i = 2; i * i < n; i++) {
            if(n % i == 0) {
                return 0;
            }
        }
        return 1;
    }
}
int main() {
    int a, b;
    cin >> a >> b;
    for(int i = a; i <= b; i++) {
        int temp = i;
        cout << i << "=";
        int flag = 0;
        while(temp != 1) {
            for(int j = 2; j <= temp; j++) {
                if(isprime(j) && temp % j == 0) {
                    temp = temp / j;
                    if(flag == 1)
                        cout << "*";
                    cout << j;
                    flag = 1;
                    break;
                }
            }
        }
        cout << endl;
    }
    return 0;
}

太妙了吧，怎么想出来的，又是膜拜大佬的一天
（答案参考：https://github.com/liuchuo/Lanqiao）