C语言混合类型运算与unsigned long溢出处理问题咨询

我需要实现一个函数，将unsigned long类型的计数器值乘以一个非负32位float因子，返回unsigned long类型结果；当结果溢出unsigned long范围时，返回其模该类型最大值的余数。相关变量定义如下：

• IntVars.Counter_ui32：unsigned long类型，目标平台范围为2^{32，最大值为`4294967295`（2}32-1）
• IntVars.Factor_fl32（及本地变量Result）：32位float类型，最大值约3.4E+38，由其他函数保证非负

目前遇到的问题：

• 无法使用math.h的fmod()函数，尝试用while循环减去最大值的方案时，遇到类型转换和跨类型溢出问题
• 疑问float是否有足够有效位表示unsigned long的值
• 测试代码的平台与目标平台类型最大值不同（测试平台unsigned long最大值为18446744073709551615），已知仅丢失小数部分的舍入误差不影响业务，想确认还有其他需要注意的问题吗？

测试代码 #

/**************************

                            Online C Compiler.
                Code, Compile, Run and Debug C program online.
Write your code in this editor and press "Run" button to compile and execute it.

***************************/

#include <stdio.h>
#include <float.h>

//#define UI32MAX          ((unsigned long)4294967295UL)
#define UI32MAX          ((unsigned long)18446744073709551615ul)

typedef struct 
{
    unsigned long counter;
    float factor;
    
}IntVars_ts;

unsigned long Function(IntVars_ts IntVars);

int main()
{
    unsigned long start = UI32MAX - 3ul;
    printf("Maximum Unsigned Long %lu\n",(unsigned long)~0);
    printf("Maximum float %f\n\n",FLT_MAX);
    
    IntVars_ts IntVars ={start,1.0f}; /* Constant factor, I only increment the counter\
                                      \in this example */
    
    while(1)
    {
            printf("The unsigned long counter value is %lu\n",IntVars.counter);
            printf("The float counter value is %f\n",(float)IntVars.counter);
            printf("The result is %lu\n\n",Function(IntVars));
            IntVars.counter++;
    }
    


    return 0;
}

unsigned long Function(IntVars_ts IntVars)
{
    float Result = (float)IntVars.counter*IntVars.factor;
    
    while(Result >= (float)UI32MAX)
    {
        printf("Result is %f\n",Result);
        printf("Max is %f\n\n",(float)UI32MAX);
        Result -= (float)UI32MAX;
    }
    
    return (unsigned long)Result;
}

代码输出 #

/tmp/z2jTCSG5B0.o
Maximum Unsigned Long 18446744073709551615
Maximum float 340282346638528859811704183484516925440.000000

The unsigned long counter value is 18446744073709551612
The float counter value is 18446744073709551616.000000
Result is 18446744073709551616.000000
Max is 18446744073709551616.000000

The result is 0

The unsigned long counter value is 18446744073709551613
The float counter value is 18446744073709551616.000000
Result is 18446744073709551616.000000
Max is 18446744073709551616.000000

The result is 0

The unsigned long counter value is 18446744073709551614
The float counter value is 18446744073709551616.000000
Result is 18446744073709551616.000000
Max is 18446744073709551616.000000

The result is 0

The unsigned long counter value is 18446744073709551615
The float counter value is 18446744073709551616.000000
Result is 18446744073709551616.000000
Max is 18446744073709551616.000000

The result is 0

The unsigned long counter value is 0
The float counter value is 0.000000
The result is 0

The unsigned long counter value is 1
The float counter value is 1.000000
The result is 1

The unsigned long counter value is 2
The float counter value is 2.000000
The result is 2

The unsigned long counter value is 3
The float counter value is 3.000000
The result is 3

1. float的精度限制问题 #

32位float只有24位有效二进制位（约7-8位十进制有效数字）：

• 对于目标平台的32位unsigned long（最大值4294967295，10位十进制）：当数值超过2^{24（16777216）后，float无法精确表示所有整数，相邻可表示的float值间隔会大于1。比如2}24到2^25之间的数值，float的最小间隔是2，所有奇数都会被舍入到相邻偶数。
• 测试代码中64位unsigned long的最大值远超float精确范围，转换为float后直接被舍入到18446744073709551616.0，这是测试结果错误的核心原因。

2. while循环的效率与溢出风险 #

• 效率问题：如果Result远大于UI32MAX，循环会执行数百甚至数千次，严重拖慢性能。
• 溢出风险：虽然float最大值很大，但极端情况下多次减法操作仍可能导致数值超出float范围，触发未定义行为。

3. 类型转换的隐性问题 #

• unsigned long转float时，若数值超出精确范围会发生舍入，丢失原始整数信息，后续模运算结果必然不准确。
• float转unsigned long时，若float值大于unsigned long最大值，C标准定义为未定义行为——即使做了减法，也必须确保最终Result在0到UI32MAX之间。

4. 高效替代实现（无需fmod） #

通过数学推导避免循环，利用浮点运算快速计算余数：

unsigned long Function(IntVars_ts IntVars) {
    const unsigned long UI32MAX = 4294967295UL;
    const float UI32MAX_FLT = (float)UI32MAX;
    
    float product = (float)IntVars.counter * IntVars.factor;
    
    // 处理乘积未溢出的情况
    if (product <= UI32MAX_FLT) {
        return (unsigned long)product;
    }
    
    // 处理无穷大的极端情况
    if (product == INFINITY) {
        return 0; // 或根据业务需求调整
    }
    
    // 计算近似商并取整
    float quotient = product / UI32MAX_FLT;
    unsigned long int_quotient = (unsigned long)quotient;
    
    // 计算余数
    float remainder = product - (float)int_quotient * UI32MAX_FLT;
    
    // 修正浮点精度误差导致的余数超限
    if (remainder > UI32MAX_FLT) {
        remainder -= UI32MAX_FLT;
    }
    
    return (unsigned long)remainder;
}

该方案仅需几次浮点运算，效率远高于循环，同时规避了多次减法的溢出风险。

额外注意点 #

• 目标平台是32位unsigned long，需确保代码中UI32MAX定义正确，避免使用测试平台的64位值。
• 若因子Factor_fl32极大，乘积可能变为INFINITY，需提前判断并处理。

内容的提问来源于stack exchange，提问作者Hi.