돌아오지 않는 숲

종종 통계물리 연구를 하다보면 종종 여러 파일에 대해서 명령어를 연속적으로 입력해 야 할때가 있습니다. 예를 들어서 

대부분의 사람들은 통계 물리에 대해 막연하게 생각하는 부분이 많다..

/* made by JEong Hyang Min  (blindfish)
 * 2008. 06. 22*/

#include<stdio.h>
#include<stdlib.h>
#include<time.h>
#include<math.h>

int main()
{
    int i,j;
    int node;
    int link;
    int num1=0;
    int num2=0;
    
    int **network;

    printf("input number of node===>");
    scanf("%d", &node);
    printf("\n");
    printf("input number of link===>");
    scanf("%d", &link);
    printf("\n");

    network = (int **)malloc(sizeof(int*) * node);

    for(i = 0; i < node; i++)
        network[i] = (int *)malloc(sizeof(int) * node);

    for(i=0; i<node; i++)
        for(j=0; j<node; j++)
            network[i][j]=0;

    srand(times(NULL));

    for(i=0; i<=link; i++)
    {
        num1=rand()%node;
        num2=rand()%node;

        if(network[num1][num2]!=0 && network[num2][num1]!=0 && num1!=num2)
            i--;
        else
        {
            network[num1][num2]=1;
            network[num2][num1]=1;
        }
    }
    for(i=0; i<node; i++)
    {
        for(j=0; j<node; j++)
            printf("%d", network[i][j]);
        printf("\n");
    }

    free(network);
    return 0;
}

/* made by JEong Hyang Min  (blindfish)
 * 2008. 07. 08*/
#include<stdio.h>
#include<stdlib.h>
#include<time.h>
#include<math.h>

int main()
{
    int i,j,k;
    int node=0;
    int link=0;
    int num1=0;
    int totsum=0;
    int *sum;
    int *prob;

    int **network;

    printf("input number of node===>");
    scanf("%d", &node);
    printf("\n");
    printf("input number of link===>");
    scanf("%d", &link);
    printf("\n");

    network = (int **)malloc(sizeof(int*) * node);
    sum = (int *)malloc(sizeof(int) * node);
    prob = (int *)malloc(sizeof(int) * node);
    srand(times(NULL));

// 처음 행렬을 0을 채우고 노드갯수만큼의 행렬을 생성한다.
    for(i=0; i<node; i++)
        network[i] = (int *)malloc(sizeof(int) * node);
    for(i=0; i<node; i++)
        for(j=0; j<node; j++)
            network[i][j]=0;
//seed값을 주고 여기서 자라게 한다.
    network[0][0]=0;
    network[1][0]=1;
    network[0][1]=1;
    network[2][1]=1;
    network[1][2]=1;
//처음에 합산과 확률분포 행렬을 초기화한다.
    for(i=0; i<node; i++)
    {
        sum[i]=0;
        prob[i]=0;
    }

    for(i=3; i<node; i++)    //노드를 하나씩 추가한다.
    {
        totsum=0;
        for(j=0; j<i; j++)
        {
            for(k=0; k<j; k++)
            {
                sum[j] += network[j][k];    //노드의 링크수를 센다.
                totsum += network[j][k];    //총 링크수를 계산한다.
            }
        }

        for(j=0; j<i; j++)
        {
            for(k=0; k<=j; k++)
                prob[j] += sum[k];    //확률 분포를 만든다..
           
            printf("%d ", prob[j]);
        }
        printf("\n");

        num1=rand();    //몬테 카를로를 쓸 수를 만든다.
        num1=num1%totsum;    //확률분포안에 정규화 시킨다.
        for(j=0; j<i; j++)
        {
            if(num1<=prob[0])
            {
                network[i][0]=1;    //해당되는 확률함수를 찻고 링크와 노드를 추가시킨다.
                network[0][i]=1;
            }
            else if(prob[j]<num1 && num1<=prob[j+1])
            {
                network[i][j]=1;
                network[j][i]=1;
            }
            else
            {
                network[i][j]=0;
                network[j][i]=0;
            }
        }

        for(j=0; j<=i; j++)
        {
            sum[j]=0;    //값을 다시 초기화한다.
            prob[j]=0;
        }
    }

    for(i=0; i<node; i++)
    {
        for(j=0; j<node; j++)
            printf("%d", network[i][j]);
        printf("\n");
    }

    free(network);
    free(sum);    //동적할당에선 꼭 메모리를 free시킨다.
    free(prob);    // 안그러면 메모리 문제로 재부팅 해야할수도 있다.
    return 0;
}

척도없는 네크워크 = 정보가 다루는 계나 네트워크의 총체적 정보가 필요함..

포텐셜 문제의 real field, pow-law

매력도==> charge

random network

전체 네트워크에서의 총체적인 정보가 없어도 된다..

== mean field thoery 와 같다..