题意
分析
- 首先一眼DP好吧,那么考虑考虑怎么设计DP呢?
- 那么我们按题意设计dp数组
真不知道怎么引领你们想出这个dp数组
解法1
- 设dp[i] 为从1号到i号石头可以组成的合法方案数
- 那么转移方程就是
d p [ i ] = ∑ j = 1 i d p [ j ] dp[i]=\sum_{j=1}^{i}dp[j] %从m到n求和 dp[i]=∑j=1idp[j]
(当 ∑ j + 1 i \sum_{j+1}^{i} %从m到n求和 ∑j+1i是质数)
解法2 (时间不大对,O(n 3))
-
设dp[i] 为从1号到i号石头可以组成的合法方案数
-
设dp[i][j] 为从i号到j号石头可以组成的合法方案数
-
那么我们想象一下现在正在处理dp[i][j],那么暴力枚举中间节点k,做一个Floyd 式dp即可。
Floyd算是个dp,这个应该知道吧 -
那么转移方程就是
d p [ i ] [ j ] = ∑ k = i j − 1 d p [ k + 1 ] [ j ] dp[i][j]=\sum_{k=i}^{j-1}dp[k+1][j] %从m到n求和 dp[i][j]=∑k=ij−1dp[k+1][j]
(其中 ∑ i k \sum_{i}^{k} %从m到n求和 ∑ik是质数) -
值得注意的是,如果i,j整个区间不处理就是一种方案,那么我们是没法累加上这个的,所以需要特判一下。
图示
废话
那么我为什么写上面的假算法呢,因为我的同学竟然把它卡过去了。他卡的方法是把long long 换成 int 我傻了。这个我是现场想了的,觉得不大行,加上解法1的小优化,想了想,应该是过不去,所以不调了。换写法了。我竟然因为没判断1的情况而少了20分。/呜呜呜
代码
解法1
#include <bits/stdc++.h>
using namespace std;
const int N = 1010;
const int M = 987654321;
int pre[6000000], vis[6000010];
int n;
int flag[N][N];
int add[N];
int f[N];
int dis[N][N] = {};
int ln[N], tot;
int ans, t[N];
long long dp[N][N];
inline int read()
{
int Num=0,f=1; char ch=getchar();
while(ch<'0'||ch>'9') {if(ch=='-') f=-1; ch=getchar();}
while(ch>='0'&&ch<='9') {Num=(Num<<1)+(Num<<3)+ch-'0'; ch=getchar();}
return Num*f;
}
void prepare()
{
vis[1]=1;
for(int i=2;i<=5000000;i++)
{
if(!vis[i]) pre[++pre[0]]=i;
for(int j=1;j<=pre[0]&&i*pre[j]<=5000000;j++)
{
vis[i*pre[j]]=1;
if(i%pre[j]==0) break;
}
}
}
int main() {
prepare();
n = read();
for (int i = 1; i <= n; i++) {
f[i] = read();
if (!vis[f[i]])
dis[i][++dis[i][0]] = i, dp[i][i] = 1;
add[i] = add[i - 1] + f[i];
}
for (int i = n; i >= 1; i--)
for (int j = i - 1; j >= 1; j--)
if (!vis[add[i] - add[j - 1]])
dis[i][++dis[i][0]] = j, dp[j][i] = 1;
for (int i = 2; i <= n; i++) {
for (int j = 1; j <= dis[i][0]; j++) dp[1][i] += dp[1][dis[i][j] - 1] % M;
}
printf("%lld\n", dp[1][n] % M);
return 0;
}
// f[i][j]从i号到j号的方案数
// 1.枚举左右端点组合
// 2.i枚举之前位点,然后相加,转移
解法2(不是我写的)
#include<bits/stdc++.h>
using namespace std;
const int mod=987654321;
inline int qr()
{
int x=0,f=1;
char ch=getchar();
while(ch>'9'||ch<'0')
{
if(ch=='-') f=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9')
{
x=x*10+ch-'0';
ch=getchar();
}
return x*f;
}
int o[1000010];
int prime[800000];
int cnt;
void zhi(int p)
{
for(int i=2;i<=p;i++)
{
if(!o[i])
{
prime[++cnt]=i;
}
for(int j=1;i*prime[j]<=p;j++)
{
o[i*prime[j]]=1;
if(i%prime[j]==0) break;
}
}
}
int n;
int f[1010][1010];
int a[1010],b[1010];
long long m;
int main()
{
zhi(1000000);
n=qr();
o[1]=1;
for(register int i=1;i<=n;i++)
{
a[i]=qr();
b[i]=b[i-1]+a[i];//前缀和
if(o[a[i]]==0)
f[i][i]=1;
}
for(register int i=2;i<=n;i++)
{
for(register int l=1;l+i-1<=n;l++)
{
int r=l+i-1;
for(register int k=l;k<r;k++)
{
if(o[b[k]-b[l-1]]==0)
{
m=(f[l][r]+f[k+1][r])%mod;
f[l][r]=m;
}
}
if(o[b[r]-b[l-1]]==0)
{
f[l][r]=f[l][r]+1;
}
}
}
printf("%d",f[1][n]%mod);
return 0;
}