[C++] Podział zakresu na k równych obszarów ale integerów

[Borneq]

To coś co ma wiele wspólnego przy powiększaniu rastra metodą Nearest przy niecałkowitej skali powiększenia.

Metoda naiwna:

Kod: Zaznacz cały

void divide_float(int n, int k)
{
    double interval = double(n)/k;
    //double
    for (int i=0; i<k; i++)
    {
        printf("%d %f %f
",i,i*interval, (i+1)*interval);
    }
}

chcę działać na integerach ale i tak mam kłopot z zaokrągleniami:

Kod: Zaznacz cały

void divide_int(int n, int k)
{
    double interval = double(n)/k;
    //int
    uint32_t intPartPos, fracPartPos;
    uint32_t intPartDelta,fracPartDelta,oldFrac;
    intPartPos = 0; fracPartPos = 0;
    intPartDelta = floor(interval);
    fracPartDelta = floor((interval - intPartDelta) * 65536 * 65536);
    for (int i = 0; i<k; i++)
    {
        intPartPos += intPartDelta;
        oldFrac = fracPartPos;
        fracPartPos += fracPartDelta;
        if (oldFrac > fracPartPos) // when LongWord overflow
            intPartPos++;
        printf("%d %d
", i, intPartPos);
    }
}

divide(10,3) zwraca 3,6,9 zamiast 3,6,10 dlatego że w ostatnim jest błąd rzędu \(\displaystyle{ 2^{-32}}\) i to wystarczy by było źle.
Jak wyeliminować liczby zmiennoprzecinkowe i dać dokładnie same inty?

[mihile]

Zamienić

Kod: Zaznacz cały

fracPartDelta = floor((interval - intPartDelta) * 65536 * 65536);

na

Kod: Zaznacz cały

fracPartDelta = ceil((interval - intPartDelta) * 65536 * 65536);

i będzie 3, 6, 10

[Borneq]

Tutaj ceil zadziałało, ale w ogólnym przypadku jest gorzej.
Ale chodzi o rozdzielanie pracy dla wątków - można spokojnie ostatniego elementu nie wyliczać tylko przyjąć ostatni = size().

[Afish]

Kod: Zaznacz cały

#include <iostream>
#include <cstdio>
using namespace std;

void divide(int n, int k){
	int result = n / k;
	int remainder = n % k;
	int bound = k - remainder;
	int current = 0;
	
	for(int i=0;i<k;++i){
		current += i < bound ? result : result + 1;
		printf("%d ", current);
	}
	
	puts("");
}

int main() {
	divide(10, 3);
	divide(10, 4);
	divide(10, 5);
	divide(11, 3);
	divide(11, 4);
	divide(11, 5);
	divide(11, 6);
	return 0;
}