九度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<stdio.h>
#include<string.h>

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;
}