史密斯数(C++实现)

若一个合数的质因数分解式逐为位相加之和等于其本身逐位相加之和,则称这个数为smith数。如4937775=3×5×5×65837,而3+5+5+6+5+8+3+7=42,4+9+3+7+7+7+5=42,所以4937775是smith数。现给定两个数L和R,求L和R之间的史密斯数。

算法分析:
从L~R的正整数开始遍历寻找smith数,若这个数是素数则舍弃,若不是素数,则分解它的质因数再逐位相加求和,如果和它的数字和相等,那么输出。

废话不多说,上代码

#define uns using namespace std;
uns
int sz(int n){
	int ans=0;
	while(n){
		ans+=n%10;
		n/=10;
	}
	return ans;
}
int zyz(int n){
	int ans=0;
	while(n>1){
		for(int i=2;i<=n;i++){
			if(n%i==0){
				ans+=sz(i);
				n/=i;
				break;
			}
		}
	}
	return ans;
}
int pss(int n){
	if(n<2){
		return 0;
	}
	for(int i=2;i<=sqrt(n);i++){
		if(n%i==0){
			return 0;
		}
	}
	return 1;
}
main(){
	int l,r,g=0;
	cin>>l>>r;
	for(int i=max(l,1);i<r+1;i++){
		if(sz(i)==zyz(i)&&!pss(i)){
			cout<<i<<"\n";
		}
	}
}
  • 2
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Smith-Waterman算法是一种用于序列比对的动态规划算法。它可以用于比对DNA、RNA、蛋白质序列等。C++是一种高效的编程语言,可以用于实现Smith-Waterman算法实现Smith-Waterman算法C++代码需要考虑以下几个方面: 1. 输入序列:需要从文件或者用户输入中读取待比对的序列。 2. 动态规划矩阵:需要创建一个二维组来存储动态规划矩阵。 3. 算法实现:需要实现算法的核心部分,包括初始化矩阵、计算得分、回溯路径等。 4. 输出结果:需要将比对结果输出到文件或者屏幕上。 以下是一个简单的Smith-Waterman算法C++实现示例: ```c++ #include <iostream> #include <fstream> #include <cstring> using namespace std; const int MAXLEN = 1000; const int GAP = -2; const int MATCH = 3; const int MISMATCH = -1; int score[MAXLEN][MAXLEN]; int max(int a, int b, int c) { int m = a; if (b > m) m = b; if (c > m) m = c; return m; } void init(int n, int m) { for (int i = 0; i <= n; i++) { score[i][0] = 0; } for (int j = 0; j <= m; j++) { score[0][j] = 0; } } void compute_score(char* s1, char* s2, int n, int m) { for (int i = 1; i <= n; i++) { for (int j = 1; j <= m; j++) { int match = score[i-1][j-1] + (s1[i-1] == s2[j-1] ? MATCH : MISMATCH); int delete_gap = score[i-1][j] + GAP; int insert_gap = score[i][j-1] + GAP; score[i][j] = max(match, delete_gap, insert_gap); } } } void traceback(char* s1, char* s2, int n, int m) { int i = n, j = m; while (i > 0 && j > 0) { int match = score[i-1][j-1] + (s1[i-1] == s2[j-1] ? MATCH : MISMATCH); int delete_gap = score[i-1][j] + GAP; int insert_gap = score[i][j-1] + GAP; if (score[i][j] == match) { cout << s1[i-1] << " " << s2[j-1] << endl; i--; j--; } else if (score[i][j] == delete_gap) { cout << s1[i-1] << " -" << endl; i--; } else { cout << "- " << s2[j-1] << endl; j--; } } } int main() { char s1[MAXLEN], s2[MAXLEN]; cin >> s1 >> s2; int n = strlen(s1), m = strlen(s2); init(n, m); compute_score(s1, s2, n, m); traceback(s1, s2, n, m); return 0; } ``` 该示例代码实现Smith-Waterman算法的核心部分,包括初始化矩阵、计算得分、回溯路径等。在输入两个待比对的序列后,程序会输出比对结果。需要注意的是,该示例代码只是一个简单的实现,实际应用中可能需要进行更多的优化和改进。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值