听说你想看题?
1.集合
关于这道题,第一个想到的就是枚举两个数,然后把符合条件的点并到一个集合中去。然而,这算法太low了,不屑一顾(手动滑稽)。
真正的解法
其实只要在使用埃氏筛法的同时,把两个数并在一起就可以了。这样的话复杂度就是并查集复杂度*n
代码
#include<iostream>
#include<vector>
#include<cstring>
#define M 100086
#define FOR(i,a,b) for(register int i=(a),i##_end_=(b);i<=i##_end_;++i)
#define DOR(i,a,b) for(register int i=(a),i##_end_=(b);i>=i##_end_;--i)
using namespace std;
int Fa[M];
int Find(int x){return x==Fa[x]?x:Find(Fa[x]);}
int a,b,p,cnt;
bool mark[M];
int main(){
cin>>a>>b>>p;
FOR(i,1,b)Fa[i]=i;
FOR(i,2,b){
if(mark[i])continue;
int flag=-1;
for(int j=i;j<=b;j+=i){
if(i>=p&&j>=a){
if(flag==-1)flag=Find(j);
else Fa[Find(j)]=flag;
}
mark[j]=1;
}
}
FOR(i,a,b)cnt+=Fa[i]==i;
cout<<cnt<<endl;
return 0;
}
2.能量项链(NOIP2006)
这道题是一个区间DP的题目,关于这个头和尾的处理,我一开始是用结构体处理的。后来啊,发现陈逸神犇的处理非常简单。就是用A[i+n]=A[i]来处理这个问题,这样的话问题就会简单许多。
接下来要做的就是枚举起始点与终点,以及两者之间拿来合并的点,
接下来就是要看子问题的决策了,对于每一个项链,可以认为这个区间的最优方案已经算出来了,然后选取区间中的一个中间点,将两个点接起来,因此可以得到如下递推式:dp[i][j]=max(dp[i][j],dp[i][k]+dp[k+1][j]+A[i]*A[j+1]*A[k+1]) (i表示起点,j表示终点,k表示中间用于连接的点);
代码:
#include<iostream>
#include<cstdio>
#include<cstring>
#define FOR(i,a,b) for(int i=(a),i##_end_=(b);i<=i##_end_;++i)
#define DOR(i,a,b) for(int i=(a),i##_end_=(b);i>=i##_end_;--i)
#define M 205
using namespace std;
int A[M];
int dp[M][M];
int n;
int main() {
cin>>n;
FOR(i,1,n) {
scanf("%d",&A[i]);
A[i+n]=A[i];
}
DOR(i,n*2+1,1)//i±íʾ×ó±ßµã
FOR(j,i+1,n*2+1)//j±íʾÓұߵã
FOR(k,i,j-1)//k±íʾÖмäµã
dp[i][j]=max(dp[i][j],dp[i][k]+dp[k+1][j]+A[i]*A[j+1]*A[k+1]);
int ans=0;
FOR(i,1,n)ans=max(ans,dp[i][n+i-1]);
cout<<ans<<endl;
return 0;
}
总结:
动态规划的题目的很重要的一点就在于子问题的解决与划分,试着去找到子问题与主问题的关系,达到化难为易的效果。