Хостинг портала RFpro.ru:
Московский хостер
Профессиональный ХОСТИНГ на базе Linux x64 и Windows x64

РАССЫЛКИ ПОРТАЛА RFPRO.RU

Лучшие эксперты данной рассылки

Гаряка Асмик Статус: Профессор Рейтинг: 6052 ∙ повысить рейтинг »

Boriss Статус: Академик Рейтинг: 2587 ∙ повысить рейтинг »

Абаянцев Юрий Леонидович aka Ayl Статус: Профессионал Рейтинг: 2205 ∙ повысить рейтинг »

Вопрос № 181258: Здравствуйте, уважаемые эксперты! Прошу Вас ответить на следующий вопрос об очередной задаче о перколяции на кубической решетке: Перколяция считается по 3 направлениям -горизонтальной, вертикальной и диагональной( ну по оси z то есть). Но есть сло...

Здравствуйте, уважаемые эксперты! Прошу Вас ответить на следующий вопрос об очередной задаче о перколяции на кубической решетке:
Перколяция считается по 3 направлениям -горизонтальной, вертикальной и диагональной( ну по оси z то есть). Но есть сложность в заполнении самой решетки : если решетка 10х10х10, то программа считает довольно быстро. Если решетка 20х20х20 то программа зависает(начинает считать КРАЙНЕ медленно) при вероятности заполнения решетки где на 23%. При решетке 50х50х50 - еще раньше.
Думаю, что проблема именно в ее заполнении.
Делала заполнение решетки, как в квадратной, но тогда программа вопще заполнять не хотела
//заполнение матрицы
int c1 = klinear*kwide*kheight*p; //число единиц
int k1 = 0;
int ofs_w, ofs_l, ofs_h;
srand (time(NULL));
while (k1<c1) {
ofs_w = rand()%kwide;
ofs_l = rand()%klinear;
ofs_h = rand()%kheight;
if (matrix[ofs_l][ofs_w][ofs_h]!=0) continue;
matrix[ofs_l][ofs_w][ofs_h] = 1;
k1++;
}

Отправлен: 12.12.2010, 22:37
Вопрос задал: Yulesik (1-й класс)
Всего ответов: 1
Страница вопроса »

