九度OJ-题目1016 火星A+B

题目描述：
读入两个不超过25位的火星正整数A和B，计算A+B。需要注意的是：在火星上，整数不是单一进制的，第n位的进制就是第n个素数。例如：地球上的10进制数2，在火星上记为“1,0”，因为火星个位数是2进制的；地球上的10进制数38，在火星上记为“1,1,1,0”，因为火星个位数是2进制的，十位数是3进制的，百位数是5进制的，千位数是7进制的……
输入：
测试输入包含若干测试用例，每个测试用例占一行，包含两个火星正整数A和B，火星整数的相邻两位数用逗号分隔，A和B之间有一个空格间隔。当A或B为0时输入结束，相应的结果不要输出。
输出：
对每个测试用例输出1行，即火星表示法的A+B的值。

样例输入：
1,0 2,1
4,2,0 1,2,0
1 10,6,4,2,1
0 0
样例输出：
1,0,1
1,1,1,0
1,0,0,0,0,0

来源：
2006年浙江大学计算机及软件工程研究生机试真题

这道题有点不懂的点是 题目说输入为不超过25位的火星整数，但是代码中定义数组大小为25是不对的，RuntimeError。看了其他博主的文章，他们也有这个疑惑，把数组改大就能AC。

#include
#include

int su[30];

void sushu(){
    int num= 2;
    int i= 0;
    int j;
    while(i< 30){
        for(j= 2; j<= num; j++){
            if(num% j== 0){
                break;
            }
        }
        if(j== num){
            su[i]= num;
            i++;
        }
        num++;
    }
}

int main(){
    char a[100], b[100];
    int ina[300], inb[300];
    int num[300], res[300];
    int i, j, la, lb, lmax;
    memset(a, 0, sizeof(a));
    memset(b, 0, sizeof(b));
    sushu();
    while(scanf("%s %s", a, b)!=EOF){
        memset(ina, 0, sizeof(ina));
        memset(inb, 0, sizeof(inb));
        if(a[0]== '0' && b[0]== '0'){
            break;
        }
        la= strlen(a);
        lb= strlen(b);
        j= 0;
        i= 0;
        while(i< la){
            while(a[i]!=',' && i< la){
                ina[j]= ina[j]*10 + (a[i]-'0');
                i++;
            }
            j++;
            i++;
        }
        la= j;
        j= 0;
        i= 0;
        while(i< lb){
            while(b[i]!=',' && i< lb){
                inb[j]= inb[j]*10 + (b[i]-'0');
                i++;
            }
            i++;
            j++;
        }
        lb= j;
        lmax= la>lb?la:lb;
        int fir[30], sec[30];
        memset(fir, 0, sizeof(fir));
        memset(sec, 0, sizeof(sec));

        for(j= 0; j< lmax- la+1; j++){
            fir[j]=0;
        }
        i= 0;
        for(j= lmax- la+1; j< lmax+1; j++){
            fir[j]= ina[i];
            i++;
        }
        for(j= 0; j< lmax- lb+1; j++){
            sec[j]=0;
        }
        i= 0;
        for(j= lmax- lb+1; j< lmax+1; j++){
            sec[j]= inb[i];
            i++;
        }
        for(i= 0; i< lmax+1; i++){
            num[i]= fir[i]+ sec[i];
        }

        j= 0;
        for(i= lmax; i>= 0; i--){
            res[i]= num[i]%su[j];
            num[i-1]+= num[i]/su[j];
            j++;
        }

        int flag= 0;
        for(j= 0; j< lmax+1; j++){
            if(flag== 0 && res[j]!= 0){
                flag= 1;
                printf("%d,", res[j]);
            }
            else if(flag== 1 && j< lmax){
                printf("%d,", res[j]);
            }
            else if(j== lmax){
                printf("%d\n", res[j]);
            }
        }
        memset(a, 0, sizeof(a));
        memset(b, 0, sizeof(b));
    }
    return 0;
}