[Programare]Suma patratelor si a cuburilor

[Hertz]

Patratele primelor 5 numere intregi pozitive sunt 1,4,9,16,25

Cuburile primelor 5 numere intregi pozitive sunt 1,8,27,125

Considerati acum numerele care pot fii scrise ca suma unui patrat si unui cub.Cateva exemple :

9=1+8=1^2+2^3

12=4+8=2^2+2^3

232=16+216=4^2+6^3

Gasiti suma tuturor numerelor distincte mai mici decat 100.000 care pot fii scrise ca suma unui patrat si unui cub.

[AlStar]

1168607 ?

#include <iostream.h>
#include <fstream.h>
#include <math.h>
void main()
{float a=1,b=1;
 long double s=0;
 ofstream f("suma.txt");
 int c=0;
 for (a=1;c==0;a++)
      for (b=1;c==0;b++)
	   if (pow(a,2)+pow(b,3)<=100000)
	       s+=pow(a,2)+pow(b,3);
	     else c=1;
 f<<s;
}

Edited December 5, 2009 by AlStar

[begood]

nu-i nevoie sa folosesti fstream aici. sau folosesti iostream, sau fstream.

[AlStar]

Stiam, dar n-am fost atent cand am postat...si mi-era lene sa mai modific..

[bogyor]

mda interesant le iei din fisier

[AlStar]

Nu. Numarul il pun in fisier.

Nu era obligatorie includerea fisierului, dar e mai usor de copiat rezultatul.

[bogyor]

#include <cstdlib>

#include <iostream>

using namespace std;

int main(int argc, char *argv[])

{

int n,a[100],pozmax,pozmin,i,l,min,max,k;

cout<<"nr de elemente=";cin>>n;

min=999;max=0;

for(i=0;i<n;i++)

{

cout<<"A["<<i<<"]=";

cin>>a;

}

for(k=0;k<n;k++)

if(a[k]<min )

{

min=a[k];

pozmin=k+1;

}

cout<<"minimul sirului este = "<<min<<endl;

cout<<"minimul sirului se alfa pe pozitia a "<<pozmin<<" a"<<endl;

for(l=0;l<n;l++)

if(a[l]>=max)

{

max=a[l];

pozmax=l+1;

}

cout<<"maximul sirului este = "<<max<<endl;

cout<<"maximul sirului se alfa pe pozitia a "<<pozmax<<" a"<<endl;

[AlStar]

ce-i cu programul ala? aci se vorbeste despre altceva...

[tw8]

@AlStar:

In primul rand, programul tau are o greseala evidenta, si anume considera doar numerele care pot fi scrise ca suma dintre 1 si un cub, din cauza conditiei de oprire .

In al doilea rand, chiar daca ii repari aceasta greseala - foarte usor de reparat - mai e ceva. Iei de mai multe ori in vedere aceleasi numar, in cazul in care poate fi scris in mai multe moduri ca suma de cub si patrat. Spre exemplu, numarul 17 poate fi scris ca :

16 + 1 = 4^2 + 1^3

9 + 8 = 3^2 + 2^3

Si, nu in ultimul rand, ai un cod foarte neeficient din punct de vedere al timpului de executie. Un cod mult mai rapid - si functional - este:

#include <iostream>
#define DIM 100000
bool trecut[DIM];
int main()
{
int P=1,C,j,s=0;
for(int i=1;P<DIM;P+=2*i+1,i++)
	for(C=j=1;P+C<DIM;C+=3*j*(j+1)+1,j++)
		if(!trecut[P+C]){trecut[P+C]=true,s+=P+C;}
std::cout<<s;
}

In concluzie ... codul de mai sus ruleaza instantaneu chiar si pentru numere de pana la 10 milioane - desigur, daca s este definit ca long long in loc de int cum am eu in program, ca sa poate intra in memorie - si in plus, da si raspunsul corect: 502682822.

Bafta!