Возвращает один и тот же результат. Найти произведение ряда

#include<iostream.h>
#include<conio.h>
#include<math.h>
#include<iomanip.h>
 
long double f(int n)
{ if(n<0) return 0;
  if (n==0) return 1;
  else return n*f(n-1);}
 
void main()
{ int x,k=0,n;
  float r;
  cout<<"Enter upper border: ";
  cin>>n;
  cout<<"Enter X: ";
  cin>>x;
  while(1<=k && k<=n) {r*=(pow((1-x),k+1)+1)/((pow(f(k-1),2)+1));}
  cout<<"Result: "<<setprecision(4)<<r;
  getch();
  }

r=1.0;  // !!!!!!!!!!!!!!!!!!!!!!!
double pk=(1.0-x)*(1.-x);
double kfact=1.0;
int k=1;
do
{
    ak=(pk+1.0)/( (kfact+1.0)*(kfact+1.0) );
    r+=ak;
 
    pk*=(1.0-x);
    kfact*=k;
    k++;
}while(k<=n);

#include<iostream.h>
#include<conio.h>
#include<math.h>
#include<iomanip.h>
 
/*long double f(int n)
{ if(n<0) return 0;
  if (n==0) return 1;
  else return n*f(n-1);}*/
 
 
 
void main()
{ int x,n;
  float r,ak;
 
  cout<<"Enter upper border: ";
  cin>>n;
  cout<<"Enter X: ";
  cin>>x;
  r=1.0;
  double pk=(1.0-x)*(1.0-x);
  double kfact=1.0;
  int k=1;
  while (k<=n)
{
  ak=(pk+1.0)/((kfact+1.0)*(kfact+1.0));
  r*=ak;
 
  pk*=(1.0-x);
  kfact*=k;
  k++;
}
  cout<<"Result: "<<setprecision(4)<<r;
  getch();
  }

void main()
{
  cout<<"Enter upper border: ";
  int n;cin>>n;
  cout<<"Enter X: ";
  int x;cin>>x;
  double r=1.0;  // !!!!!!!!!!!!!!!!!!!!!!!
  double pk=(1.0-x)*(1.-x);
  double kfact=1.0;
  int k=1;
  do
  {
    ak=(pk+1.0)/( (kfact+1.0)*(kfact+1.0) );
    r+=ak;
 
    pk*=(1.0-x);
    kfact*=k;
    k++;
  }while(k<=n);
  cout<<"Result: "<<setprecision(4)<<r;
  getch();
}