Matematyka.pl

Kod: Zaznacz cały

#include <stdio.h>
int tablica[30001];
int ok(int i)
{
    int t[5],w=0,s=i,k,suma=0,iloczyn=1;
    while(s>0)
    {
        t[w]=s%10;
        s/=10;
        w++;
    }
    for(k=0;k<w;k++)
    {
        suma+=t[k];
        iloczyn*=t[k];
    }
    return (tablica[suma]&&tablica[iloczyn]);
}
int main()
{
	int i,j,t;
	for(i=1;i<=30000;i++) tablica[i]=1;
	for (i=2; i<=30000; i++)
	{
        if (tablica[i])
   		{
      	  	j = 2*i;
          	while (j<=30000)
      	  	{
         		tablica[j] = 0;
         		j += i;
          	}
   		}
	}
	for(i=30;i<30001;i++)
        if(tablica[i]&&ok(i)) printf("%d
",i);
	return 0;
}

Kod: Zaznacz cały

   for(i=3;i<=30000;i++) tablica[i]=i & 1;
   tablica[2]=1; tablica[0]=tablica[1]=0;
   for (i=2; i<=173; i++)
        if (tablica[i])    {
             j = i*i;
             while (j<=30000) {
               tablica[j] = 0;
               j += i;
             }
         }

Kod: Zaznacz cały

#include <stdio.h>

int isprime(unsigned long long int n){    //  poprawka, dodałem "long long"
	unsigned int d;
	if (n<2) return 0;
	if (n<4) return 1;
	if (n%2==0 || n%3==0) return 0;
	d=5;
	while((unsigned long long int)d*d<=n){
		if (n%d==0) return 0;
		d+=2;
		if (n%d==0) return 0;
		d+=4;
	}
	return 1;
}

void g(int i, int p){
	unsigned long long s;
	int k,j;
	if(isprime(p+i)) 
		for(j=0; j<=i; j++) {
			s=0;
			for(k=i; k>=0; k--)
				if(k==j) s=s*10+p;
				else s=s*10+1;
			if (/*s<=30000 &&*/ isprime(s))
				printf("%I64d ", s);
		}
}

int main(){
	int i;
	for(i=0; i<18; i++) {
		if(i==0) g(0,2);
		if(i&1)  g(i, 2);
		else {
			g(i, 3);
			g(i, 5);
			g(i, 7);
		}
	}
	while(1);
}

Kod: Zaznacz cały

#include <cstdlib>
#include <iostream>
#include <ctime>

using namespace std;

const int nieparzyste[] = {1, 3, 5, 7};

bool czyPierwsza(int liczba)
{
    for(int i = 2; i*i <= liczba; i++)
        if(liczba % i == 0)
            return false;

    return true;
}


int main(int argc, char *argv[])
{
    long czas = clock();
    int licznik = 0;

    for(int i = 2; i < 20; i++)
        for(int j = 0; j < i; j++)
            if(i%2 == 0)
            {
                long long liczba = 0;
                for(int l = 0; l < i; l++)
                {
                    liczba *= 10;
                    if(l == j) liczba += 2;
                    else liczba += 1;
                } 
                if(czyPierwsza(i+1) && czyPierwsza(liczba))
                {
                    cout << liczba << endl; 
                    licznik++;
                }
            } else
                for(int k = 0; k < 4; k++)
                {
                    long long liczba = 0;
                    for(int l = 0; l < i; l++)
                    {
                        liczba *= 10;
                        if(l == j) liczba += nieparzyste[k];
                        else liczba += 1;
                    } 
                    if(czyPierwsza(i - 1 + nieparzyste[k]) && czyPierwsza(liczba))
                    {
                        cout << liczba << endl;
                        licznik++;
                    }
               }

    czas = clock() - czas;
    float s = static_cast<float>(czas) / CLOCKS_PER_SEC;
    cout << "Znalazlem " << licznik << " liczb" << endl;
    cout << "Zejelo mi to " << s << " sekund" << endl;
    system("PAUSE>>null");
    return EXIT_SUCCESS;
}

Kod: Zaznacz cały

uint64 mul_mod1(uint64 a, uint64 b, uint64 m) {
  uint64 y = (uint64)((float64)a*(float64)b/m+(float64)1/2); // floor(a*b/m)
  y = y * m; // m*floor(a*b/m) mod z
  uint64 x = a * b; // a*b mod z
  uint64 r = x - y; // a*b mod z - m*floor(a*b/m) mod z
  if ( (int64)r < 0 ) // normalization needed ?
  {
    r = r + m;
    y = y - 1; // (a*b)/m quotient, omit line if not needed
  }
  return r; // (a*b)%m remnant
}

Kod: Zaznacz cały

void g(int i, uint64 j, uint64 p){
	uint64 w;
	int k;
	if(isprime(p+i))	{
		k=i; w=p;
		while(k>0) { 
			w*=10; k--;
		}
		p--;
		while(w>p) {
			if(isprime(p+j)) {
                printf("%19I64d ", p+j);
				m++;
			}
			p*= 10;
		}
	}
}

int main(){
	unsigned int i,c;
	uint64 j;
	c=clock(); m=0;
	j=1;
	for(i=0; i<10; i++) {
		//printf("
(%d) ",i);
		if(i==0) g(0, j, 2);
		if(i &1) g(i, j, 2);
		else {
			g(i, j, 3);
			g(i, j, 5);
			g(i, j, 7);
		}
		j=j*10+1;
	}
	c=clock()-c;
	printf("
koniec
%d sztuk
%4d.%03d sekund",m,c/1000,c%100);while(1);
}

Kod: Zaznacz cały

bool czyPierwsza(long long int liczba)
{
    for(int i = 2; (unsigned long long int)i*i <= liczba; i++)
        if(liczba % i == 0)
            return false;

    return true;
}