//#include "stdafx.h"
#include<conio.h>
#include<iostream.h>
#include<fstream.h>
#include<process.h>

template<class TIP> class Vector
{ 
 TIP* pv;
   int dim;
public:  
	  Vector() {pv=NULL;dim=0;}         //constructori de clasa
	  Vector(int);
	  Vector(TIP,int);
	  Vector(TIP*,int);
	  Vector(Vector &);                 //constr.de copiere
	  ~Vector() {delete[] pv;}          //destructor
	  Vector& operator=(Vector);
	  int operator!() {if(dim) return 1; //test vector vid
						else return 0;}
	  Vector operator+(Vector);         //concatenare a 2 vectori
	  Vector& operator+=(Vector);    
	  int operator==(Vector);           //test de egalitate a 2 vectori
	  operator TIP() {return pv[0];}    //cast catre TIP
	  TIP& operator[](int i) {return pv[i];}  //referirea unui element

bool cresc(TIP,TIP);
bool desc(TIP,TIP);
void comuta(TIP&a,TIP&b) {TIP aux;aux=a;a=b;b=aux;}
void sortare(bool(*sens)(TIP,TIP));         //sortarea dupa o directie primita ca pointer de functie
int regresie(Vector&,Vector&,double&,double&);  //coeficientii dreptei de regresie
      
	  friend ostream& operator<<(ostream&,Vector);  
	  friend istream& operator>>(istream&,Vector&);
	  friend ifstream& operator>>(ifstream&,Vector&);
};


template<class TIP>Vector<TIP>::Vector(int n)
{ pv=new TIP[n];
  dim=n;
}

template<class TIP>Vector<TIP>::Vector(TIP in,int n)
{ pv=new TIP[n];
  dim=n;
  for(int i=0;i<dim;i++)
     pv[i]=in;     						
}

template<class TIP>Vector<TIP>::Vector(TIP*ini,int n)
{ pv=new TIP[n];
  dim=n;
  for(int i=0;i<dim;i++)
     pv[i]=ini[i];     						
}

template<class TIP>Vector<TIP>::Vector(Vector &ve)
{ if(ve.pv)
	{ dim=ve.dim;
	pv=new TIP[dim];
  for(int i=0;i<dim;i++) pv[i]=ve.pv[i];
	}
  else pv=NULL;
}

template<class TIP>Vector<TIP>& Vector<TIP>::operator=(Vector ve)
{if((pv!=ve.pv)&&(ve.pv!=NULL))
	{delete[] pv;
    dim=ve.dim;
	pv=new TIP[dim];
  	for(int i=0;i<dim;i++) pv[i]=ve.pv[i];
	}   
return *this;
}

template<class TIP>Vector<TIP> Vector<TIP>::operator+(Vector ve)
{ Vector<TIP> re(dim+ve.dim);
  re.dim=dim+ve.dim;
  re.pv=new TIP[re.dim];
	if(pv==NULL) re=ve;
	 else if(ve.pv==NULL) re=*this;
	  else
	  {
       for(int i=0;i<dim;i++) re.pv[i]=pv[i];
       for(i=dim;i<re.dim;i++) re.pv[i]=ve.pv[i-dim];
	  }
return re;
}

template<class TIP>Vector<TIP>& Vector<TIP>::operator+=(Vector ve)
{
	if(ve.dim==1)
	{  Vector aux(dim);
	   aux=*this;
	   *this=*this+ve;
	   return *this;
	}
}

template<class TIP>int Vector<TIP>::operator==(Vector v)
{if(dim!=v.dim) return 0;
   else 
   { int cont=0;
     for(int i=0;i<dim;i++) if(pv[i]==v.pv[i]) cont++;
     if(cont==dim) return 1;
                 else return 0;
   }
}

template<class TIP>bool cresc(TIP a,TIP b)
{
	return a<b?1:0;
}

template<class TIP>bool desc(TIP a,TIP b)
{
	return a>b?1:0;
}

template<class TIP>void Vector<TIP>::sortare(bool(*sens)(TIP,TIP))
{ int i,j;
for(i=0;i<dim-1;i++)
 for(j=i+1;j<dim;j++)
   if(!(*sens)(pv[i],pv[j])) comuta(pv[i],pv[j]);
}

template<class TIP>int Vector<TIP>::regresie(Vector<TIP>&v1,Vector<TIP>&v2,double&a0,double&a1)
{ 	double mx,my;
    mx=my=0;
   for(int i=0;i<v1.dim;i++)
	{
	  mx+=v1.pv[i];
      my+=v2.pv[i];
	}
    int rez=1;
    if(v1.dim)
    {
	  mx/=v1.dim;
      my/=v2.dim;
    }
    else rez=0;
   double cov,disp;
   cov=disp=0;
   for(i=0;i<v1.dim;i++)
    {
	 cov+=(v2.pv[i]-my)*(v1.pv[i]-mx);
     disp+=(v1.pv[i]-mx)*(v1.pv[i]-mx);
    }
    if (disp)
	{
	 a1=cov/disp;
     a0=my-a1*mx;
    }
    else rez=0;

return rez;
}

