Критерий сериальной корреляции

[ilyakonst]

Сообщение от f.hump

Аватар прав.
и в аккумуляторах делай каст (double)U, (double)V

что это значит? не понимаю

[f.hump]

значит, что U, V должны были остаться целочисленными.

аккумуляторы у тебя c1, c2, c3, z1, z2, z3, z4

ну, и нужно, чтобы при обновлении аккумулятора выполнялась не целочисленная операция а действительная, компилятор конечно может и сам инициативу проявить, но лучше обозначить, что целые должны быть сконвертированы в дабл (прокастить)

c1 += (double)U*(double)V;
.....

фарштанден?

[ilyakonst]

Сообщение от f.hump

значит, что U, V должны были остаться целочисленными.

аккумуляторы у тебя c1, c2, c3, z1, z2, z3, z4

ну, и нужно, чтобы при обновлении аккумулятора выполнялась не целочисленная операция а действительная, компилятор конечно может и сам инициативу проявить, но лучше обозначить, что целые должны быть сконвертированы в дабл (прокастить)

c1 += (double)U*(double)V;
.....

фарштанден?

void korel()
{
int i,n=11;
double C=0, koren;
double c1 = 0, c2 = 0, c3 = 0, z1 = 0, z2 = 0, z3 = 0, z4 = 0, chisl, znam; int U, V;
U = 31;
V = 57;
for (i = 0; i < n; i++) {

U = (40664517 * U + 1) % 4294967296;
//printf("%d\n", U);
V = (40664517 * V + 1) % 4294967296;
//printf("%d\n", V);
c1 += (double)U*(double)V;
c2 += (double)U;
c3 += (double)V;
z1 += (double)U*(double)U;
z2 += (double)U;
z3 += (double)V*(double)V;
z4 += (double)V;

}
//printf("%d\n", z1); printf("%d\n", z2); printf("%d\n", z3); printf("%d\n", z4);
chisl = n*c1 - c2*c3;
znam = (n*z1 - z2*z2)*(n*z3 - z4*z4);
printf("chislitel:%lld\n", chisl);
printf("znamenatek:%lld\n", znam);
koren = sqrt(double(znam));
printf("koren iz znam:%3.lf\n", koren);
C = chisl / koren;
printf("koefficent: %.3f\n",C);

}

вот так получается? вроде работает

[f.hump]

это конечно прогресс.
поздравляю.

но U, V были long long и должны остаться long long, иначе последовательности окажутся не такими как планировалось.

можешь оставить int, но сделай его unsigned int и замени

U = (40664517 * U + 1) % 4294967296;
//printf("%d\n", U);
V = (40664517 * V + 1) % 4294967296;

на

U = (40664517 * U + 1);
//printf("%d\n", U);
V = (40664517 * V + 1);

[ilyakonst]

Почему я должен откинуть нахождение остатка? это же сама суть лкг

[f.hump]

ну, потому что целочисленные операции на unsigned int устроены таким образом, что результат любой операции будет равен остатку от деления на 4294967296 результата полноемкостной операции.