SGU116
来源:互联网 发布:php学java要多久 编辑:程序博客网 时间:2024/06/07 05:24
题目大意:
在一个由素数组成的数列 2,3,5,7......
定义超级素数,数列中的第k项为超级素数当且仅当k也为素数
给定n
如果n能被超级素数的和所表示
输出最少要多少个超级素数,以及这些超级素数
否则输出impossible
暴力找出超级素数
然后多重背包
//Lib#include<cstdio>#include<cstring>#include<cstdlib>#include<cmath>#include<ctime>#include<iostream>#include<algorithm>#include<vector>#include<string>#include<queue>using namespace std;//Macro#define rep(i,a,b) for(int i=a,tt=b;i<=tt;++i)#define rrep(i,a,b) for(int i=a,tt=b;i>=tt;--i)#define erep(i,e,x) for(int i=x;i;i=e[i].next)#define irep(i,x) for(__typedef(x.begin()) i=x.begin();i!=x.end();i++)#define read() (strtol(ipos,&ipos,10))#define sqr(x) ((x)*(x))#define pb push_back#define PS system("pause");typedef long long ll;typedef pair<int,int> pii;const int oo=~0U>>1;const double inf=1e20;const double eps=1e-6;string name="",in=".in",out=".out";//Varint f[10008],g[10008],p[10008],prime[10008],cnt,pos,n;bool isPrime(int n){if(n<2)return false;rep(i,2,sqrt((double)n))if(n%i==0)return false;return true;}void Pre(){rep(i,2,n)if(isPrime(i))prime[++cnt]=i;int cnt2=cnt;cnt=0;rep(i,1,cnt2){if(isPrime(i))p[++cnt]=prime[i];}}void Work(){scanf("%d",&n);Pre();rep(i,1,n)f[i]=(oo>>1);rep(i,1,cnt)rep(j,p[i],n){if(f[j]>f[j-p[i]]+1){f[j]=f[j-p[i]]+1;g[j]=j-p[i];}}if(f[n]==(oo>>1)){cout<<0<<endl;return;}cout<<f[n]<<endl;pos=n;rep(i,1,f[n]-1)cout<<pos-g[pos]<<' ',pos=g[pos];cout<<pos<<endl;}int main(){//freopen((name+in).c_str(),"r",stdin);//freopen((name+out).c_str(),"w",stdout);//Init();Work();//PS;return 0;}
- SGU116
- sgu116
- sgu116
- SGU116 动态规划 DP
- sgu116:Index of super-prime
- SGU116 Index of super-prime
- SGU 水题集3 SGU116 SGU117 SGU118 SGU119
- oracle中如何移动表空间
- 内建函数解析
- C++程序设计之四书五经/C++学习书籍介绍
- 去掉window关闭时候的弹出框
- SQL Server 2005“备份集中的数据库备份与现有的数据库不同”解决方法
- SGU116
- html5 和 flash的性能比较
- java文件的分割
- UC故事2011/11/17
- 六世达赖喇嘛-仓央嘉措诗集
- web工程部署jsp,servlet
- 接口设计
- 55days
- 在CListView中加入右键菜单!...