大整数相除取余数

写了一个大整数相除取余数的算法，没有在OJ上测试过，有机会测试一下。
题目：给定两个字符串s1和s2，求s1除以s2的余数。

#include <iostream>
#include <string>
#include <cstring>
#include <vector>
#include <stack>
#include <queue>
#include <map>
#include <set>
#include <math.h>
using namespace std;

bool cmp(string s1, string s2)//比较s1是否大于等于s2
{
    int len1 = s1.size(), len2 = s2.size();
    if (len1 > len2)
        return true;
    else if (len1 < len2)
        return false;
    else
    {
        for (int i = 0; i < len1; i++)
            if (s1[i] != s2[i])
                return s1[i]>s2[i];
    }
    return true;
}

string BigSub(string s1, string s2)//大整数相减
{
    int len1 = s1.size(), len2 = s2.size();
    string res; int aux = 0;
    for (int i = 1; i <= len2; i++)
    {
        int d1 = s1[len1 - i] - '0' - aux;
        int d2 = s2[len2 - i] - '0';
        if (d1 < d2)
        {
            d1 += 10;
            aux = 1;
        }   
        else
            aux = 0;
        res += d1 - d2 + '0';
    }
    for (int i = len2+1; i <= len1; i++)
    {
        int d1 = s1[len1 - i] - '0' - aux;
        if (d1 < 0)
        {
            d1 += 10;
            aux = 1;
        }
        else
            aux = 0;
        res += d1 + '0';
    }
    reverse(res.begin(),res.end());
    int pos = res.find_first_not_of("0");
    return pos==string::npos? "0":res.substr(pos);
}
int main()
{
    string s1, s2;
    while (cin >> s1 >> s2)
    {
        while (cmp(s1, s2))//当s1>=s2时做减法
        {
            int len1 = s1.size(), len2 = s2.size();
            if (len1 - len2 > 1)//s1的长度比s2的长度大于等于2时，可以将s2*10的幂次方
                s1 = BigSub(s1, s2.append(string(len1 - len2 - 1, 0)));
            else
                s1 = BigSub(s1, s2);
        }
        cout << s1 << endl;
    }
    return 0;
}