rand() бьет в один из интервалов в два раза чаще

[Symba]

Всем привет!
ниже полностью код, под ним описание проблемы.

#include <iostream>
#include <windows.h>
#include <conio.h>
#include <stdio.h>
#include <time.h>
#include <stdlib.h>

using namespace std;
int main()
{
	srand((unsigned)time(NULL));

	int x = 100;

	int t = 8760*x;

	int tvn1 = 33000,
		tvn2 = 29000,
		tvn3 = 34000,
		tvn4 = 35000,
		tvn5 = 19000,
		tvn6 = 7000,
		tvn7 = 7000,
		tvn8 = 3500;

	int tv1=29,
		tv2=22,
		tv3=24,
		tv4=25,
		tv5=37,
		tv6=60,
		tv7=47,
		tv8=135;

	double	notk1 = double(t)/double(tvn1+tv1),
			notk2 = double(t)/double(tvn2+tv2),
			notk3 = double(t)/double(tvn3+tv3),
			notk4 = double(t)/double(tvn4+tv4),
			notk5 = double(t)/double(tvn5+tv5),
			notk6 = double(t)/double(tvn6+tv6),
			notk7 = double(t)/double(tvn7+tv7),
			notk8 = double(t)/double(tvn8+tv8);

	double	votk1 = notk1/t,
			votk2 = notk2/t,
			votk3 = notk3/t,
			votk4 = notk4/t,
			votk5 = notk5/t,
			votk6 = notk6/t,
			votk7 = notk7/t,
			votk8 =	notk8/t;

	double	interval1 = votk1,
			interval2 = interval1+votk2,
			interval3 = interval2+votk3,
			interval4 = interval3+votk4,
			interval5 = interval4+votk5,
			interval6 = interval5+votk6,
			interval7 = interval6+votk7,
			interval8 = interval7+votk8,
			interval9 = 1.00000000;
	double prn;
	
	
	 do{
	 int	     i1 = 0,
				 i2 = 0,
				 i3 = 0,
				 i4 = 0,
				 i5 = 0,
				 i6 = 0,
				 i7 = 0,
				 i8 = 0,
				 i9 = 0;
	  for (int i=1;i<=t;i++){
		  prn = (double)rand() / RAND_MAX;

		  if (prn<=interval1){
			  i=i+tv1;
			  i1=i1+tv1;
		  }else if (prn<=interval2){
			  i=i+tv2;
			  i2=i2+tv2;
		  }else if (prn<=interval3){
			  i=i+tv3;
			  i3=i3+tv3;
		  }else if (prn<=interval4){
			  i=i+tv4;
			  i4=i4+tv4;
		  }else if (prn<=interval5){
			  i=i+tv5;
			  i5=i5+tv5;
		  }else if (prn<=interval6){
			  i=i+tv6;
			  i6=i6+tv6;
		  }else if (prn<=interval7){
			  i=i+tv7;
			  i7=i7+tv7;
		  }else if (prn<=interval8){
			  i=i+tv8;
			  i8=i8+tv8;
		  }else if (prn<=interval9){
			  i9++;
		  }
	  }
	  double	v1 = (double)i1/t,
				v2 = (double)i2/t,
				v3 = (double)i3/t,
				v4 = (double)i4/t,
				v5 = (double)i5/t,
				v6 = (double)i6/t,
				v7 = (double)i7/t,
				v8 = (double)i8/t,
				v9 = (double)i9/t;

				cout << "===============" << endl;
				cout << v1 << " " << (double)i1/x << endl;
				cout << v2 << " " << (double)i2/x << endl;
				cout << v3 << " " << (double)i3/x << endl;
				cout << v4 << " " << (double)i4/x << endl;
				cout << v5 << " " << (double)i5/x << endl;
				cout << v6 << " " << (double)i6/x << endl;
				cout << v7 << " " << (double)i7/x << endl;
				cout << v8 << " " << (double)i8/x << endl;
				cout << v9 << " " << (double)i9/x << endl;
				cout << "===============" << endl;
	  } while (getch()!=0x1B);
}

проблема в том, что ПСЧ попадает во второй интервал чаще в два раза, чем задумано, а в пятый реже в два раза, чем задумано.
Это каким способом решить?
СПАСИБО!
P.S. что это реализуемо через массивы, мне известно, код нужен именно в таком виде =)

[Abstraction]

проблема в том, что ПСЧ попадает во второй интервал чаще в два раза, чем задумано, а в пятый реже в два раза, чем задумано.
Это каким способом решить?

Использовать другой генератор случайных чисел. Относительно этого никто и не обещал, что он будет работать нормально.

[Symba]

Сообщение от Abstraction

Использовать другой генератор случайных чисел. Относительно этого никто и не обещал, что он будет работать нормально.

пример, пожалуйста, или ссылку.

[waleri]

Попробуйте
rand() % range
хотя это тоже не очень правильный способ.
Еще лучше:
rand() * range / RAND_MAX

[Symba]

Сообщение от waleri

Попробуйте
rand() % range
хотя это тоже не очень правильный способ.
Еще лучше:
rand() * range / RAND_MAX

как же мне это реализовать, если мне на промежутке [0,1] ПСЧ нужно?
я так понял range - это интервал?
и по ходу вся проблема в компиляторе...компилю код в VS C++ Express...
только что попробовал в Linux с++ компилировать и QtCreator' ом - все работает как надо....

под Windows 7 и MinGW работает криво....не подскажите компилятор под Win ? =]

[Igor95]

Visual Studio 2012

[Пепел Феникса]

тут не компиль виноват а реализация ГПСЧ.
возьмите стороннюю реализация, например с буста.

[_Bers]

Не гони волну на студийный компилятор, и ищи косяк в своем коде.

Стандартный виндузятный rand() дает равномерное покрытие по всему диапазону (сам лично проверял как работает студийный, когда изобретал свой собственный рандом).

http://codepad.org/HvSIuKpX