问题描述:
某人最喜欢的数字是1!所以他经常会使用一些手段,把一些非1的数字变 成1,并为此得意不已。他会且仅会的两种手段是:
1.把某个数m除以某个质数p——当然p必须能整除这个数,即m=m/p
2.把某个数m减1,即m=m-1
有一天他突发奇想,想把[a,b]区间中所有的数一个一个地变成1,这是一个巨大的无聊的工程,所以他想知道他最少得花多少操作才能达到目 的。
Input
输入包含多组数据(1000组数据),EOF结束。
每组数据以两个整数开头:a,b(0<a<=b<=100000),意义如题意描述。
Output
每组数据输出一行,最少操作数。
Sample Input
2 3
3 5
11 12
Sample Output
2
4
3
算法:
1,这个题目的意思是求某一个数具体操作的最小值,也就是从1到1000000这1000000个数,每个数都对应一个值,最后按区间相加即可。
2,有两个操作,+1和乘以质数。所以得到一个数的值,(我定义了一个数组c[1000000]来存放这个最小值)c[i]可能有多种情况得到,可能是c[i-1]+1或者是c[j]+1其中j成一个质数等于i,所以c[i]只需要取这些值中最小的了。这样就得到了最优子结构,明显是动态规划DP
C[i]=MIN{C[i-1]+1,c[j]+1}其中j乘以一个质数为i
3,所以我的程序用动态规划的思想,自底向上来构造c[i],也就是从c[1]开始
首先,求出1到1000000所有的质数存放到一个vector里面
然后用一个队列开始动态规划具体代码如下
#include<iostream>
#include<vector>
#include<queue>
using namespace std;
#define num 100001
int main()
{
vector<unsigned int> zhi;
unsigned int zhishu[num]={0};
unsigned int c[num]={0};
queue<unsigned> qu;
for(unsigned int i=2;i<num/2;i++)
{
if(zhishu[i]==0)
{
for(unsigned int j=2*i;j<num;j+=i)
{
zhishu[j]=1;
}
}//end if
}
//把所有质数存到vector里面
for(unsigned int i=2;i<num;i++)
{
if(zhishu[i]==0)
zhi.push_back(i);
}
qu.push(1);
while(!qu.empty())
{
unsigned int x=qu.front();
qu.pop();
//加一操作
if(x+1<num&&c[x+1]==0)
{
c[x+1]=c[x]+1;
qu.push(x+1);
}
//乘以质数操作
vector<unsigned>::iterator ite;
for(ite=zhi.begin();ite!=zhi.end();ite++)
{
unsigned long y=x*(*ite);
if(y>=num)//这里一定要优化,否则很慢
break;
if(y<num&&c[y]==0)
{
c[y]=c[x]+1;
qu.push(y);
}
}
}
//输入的数,直接可以读取
unsigned int a,b;
while(cin>>a>>b)
{
unsigned int sum=0;
for(unsigned i=a;i<=b;i++)
sum+=c[i];
cout<<sum<<endl;
}
}