Combination list

 Crv ( k, v, d ) = {
    if( d == k,
        print ( vecextract( v , "2..-2" ) )
    ,   

        for( i = v[ d + 1 ] + 1, #v,
            v[ d + 2 ] = i;
            Crv( k, v, d + 1 ) ));
} 
combRV( n, k ) = Crv ( k, vector( n, X, X-1), 0 );
\\--------------------------------------------------------------

Cr ( c, z, b, n, k ) = {
    if( z < b,  print1( c, " " );
        if( n>0, Cr( c+1, z  , b*  k  \n, n-1, k - 1 ))
    ,
        if( n>0, Cr( c+1, z-b, b*(n-k)\n, n-1, k     ))
    );
}

combR( n, k ) = { 
    local( 
        bnk = binomial( n,   k ),
        b11 = bnk * k \ n );          \\binomial( n-1, k-1 ) );
    for( z = 0,  bnk - 1,
        Cr( 1, z, b11, n-1, k-1 );
        print
    );
}
\\--------------------------------------------------------------

Ci( z, b, n, k ) = { local( c = 1 );
    n--; k--;
    while( k >= 0 ,
        if( z < b,
            print1(c, "-");
            c++;
            if( n > 0, 
                b = b*k \ n);
            n--; k--;
        ,
            c++;
            z -= b;
            b = b*(n-k)\n;
            n--
        )
    );
}

combI( n, k ) = { 
    local(  bnk = binomial( n, k ), 
            b11 = bnk * k \ n );            \\ binomial( n-1, k-1 )
    for( z = 0, bnk - 1,
        Ci(z,   b11,   n, k );
        print
    );
}

Dzielniki liczby (naiwnie)

 sumaD=0; liczbaD = 0; 

dzielniki( n ) = { local( j, d=2, g=n );

sumaD = 1; liczbaD = 1;

while( d < g, 

j = n\d; 

if( d * j == n,

sumaD   += d + j; 

liczbaD += 2 );

d++ ;  g = j

);

if( g*j == n,  \\ n == x^2

liczbaD--; sumaD -= j )

};

Piękna myśl

Jeśli chcesz się nauczyć - czytaj i słuchaj,
Jeśli chcesz coś zrozumieć - pisz o tym,
Jeśli chcesz osiągnąć mistrzostwo - nauczaj o tym.

piękna myśl na dziś.

ludzie podchodzą do tego od wszystkich czterech możliwych stron
tylko czterech?

jedno chciałem zrozumieć, a tylko słuchałem i czytałem
czegoś chciałem się nauczyć, tylko pisałem o tym i w tym
chciałem być mistrzem, taaa nauczałem,

nie silnia ale klika (grafy)

można by zapytać OEIS ile jest dróg, ,,,,,,,

najmłodsza cyfra silni

int ld(int n){
        int q, t, x, z, ai;
        q=0; t=0; x=0;
        if( ((n%5)&1)==0 ) t=n%5;
        n /= 5;
        while(n){
                ai=n%5; n/=5;
                q+=ai; x+=q;
                if((ai&1)==0)
                        t += ai; }
        z = (x+t/2)%4; 
        if( z==0 ) return 6;
        else return 1<<z;}

binomial

unsigned long binomial (int n, int m) { 
        unsigned long cnm = 1UL; 
        int i, f; 
        if (m*2 >n) 
                m = n-m; 
        for (i=1 ; i <= m; n--, i++) { 
                if ((f=n) % i == 0) 
                        f /= i; 
                else 
                        cnm /= i; 
                cnm *= f; } 
        return cnm; }

Test youtube

Sito Eratostenesa

#include <iostream>
#include <bitset>
using namespace std;
 
#define N 590000000
bitset <N/2+1> t;
 
int main(){
        int i,j,k;
        printf("N=%d\n", N);
        for( i=1,k=3; k*k<N; k+=2,i++) 
                if( !t.test(i) )
                        for( j=(i*i+i)*2; j+j<N; j += 2*i+1 ) 
                                t.set(j);
        cout << t.size()-t.count()-1;
        return 0;
}