Декартово произведение

#include "stdafx.h"
#include <iostream>
using namespace std;
int _tmain(int argc, _TCHAR* argv[])
{
    setlocale (LC_ALL, "rus");
    int **dArray;
    int k, n, l, i, j, p;
    cout<<"Введите количество множеств: "; cin>>k;
    int *cArray=new int[k];
    dArray=new int*[k];
    l=0;
    for (i=0; i<k; i++)
    {
        cout<<"Количество элементов в "<<i+1<<" множестве: "; cin>>n;
        cArray[l]=n;
        l++;
        dArray[i]=new int[n];
    }
    cout<<endl<<endl<<endl;
    l=0;
    for (i=0; i<k; i++)
    {
        cout << "\nВведите "<<i+1<<" множество состоящее из "<<cArray[l]<<" неповторяющихся элементов!"<<endl;
        for (j=0; j<cArray[l]; j++)
        {
            cout<<j+1<<"-й элемент множества "<<i+1<<": "; cin>>p;
            dArray[i][j]=p;
        }
        l++;
    }
    l=0;
    cout<<endl<<endl;
    for (i=0; i<k; i++)
    {
        cout<<i+1<<" множество: A"<<i+1<<"= {";
        for (j=0; j<cArray[l]; j++)
        {
            if(j==cArray[l]-1)
            {
                cout<<dArray[i][j]<<"}"<<endl;
                break;
            }
            cout<<dArray[i][j]<<" , ";
        }
        l++;
        cout<<endl<<endl;
    }
    cout<<endl<<"Декартово (прямое) произведение множеств: D={ "; 
    // здесь должен быть вывод элементов в виде декартового произведения
    cout<<" }"<<endl<<endl;
    system("pause");
    return 0;
}