石头

题意

在这里插入图片描述

分析

  • 首先一眼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=ij1dp[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;
 }

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值