__kernel void combination(const int n, const int k, __global int *out)
{
   size_t gid = get_global_id(0);
   int a[k];
   int index = gid*100;

   for(int i = 0; i < k; i++)
      a[i] = gid+i;

   while(index != -1)
   {
      for(int i = 0; i < k; i++)
         out[index++] = a[i];

      int j = k - 1;
      if(a[j] < n - 1)
      {
         ++a[j];
         continue;
      }

      while(a[j] - a[j-1] == 1)
         --j;

      if(j == 1)
      {
         index = -1;
         continue;
      }

      int z = ++a[j-1];

      while(j < k)
      {
         a[j] = ++z;
         ++j;
      }
   }
}