Приложение:
#include "stdafx.h" #include <time.h> #include <conio.h> #include <iostream> #include <assert.h> #include <stdlib.h> using namespace std; const char FlagX=2, FlagY=3, FlagZ=4; struct State // структура для хранения рекурсивного стека { State(int x, int y, int z, char *pos, State *prev) : x(x),y(y),z(z),Prev(prev),Pos(pos) {Stage=0;} int x, y, z; int Stage; char *Pos; State *Prev; }; int X, Y, Z;//длины размерностей int yStep, zStep;//Шаг в байтах между элементами этой размерности char *Matrix; //Обслуживание создания/удаления элементов стека State *WaitList = NULL; inline State *newState(int x, int y, int z, char *pos, State* prev) { if(!WaitList) return new State(x,y,z,pos,prev); else { State *S=WaitList; WaitList=WaitList->Prev; S->x=x; S->y=y; S->z=z; S->Pos=pos; S->Prev=prev; S->Stage=0; return S; } } inline void deleteState(State *S) { S->Prev=WaitList; WaitList=S; } #define Del(S) while(S){ State *T=S; S=S->Prev; deleteState(T);} //сама программа bool FindXPercolation(); bool FindYPercolation(); bool FindZPercolation(); int main() { setlocale( LC_ALL,"Russian" ); int cntAll = 0, cntSome = 0, cntX = 0, cntY = 0, cntZ = 0; // количество перколяций float pv1,pg1,pd1, pob1,pod1; // вероятности возникновения перколяций cout<<"Vvedite dlinu, vysotu and glubinu "<<endl; cin >> X >> Y >>Z ; // m = строки, n = стобцы yStep = X; zStep = X*Y; // создаю matrix размерностью [X][Y][Z] int size = X*Y*Z; Matrix = new char[size]; if(!Matrix) { cout <<"Слишком большие размеры!"; return 1; } srand( (unsigned)time( NULL ) ); int randstep = RAND_MAX / 100; for (int p = randstep; p <= randstep*100; p += randstep) { int k=0; int cntAll = 0, cntSome = 0, cntX = 0, cntY = 0, cntZ = 0; // вывод вероятности cout<<"вероятность - "<< p/randstep*0.01 << endl; // счетчик прогонов for (int s = 0; s < 10; s++) { char *b = Matrix; for (int i = 0; i<size; i++) //Заполняем матрицу *(b++) = (rand() <= p ) ? 1 : 0; bool HaveX = FindXPercolation(); cout << "X " << (HaveX ? "да":"нет"); bool HaveY = FindYPercolation(); cout << "\tY " << (HaveY ? "да":"нет"); bool HaveZ = FindZPercolation(); cout << "\tZ " << (HaveZ ? "да":"нет") << endl; if (HaveX) cntX++; if (HaveY) cntY++; if (HaveZ) cntZ++; if (HaveX && HaveY && HaveZ) cntAll++; if (HaveX || HaveY || HaveZ) cntSome++; k++; } cout<<"Кол-во перколяций по ширине: "<< cntX<<";\nКол-во перколяций по длине: " <<cntY<<";\nКол-во перколяций по высоте: "<<cntZ<<";\nКол-во всех перколяций: "<<cntAll<<";\nХотя бы одна перколяция: "<<cntSome<<endl; cout<< "kol-vo progonov dlya veroyatnosti " << p/randstep*0.01 <<" - "<<k<<endl<<endl; pv1=(float)cntX/k; pg1=(float)cntY/k; pd1=(float)cntZ/k; pob1=(float)cntAll/k; pod1=(float)cntSome/k; cout<<"Вероятность возникновения перколяции по ширине "<<pv1<<";\nВероятность возникновения перколяции по длине "<<pg1<<";\nVВероятность возникновения перколяции по высоте "<<pd1<<";\nВероятность возникновения всех перколяций "<<pob1<<";\nВероятность возникновения хотя бы одной "<<pod1<<endl<<endl; cout<<"-------------------------------------------------------------"<<endl; } int st=0; while(WaitList) { State *S=WaitList; WaitList=WaitList->Prev; delete S; st++; } cout << "Максимальный путь: " << st << endl; getch(); return 0; } //Эти функции пытаются сделать шаг в соответствующем направлении inline State* Zup(State *S, char Flag) { if(S->z < Z-1) { char *t=S->Pos+zStep; if(*t!=0 && *t!=Flag) S = newState(S->x, S->y, S->z+1, t, S); } return S; } inline State* Zdown(State *S, char Flag) { if(S->z > 0) { char *t=S->Pos-zStep; if(*t!=0 && *t!=Flag) S = newState(S->x, S->y, S->z-1, t, S); } return S; } inline State* Yup(State *S, char Flag) { if(S->y < Y-1) { char *t=S->Pos+yStep; if(*t!=0 && *t!=Flag) S = newState(S->x, S->y+1, S->z, t, S); } return S; } inline State* Ydown(State *S, char Flag) { if(S->y > 0) { char *t=S->Pos-yStep; if(*t!=0 && *t!=Flag) S = newState(S->x, S->y-1, S->z, t, S); } return S; } inline State* Xup(State *S, char Flag) { if(S->x < X-1) { char *t=S->Pos+1; if(*t!=0 && *t!=Flag) S = newState(S->x+1, S->y, S->z, t, S); } return S; } inline State* Xdown(State *S, char Flag) { if(S->x > 0) { char *t=S->Pos-1; if(*t!=0 && *t!=Flag) S = newState(S->x-1, S->y, S->z, t, S); } return S; } //Эти функции перебирают возможные шаги (сам поиск) bool TravZ(int x, int y, int z, char *pos) { State *S = newState(x,y,z,pos,0); char Flag = FlagZ; while(S!=0) { switch(S->Stage) { case 0://Вход в функцию if(S->z == Z - 1) { Del(S); return true;} *(S->Pos) = Flag; case 1: S->Stage=2; S = Zup(S, Flag); break; case 2: S->Stage=3; S = Yup(S, Flag); break; case 3: S->Stage=4; S = Xup(S, Flag); break; case 4: S->Stage=5; S = Xdown(S, Flag); break; case 5: S->Stage=6; S = Ydown(S, Flag); break; case 6: S->Stage=7; S = Zdown(S, Flag); break; case 7://ячейка обработана State *t=S; S=S->Prev; deleteState(t); } } return false; } bool TravY(int x, int y, int z, char *pos) { State *S = newState(x,y,z,pos,0); char Flag = FlagY; while(S!=0) { switch(S->Stage) { case 0://Вход в функцию if(S->y == Y - 1) { Del(S); return true;} *(S->Pos) = Flag; case 1: S->Stage=2; S = Yup(S, Flag); break; case 2: S->Stage=3; S = Zup(S, Flag); break; case 3: S->Stage=4; S = Xup(S, Flag); break; case 4: S->Stage=5; S = Xdown(S, Flag); break; case 5: S->Stage=6; S = Zdown(S, Flag); break; case 6: S->Stage=7; S = Ydown(S, Flag); break; case 7://ячейка обработана State *t=S; S=S->Prev; deleteState(t); } } return false; } bool TravX(int x, int y, int z, char *pos) { State *S = newState(x,y,z,pos,0); char Flag = FlagX; while(S!=0) { switch(S->Stage) { case 0://Вход в функцию if(S->x == X - 1) { Del(S); return true;} *(S->Pos) = Flag; case 1: S->Stage=2; S = Xup(S, Flag); break; case 2: S->Stage=3; S = Yup(S, Flag); break; case 3: S->Stage=4; S = Zup(S, Flag); break; case 4: S->Stage=5; S = Zdown(S, Flag); break; case 5: S->Stage=6; S = Ydown(S, Flag); break; case 6: S->Stage=7; S = Xdown(S, Flag); break; case 7://ячейка обработана State *t=S; S=S->Prev; deleteState(t); } } return false; } //Эти функции обходят плоскость и запускают поиск из найденых узлов bool FindZPercolation() //Z=0 { char *b = Matrix; for(int y=0; y<Y; y++) for(int x=0; x<X; x++) { if(*b!=0 && *b!=FlagZ) if(TravZ(x, y, 0, b)) return true; b++; } return false; } bool FindYPercolation() //Y=0 { char *b; for(int z=0; z<Z; z++) { b = Matrix + z*zStep; for(int x=0; x<X; x++) { if(*b!=0 && *b!=FlagY) if(TravY(x, 0, z, b)) return true; b++; } } return false; } bool FindXPercolation() //Y=0 { char *b; for(int z=0; z<Z; z++) { b = Matrix + z*zStep; for(int y=0; y<Y; y++) { if(*b!=0 && *b!=FlagX) if(TravX(0, y, z, b)) return true; b += yStep; } } return false; }

