Корни квадратного уравнения

if ( (корень - static_cast<int>(корень) ) == 0 )
     да
else нет

#include<iostream>
#include<iomanip>
#include<cmath>
using namespace std;
//include<
int main()
{
    //x^2 - p[i]x + q = 0;
    double D;
    const int n=7;
    double root1;
    double root2;
    const double q=4.0;
    double p[n] = { 2.6, 3.3, 0.8, 5.6, 0.5, -2.8, -4.2 };
    double r[n] = {0};
    int i=0, j=0;
    cout<<"x^2-p[i]x+q = 0\nq="<<q<<"\n\n";
    cout<<"Vector p:\n";
    for ( int i=0; i<n; i++ )
        cout<<p[i]<<"   ";
    cout<<"\n\nVector r:\n";
    for ( i=0, j=0; i<n; i++ )
    {
        D = pow( p[i], 2.0 ) - (4.0*q);
        if ( D < 0 )
            continue;
        root1 = ( p[i]+sqrt(D) ) / 2.0;
        root2 = ( p[i]-sqrt(D) ) / 2.0;
        if ( ( root1!=root2 ) )
        {
            r[j] = p[i];
            cout<<r[j]<<"   ";
            j++;            
        }
    }
    cout<<endl;
    system("pause");
    return 0; // ретурн отакэ
}