Найти пары простых чисел

#include <QCoreApplication>
#include <iostream>
 
using namespace std;
 
int main()
{
    unsigned long j,p,q, result;
 
    for (p=q=j=1000; j<10000; j++)
        if (p%j!=0 && q%j!=0 && (p+1)/q!=0)
        {
            p+=j;
            q+=j;
            result = (p+1)/q;
            cout << p <<" " << q << "Result: " << result << endl;
        }
 
    return 0;
}

#include <iostream>
#include <cmath>
using namespace std;
 
bool is_prime(unsigned long p)
{
    if (p == 2)
        return true;
    if (p % 2 == 0)
        return false;
 
    unsigned long lim = sqrt(p);
    for (int delim = 3; delim < lim; delim += 2)
    {
        if (p % delim == 0)
            return false;
    }
    return true;
}
 
int main()
{
    unsigned long lower_bound = 1000, upper_bound = 10000;
    unsigned long p, q;
    for (p = lower_bound; p < upper_bound; ++p)
    {
        if (is_prime(p))
        {
            for (q = lower_bound; q < upper_bound; ++q)
            {
                if (((p + 1) % q == 0) && is_prime(q))
                {
                    cout << p << " " << q << endl;
                    break;
                }
            }
        }
    }
    return 0;
}