Ответ отправил: Evgenijm (10-й класс)
Ответ отправлен: 17.12.2010, 04:35
Номер ответа: 264832

Вам помог ответ? Пожалуйста, поблагодарите эксперта за это! Как сказать этому эксперту "спасибо"?

Отправить SMS #thank 264832 на номер 1151 (Россия) | Еще номера »

Отправить WebMoney: Неплохо (30 руб.) Хорошо (40 руб.) Отлично (55 руб.) БЛЕСТЯЩЕ! (80 руб.)

Здравствуйте, уважаемые эксперты! Помогите пожалуйста решить такую задачу: Написать программу нахождения в графе мостов и точек сочленения.

Отправлен: 13.12.2010, 18:47
Вопрос задал: Николаев Александр Валерьевич (Посетитель)
Всего ответов: 1
Страница вопроса »

Приложение:
/* * Поиск мостов и двусвязных компонент. * Даниил Швед, 2008. МФТИ. * danshved [no-spam] gmail.com */ #include <stdio.h> #include <vector> #include <algorithm> #include <conio.h> using namespace std; typedef vector<int> VInt; typedef vector<VInt> VVInt; typedef VInt::iterator VIter; typedef pair<int, int> PInt; typedef vector<PInt> VPInt; typedef vector<VPInt> VVPInt; typedef VPInt::iterator VPIter; VVInt graph; VInt colors, parents, enter, leave, low, bcc; int myTime = 0; int newIndex = 0; /* * Поиск в глубину, выполняющий вычисление enter, leave и low */ void visitLow(int u) { colors[u] = 1; low[u] = enter[u] = ++myTime; printf("visitLow(%d)\n",u); for(VIter it = graph[u].begin(); it != graph[u].end(); it++)//проходим по всем вершинам, связанным с текущей if(colors[*it] == 0) { parents[*it] = u; printf("parents[%d]=%d\n",*it,u); visitLow(*it); low[u] = min(low[u], low[*it]); printf("low[%d]=%d\n",u,low[u]); } else if(colors[*it] == 1 && *it != parents[u]) { low[u] = min(low[u], enter[*it]); printf("low[%d]=%d\n",u,low[u]); } colors[u] = 2; leave[u] = ++myTime; printf("enter[%d]=%d\n",u,enter[u]); } /* * Второй поиск в глубину, присвающий идентификаторы bcc */ void visitBCC(int u) { for(VIter it = graph[u].begin(); it != graph[u].end(); it++) if(parents[*it] == u) { bcc[*it] = (low[*it] < enter[u]) ? bcc[u] : // Дочернее ребро эквивалентно текущему (low[*it] > enter[u]) ? -1 : // Дочернее ребро есть мост (newIndex++); // Дочернее ребро лежит в новой компоненте visitBCC(*it); } } /* * Получение номера BCC, которому принадлежит ребро, или -1, если это мост * (u, v) действительно должно быть ребром, иначе возвращенный результат не имеет смысла */ int getBCC(int u, int v) { return bcc[(enter[u] > enter[v]) ? u : v]; } /* * На входе: числа n и m, затем m описаний ребер * На выходе: список мостов, затем списки ребер в каждой компоненте */ int main() { int n, m, i; FILE *input= fopen("input.txt", "r+"); // Прочитаем граф fscanf(input,"%d%d", &n, &m); graph.resize(n); while(m--) { int from, to; fscanf(input,"%d%d", &from, &to); graph[from].push_back(to); graph[to].push_back(from); } // Запустим первый поиск (вычислим enter и low) colors.assign(n, 0); parents.assign(n, -1); enter.resize(n); leave.resize(n); low.resize(n); for(i = 0; i < n; i++) if(colors[i] == 0) visitLow(i); // Запустим второй поиск (определение идентификаторов bcc) // Второй поиск запускается "по следам" первого, // то есть проходит по уже найденному дереву bcc.assign(n, -1); for(i = 0; i < n; i++) if(parents[i] == -1) visitBCC(i); if(graph[0].size()>1) printf("Tochka sochleneniya %d ee vremya dostupa %d",0,low[i]); for(i = 1; i < n; i++) { //if(low[i]==enter[i]) if(bcc[i]==-1) printf("Tochka sochleneniya %d ee vremya dostupa %d",i,low[i]); } // Выведем результат VPInt bridges; VVPInt comps(newIndex); for(i = 0; i < n; i++) for(VIter it = graph[i].begin(); it != graph[i].end(); it++) if(i < *it) { int id = getBCC(i, *it); ((id == -1) ? bridges : comps[id]).push_back(PInt(i, *it)); } printf("Bridges: "); for(VPIter bridge = bridges.begin(); bridge != bridges.end(); bridge++) printf("(%d, %d) ", bridge->first, bridge->second); printf("\n"); for(i = 0; i < newIndex; i++) { printf("Component %d: ", i); for(VPIter edge = comps[i].begin(); edge != comps[i].end(); edge++) printf("(%d, %d) ", edge->first, edge->second); printf("\n"); } getch(); }

