#include<iostream>
#include<cmath>
using namespace std;

int d;
double v[16];

double interpolate(int n,int e);
bool outlier(int e);

int main(){
  while(cin>>d && d){
    for(int i=0;i<d+3;++i) cin >>v[i];
    for(int i=0;i<d+3;++i){
      if(outlier(i)){
	cout<<i<<endl;
	break;
      }
    }
  }
}

double interpolate(int n,int e){
  double ln,numer,denom;
  ln=0.0;
  for(int i=0;i<d+3;i++){
    numer=1.0;
    denom=1.0;
    if(i!=n&&i!=e){
      for(int j=0;j<d+3;j++){
	if(j!=n&&j!=e&&j!=i){
	  numer*=n-j;
	  denom*=i-j;
	}
      }
      ln+=numer*v[i]/denom;
    }
  }
  return ln;
}

bool outlier(int e){
  for(int i=0;i<d+3;i++){
    if(i==e) continue;
    double p = interpolate(i,e);
    if(fabs(p-v[i])>0.1) return false;
  }
  return true;
}