Вывести последовательность элементов, сумма которых будет простым числом.

int sk; // Используется как битовая шкала:i-тый бит= 1 - M[i] входит в сумму
int Max = (1<<n); 
for(sk=1; sk<Max; sk++) {
  for(i=sum=0; i<n; i++)
    if (sk&(1<<i)) sum+=M[i];
  if (Prostoe(sum)) {
   for(i=sum=0; i<n; i++)
     if (sk&(1<<i)) cout<<"+" <<M[i];
   cout<< "="<<sum<<endl;
  }
}

bool isProstoe(int number)
{  int M;
  if (number<=1) return false;
  if (number==2) return true;
  M = sqrt(n);
  for (int i=3; i<M; i+=2)
    if (number % i == 0) return true; 
  return false;
}

#include <stdio.h>
 
#define N  8
 
int buf[1000];
int bufp;
int sum;
 
void process_term(int pos, int *arr, int size);
 
int main()
{
    /* 
     * There is a one trick: the sequence is preceded by zero,
     * fake-element, which are the other terms processed by. 
     */
    int arr[N] = {0, 1, 2, 2, 3, 2, 4, 1};
    int i, b;
    
    printf("Sequence:");
    for (i = 1; i < N; i++)
        printf(" %d", arr[i]);
    putchar('\n');
 
    process_term(0, arr, N);
    
    for (i = 0; i < bufp; i++) {
        for (b = 1; b < N; b++) {
            if (buf[i] & (1 << b)) 
                printf("%d ", arr[b]);
        }
        printf("\n");
    }
    return 0;
}
 
int isprime(int);
 
void process_term(int pos, int *arr, int size)
{
    int localbuf, localsum;
 
    buf[bufp] |= 1 << pos;
    sum += arr[pos];
 
    localbuf = buf[bufp];
    localsum = sum;
 
    if (isprime(sum) && sum != arr[pos]) {
        bufp++;
        buf[bufp] = buf[bufp - 1];
    }
    for (pos += 1; pos < size; pos++) {
        process_term(pos, arr, size);
        sum = localsum;
        buf[bufp] = localbuf;
    }   
}
 
int isprime(int n)
{
    int i;
 
    for (i = 2; i < n; i++)
        if (n % i == 0)
            return 0;
    return 1;
}

bool isProstoe(int number)
{  int M;
  if (number<=1) return false;
  if (number==2) return true;
  M = sqrt(number);
  for (int i=3; i<=M; i+=2)
    if (number % i == 0) return false; 
  return true;
}