To coś co ma wiele wspólnego przy powiększaniu rastra metodą Nearest
przy niecałkowitym powiększeniu.

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

chcę działać na integerach ale i tak mam kłopot z zaokrągleniami:
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\n", i, intPartPos);
}
}

divide(10, 3) zwraca 3,6,9 zamiast 3,6,10 dlatego że w ostatnim jest
błąd rzędu 2^-32 i to wystarczy by było źle.