不清楚有没有,用scilab算很快,大体流程是
1．判断奇偶（2除外）,偶数放弃
2．奇数加1,取加1后的一半,然后判断是否整除
3．用循环模式
即可