Ответ отправил: Гаряка Асмик (Профессор)
Ответ отправлен: 18.12.2010, 18:40
Номер ответа: 264878

Вам помог ответ? Пожалуйста, поблагодарите эксперта за это! Как сказать этому эксперту "спасибо"?

Отправить SMS #thank 264878 на номер 1151 (Россия) | Еще номера »

Отправить WebMoney: Неплохо (30 руб.) Хорошо (40 руб.) Отлично (55 руб.) БЛЕСТЯЩЕ! (80 руб.)

Здравствуйте, уважаемые эксперты! Можете пожалуйста помочь с разработкой программы ?

Разработать программу для обработки растровых изображений, хранящихся в файлах формата BMP (BitMaP) :

Изменить яркость изображения. Самый простой способ - изменить значения ком-понентов всех пикселей растрового изображения на некоторое значение V, лежа-щее в диапазоне от -255 до 255 (отрицательные значения уменьшают яркость, по-ложительные - увеличивают) и вводимое из текстового файла (но лучше, если из командной строки).

Далее условия выполнения (и есть еще структура заголовка которая приведена в таблице)

(РАСТРОВЫЕ ИЗОБРАЖЕНИЯ

Растровое изображение представляет собой матрицу пикселей, каждый из которых имеет свой цвет. Основными характеристиками растрового изображения являются разре-шение и глубина цвета. Разрешение определяет размеры изображения: его длину (количе-ство пикселей в строке) и высоту (количество строк). Глубина цвета определяет количест-во используемых в изображение цветов и оттенков. Глубина цвета зависит от выбранной цветовой модели, а также количества бит, отводимых для хранения цвета одного пикселя изображения. Наиболее распространенной цветовой моделью является цветовая модель RGB (Red, Green, Blue), в соответствии с которой цвет пикселя изображения описывается последовательностью из трех значений: интенсивность красного, зеленого и синего цве-тов.

ФОРМАТ BMP ФАЙЛА

В файле формата bmp можно выделить несколько частей:
∙ заголовок;
∙ палитра (присутствует только в палитровых изображениях, их рассматривать не будем);
∙ массив пикселей изображения.
Заголовок описывает основные характеристики хранимого изображения. Он располо-жен в самом начале файла и всегда занимает 54 байта. )