template<class TIP>ostream& operator<<(ostream& ies,Vector<TIP> v)
{  if((v.dim!=0)&&(v.pv!=NULL)) 
	{ ies<<"dim.vect: ";ies<<v.dim<<"  ";
      ies<<" Vectorul: ";
      for(int i=0;i<v.dim;i++) ies<<v.pv[i]<<" ";
      ies<<endl;
         return ies;
	}
 //if(v.dim==0)
 else
 { ies<<"dim vect: ";ies<<v.dim<<"  ";
   for(int i=0;i<v.dim;i++) ies<<v.pv[i]<<" ";
   ies<<endl;
       return ies;
 }
}

template<class TIP>istream& operator>>(istream& intr,Vector<TIP> &v)
{
    if(v.dim==0) 
	{cout<<"dim.vect: ";intr>>v.dim;
	 v.pv=new TIP[v.dim];
	 for(int i=0;i<v.dim;i++) intr>>v.pv[i];
	   return intr;
	}
      else {
		  for(int i=0;i<v.dim;i++) intr>>v.pv[i];
          return intr;
			}
}

template<class TIP>ifstream& operator>>(ifstream& intr,Vector<TIP>& v)
{ //if((v.dim!=0)&&(v.pv!=NULL))
	 for(int i=0;i<v.dim;i++) intr>>v.pv[i];
      return intr;
	
 }

int main(int argc, char* argv[])
{
	Vector<int>vi(3,2),ve,vu,va;
	cout<<endl<<"Pentru testarea lucrului cu vectori de intregi"<<endl<<endl;
    cin>>ve;
    cout<<"Deci vectorul citit este:  "<<ve;
	va=vi;
	cout<<"Un alt vector este   "<<va;
	vu=ve+vi;
	cout<<"Operatorul +(concatenare)  "<<vu;
	cout<<"Referirea unui element din vectorul introdus(1.elem)  "<<ve[0]<<endl;
	cout<<"Cast catre int(primul elem.din vector)  "<<int(ve)<<endl;
    if(ve==va) cout<<"Cei doi vectori sunt identici!"<<endl;
	        else cout<<"Cei doi vectori nu sunt egali!"<<endl;
    
 
    cout<<endl<<"Pentru testarea lucrului cu vectori de float"<<endl<<endl;
	float v[]={-1,0,1.1};
	Vector<float>vf(1),vh(2),vp,vl,vr(1,1),vg(v,3),vd,vk;

	cout<<"Introduceti vectorul cu variabila cauza:"<<endl;
	cin>>vp;
	cout<<"Introduceti vectorul cu variabila efect:"<<endl;
	cin>>vd;
	int er,ereg;
	double re0,re1;
	er=vp.regresie(vp,vd,re0,re1);
	if(er)
	  cout<<"Dreapta de regresie este  Y^="<<re0<<"+"<<re1<<"*X"<<endl<<endl;
	  else cout<<"Eroare!"<<endl;

	vl=vg;
	cout<<"Un alt vector este  "<<vl;
	cout<<"Iar celalalt este  "<<vg;
	ereg=vl==vg;
	if(ereg) cout<<"Acesti vectori sunt identici!"<<endl;
	  else cout<<"Acesti vectori nu sunt egali!"<<endl;

    cout<<"Operatorul + ptr.acesti 2 vectori "<<vl+vg;

    	float t[]={1.2,2.3};
	Vector<float>vo;
    vk=vr+vo+vp;
    vp+=vr;
	cout<<"Op. += ptr.vect.cauza si un vect.de 1 elem ";cout<<vp;
 
	cout<<"Cast catre float(1.elem din acest vect)  "<<float(vp);
	if(!vp) cout<<"   si in plus vectorul NU este vid!"<<endl<<endl;
	  else cout<<"In plus vectorul este VID!"<<endl;
    

//scriere in fisier pe disc
ofstream fisout("VECT.DAT");
	if(!fisout) {cout<<"Eroare desch.fis!";exit(1);}
	fisout<<vr<<vg<<vp;
	fisout.close();

//citire obiect din fisier de pe disc
	ifstream fisin("VECT.DAT");
	if (!fisin) {cout<<"Fis negasit!";exit(2);}
	 cout<<"Citire din fisier:"<<endl;
	 fisin>>vr;
	 cout<<vr;
	 fisin>>vg;
	 cout<<vg;
	fisin.close();

//sortare dupa o anumita directie 
	bool (*pf)(float,float);
	pf=desc;
	cout<<"Acest vector sortat descresc.este  ";
	vg.sortare(pf);
	cout<<vg;

	return 0;
}