example_funcmatrix_indexsum.cpp

 
#include <Pds/Ra>
    
int main(void)
{
    std::cout<<"\n((3,3)): "<<Pds::NmultichooseK(3,3)<<std::endl;
    Pds::Matrix J = Pds::MultisetIndexSum(3,3);
    J.Print("J=MultisetIndexSum(3,3):\n");
    
    std::cout<<"\n((3,2)): "<<Pds::NmultichooseK(3,2)<<std::endl;
    Pds::Matrix K = Pds::MultisetIndexSum(3,2);
    K.Print("K=MultisetIndexSum(3,2):\n");
    
    std::cout<<"\n((3,1)): "<<Pds::NmultichooseK(3,1)<<std::endl;
    Pds::Matrix L = Pds::MultisetIndexSum(3,1);
    L.Print("L=MultisetIndexSum(3,1):\n");
    
    std::cout<<std::endl;
    
    std::cout<<"\n(3,3): "<<Pds::NchooseK(3,3)<<std::endl;
    Pds::Matrix A = Pds::SetIndexSum(3,3);
    A.Print("A=SetIndexSum(3,3):\n");
    
    std::cout<<"\n(3,2): "<<Pds::NchooseK(3,2)<<std::endl;
    Pds::Matrix B = Pds::SetIndexSum(3,2);
    B.Print("B=SetIndexSum(3,2):\n");
    
    std::cout<<"\n(3,1): "<<Pds::NchooseK(3,1)<<std::endl;
    Pds::Matrix C = Pds::SetIndexSum(3,1);
    C.Print("C=SetIndexSum(3,1):\n");
    
    std::cout<<"\n(4,2): "<<Pds::NchooseK(4,2)<<std::endl;
    Pds::Matrix D = Pds::SetIndexSum(4,2);
    D.Print("D=SetIndexSum(4,2):\n");
    
    return 0;
}