NOIP2007 矩阵取数 DP+高精度

最新推荐文章于 2019-06-01 16:58:43 发布

ztzshiwo_cnyali

最新推荐文章于 2019-06-01 16:58:43 发布

阅读量278

点赞数

文章标签： NOIP2007 高精度 dp

本文链接：https://blog.csdn.net/ztzshiwo_cnyali/article/details/52185947

版权

这个题并不难想……但是高精度和DP的组合有点恶心……

不说了都是QAQ

代码如下：

我用的重定义运算符。。。乘法小小的偷个懒。。。因为只有*2...

/* 
手动测NOIP数据全过
但是评测30分。。。
不做了。。。
另:A了
是比较运算符的问题 用刘汝佳的过了。。。 
*/ 
#include<cmath>
#include<cstdio>
#include<vector>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
struct bigint
{
	vector<long long>s;
	void clean()
	{
		s.clear();
		s.push_back(0);
	}
	bigint operator + (const bigint& b)const
	{
		int i;
		bigint c;
		c.s.clear();
		int maxs=s.size(),maxb=b.s.size();
		for(int i=0,g=0;;i++)
		{
			if(g==0&&i>=maxs&&i>=maxb)break;
			int x=g;
			if(i<maxs)x+=s[i];
			if(i<maxb)x+=b.s[i];
			c.s.push_back(x%10000);
			g=x/10000;
		}
		return c;
	}
	bool operator < (const bigint& b)const
	{
		if(s.size()!=b.s.size())return s.size()<b.s.size();
		for(int i=s.size()-1;i>=0;i--)
			if(s[i]!=b.s[i])return s[i]<b.s[i];
		return false;
	}
}ans,f[200][200],mp[200];
bigint max1(bigint a,bigint b)
{
	if(a<b)return b;
	return a;
}
bigint mult(bigint a)
{
	return a+a;
}	
void fout(long long x)
{
	if(x<10)printf("000%lld",x);
	else
	if(x<100)printf("00%lld",x);
	else
	if(x<1000)printf("0%lld",x);
	else
	printf("%lld",x);
	return;
}
void putout(bigint a)
{
	int flag=1;
	for(int i=a.s.size()-1;i>=0;i--)
		if(a.s[i]==0&&flag==1)continue;
		else
		{
			if(flag==1)printf("%lld",a.s[i]);
			else fout(a.s[i]);
			flag=0;
		}
	if(flag==1)printf("0\n"); 
	return;
}
int n,m;
int main()
{
	int x;
	ans.clean();
	scanf("%d%d",&n,&m);
	for(int i=0;i<=100;i++)
		for(int j=0;j<=100;j++)
			f[i][j].clean();
	for(int i=1;i<=n;i++)
	{
		bigint sum;
		sum.clean();
		for(int j=1;j<=m;j++)
		{
			mp[j].s.clear();
			scanf("%d",&x);
			mp[j].s.push_back(x);
		}
		for(int x=m;x>=1;x--)
			for(int y=x;y<=m;y++)
				f[x][y]=mult(max1(f[x+1][y]+mp[x],f[x][y-1]+mp[y]));
		ans=ans+f[1][m];
	}
	putout(ans);
	return 0;
}