求1到10000之间的所有满足个位数字的立方和等于它本身的数（即水仙花数）的代码

我知道怎么写啦
#include<stdio.h>
void main()
{
  int a,b,c,sum;
  for(a=1;a<=10000;a++)
  {
    sum=0;
    c=a;
    do{
       b=c%10;
       c=c/10;
       sum=sum+b*b*b;
       }while(c!=0);
       if(sum==a)
       printf("%d\n",a);
    }
}