现的从1到N的出现M的次数算法

1. #include <stdio.h>
2. unsigned long func(int a, int n)//求组合c(a,n) 
3. {
4.     unsigned long result = 1, temp = 1;
5.     int i;
6.     for (i = n;i > n-a; i--)
7.     {
8.         result *= i;
9.     }
10.     
11.     for (;a > 0; a--)
12.     {
13.         temp *= a;
14.     }
15.     if (temp != 0)
16.         result /= temp;
17.     return result;
18. }
20. unsigned long fun(int a, int n)//求a^n
21. {
22.     int result;
23.     if (n == 0)
24.     {
25.         result = 1;
26.     }
27.     else
28.     {
29.         result = a;
30.     }
31.     
32.     for (;n > 1;n--)
33.     {
34.         result *= a;
35.     }
36.     return result;
37. }
39. unsigned long _findm(unsigned long a)//a位数，包含m的情况有多少
40. {
41.     unsigned long i, result = 0;
42.     for (i = 1; i <= a; i++)
43.     {//有i个m的情况
44.         result += (func(i, a)*fun(9, a-i));
45.     }
46.     return result;
47. }
49. unsigned long findm(unsigned long arg, unsigned long m)//0～arg之间有多少个含有m的整数
50. {
51.     unsigned long n = 0, a = 1;//arg是n+1位数，第一位数字是a，其余位都是0
52.     unsigned long result;
53.     a = arg;
54.     while (a > 10)
55.     {
56.         a /= 10;
57.         n++;
58.     }
59.     result = _findm(n);
60.     if (a < m)
61.     {
62.         result *= a;
63.     }
64.     else 
65.     {
66.         result *= (a-1);
67.         result += (arg/a);
68.     }
69.     return result;
70. }
72. unsigned long group(unsigned long n, unsigned long a)//将n分组，起点是a，返回终点
73. {
74.     unsigned long result;
75.     int i = 0;
76.     n -= a;
77.     while (n > 10)
78.     {
79.         n /= 10;
80.         i++;
81.     }
82.     result = fun(10, i);
83.     result *= n;
84.     return result;
85. }
87. unsigned long calc(unsigned long n, unsigned long m)//计算从1到N之间的整数有多少个含有数字m
88. {
89.     unsigned long a = 0, b, result = 0;
90.     do
91.     {
92.         b = group(n, a);
93.         result += findm(b, m);
94.         a += b;
95.     } while (b != 0);
96.     return result;
97. }
99. void main()
100. {
101.     unsigned long result, n, m;
102.     //    printf("Input N:");scanf("%d", &n);
103.     //    printf("Input M:");scanf("%d", &m);
104.     n = 2008080808;
105.     m = 7;
106.         result = calc(n, m);
107.     if (m == 0)
108.         result--;
109.     printf("共有%d个整数含有7/n", result);
110. }