Отправлен: 14.12.2010, 01:22
Вопрос задал: Алексей (Посетитель)
Всего ответов: 1
Страница вопроса »

Ответ отправил: Павел Юрьевич (6-й класс)
Ответ отправлен: 14.12.2010, 12:05
Номер ответа: 264762

Оценка ответа: 5

Вам помог ответ? Пожалуйста, поблагодарите эксперта за это! Как сказать этому эксперту "спасибо"?

Отправить SMS #thank 264762 на номер 1151 (Россия) | Еще номера »

Отправить WebMoney: Неплохо (30 руб.) Хорошо (40 руб.) Отлично (55 руб.) БЛЕСТЯЩЕ! (80 руб.)

Здравствуйте, уважаемые эксперты! Прошу Вас помочь с переделыванием структуры программы в приложении, написанной на code blocks, на работу с файлами.
Хочется в готовой программе увидеть следующий вывод: http://s43.radikal.ru/i102/1011/cc/79e92bc876cb.png
Спасибо.

Отправлен: 14.12.2010, 12:07
Вопрос задал: Aptemis (Посетитель)
Всего ответов: 1
Страница вопроса »

Приложение:
#include <stdio.h> #include<string.h> #define SIZE 1//4 #define MAX_BALLS 4 #define INPUT "input.txt" //входной файл #define OUTPUT "output.txt" //выходной файл FILE* f; //переменная файла, если не объявлять её глобально, то надо передавать её как параметр в функцию вывода typedef struct { int year; int month; int day; } Born; typedef struct student { char FIO[20]; char Gender[2]; Born b; char City[30]; int Balls[MAX_BALLS]; } Student; void Print_Group(Student MyMas[], int quant); void Izm_Group(Student MyMas[], Student MyMas1[], int quant, int *quant1); void Sort_Group(Student MyMas[], int quant); int main() { int i; int j; int n1; Student GROUP[SIZE], GROUP1[SIZE]; f = fopen (INPUT, "rt"); for ( i = 0; i<SIZE; i++) { //printf("\nNAME: "); fscanf(f, "%s", GROUP[i].FIO); //printf("Gender(m/f): "); fscanf(f, "%s", GROUP[i].Gender); //printf("City : "); fscanf(f, "%s", GROUP[i].City); //printf("Marks[%d] : ", MAX_BALLS); for (j=0;j<MAX_BALLS;j++) fscanf(f, "%d", &GROUP[i].Balls[j]) ; } fclose (f); f = fopen (OUTPUT, "wt"); fprintf(f, "______________________\n"); fprintf(f, "\nVvednniy spisok[%i]\n",SIZE); fprintf(f, "______________________\n"); Print_Group(GROUP, SIZE); fprintf(f, "\n"); Izm_Group(GROUP,GROUP1,SIZE, &n1); fprintf(f, "______________________\n"); fprintf(f, "Sokrazhennyi spisok[%i]\n",n1); fprintf(f,"______________________\n"); Print_Group(GROUP1,n1); fprintf(f,"\n"); Sort_Group(GROUP1,n1); fprintf(f,"______________________\n"); fprintf(f,"Otsortirovannyi spisok[%i]\n",n1); fprintf(f, "______________________\n"); Print_Group(GROUP1,n1); fclose (f); return 0; } void Print_Group(Student MyMas[], int quant) { int i,j; double sum; fprintf(f,"\Name\tGender\tCity\tMarks\tSredniy ball\n"); for( i=0, sum=0.0; i<quant; i++) { fprintf(f,"%s\t%s\t%s\t", MyMas[i].FIO, MyMas[i].Gender, MyMas[i].City); for(j=0, sum=0.0; j<MAX_BALLS; j++) { sum+=MyMas[i].Balls[j]; fprintf(f,"%d ", MyMas[i].Balls[j]); } fprintf(f, "%4.2f\n", sum/MAX_BALLS); } } void Izm_Group(Student MyMas[], Student MyMas1[], int quant, int *quant1) { int i; for(i=0, *quant1=0; i<quant; i++) if((strcmp(MyMas[i].Gender,"f")==0)&&(strcmp(MyMas[i].City,"SPb")==0)) { MyMas1[(*quant1)++]=MyMas[i]; } } void Sort_Group(Student MyMas[],int quant) { int i, j, k, k1; int sum, sum1; Student buf; for(i=1; i<quant;i++) for(j=quant-1; j>=i;j--) { for (k=0, sum=0; k<MAX_BALLS; k++) sum+=MyMas[j].Balls[k ]; for (k1=0, sum1=0; k1<MAX_BALLS; k1++) sum1+=MyMas[j-1].Balls[k1]; if(sum1>sum) { buf=MyMas[j]; MyMas[j]=MyMas[j-1]; MyMas[j-1]=buf; } } }
-----
Эта история - не для истории, понимаешь?

Ответ отправил: Verena (Профессор)
Ответ отправлен: 18.12.2010, 01:19
Номер ответа: 264863

Вам помог ответ? Пожалуйста, поблагодарите эксперта за это! Как сказать этому эксперту "спасибо"?

Отправить SMS #thank 264863 на номер 1151 (Россия) | Еще номера »

Отправить WebMoney: Неплохо (30 руб.) Хорошо (40 руб.) Отлично (55 руб.) БЛЕСТЯЩЕ! (80 руб.)

Оценить выпуск »
Нам очень важно Ваше мнение об этом выпуске рассылки!

Задать вопрос экспертам этой рассылки »