MyException - 我的异常网
当前位置:我的异常网» 编程 » 几个分组双肩包题 zoj 3450 hdu 4341 hdu 4345

几个分组双肩包题 zoj 3450 hdu 4341 hdu 4345

www.MyException.Cn  网友分享于:2013-09-07  浏览:18次
几个分组背包题 zoj 3450 hdu 4341 hdu 4345

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3450

http://acm.hdu.edu.cn/showproblem.php?pid=4345

http://acm.hdu.edu.cn/showproblem.php?pid=4341


 zoj 3450 和 hdu 4341 是同样的题,- -!,真不知道多校审题是怎么审的

可以直接极角排序,对同一条线上的点处理一下,

即要么区第一个,要么取第一二两个。。。。

一条线就是一组物品,每组物品最多选一个,典型的分组背包,如果用一维数组的话要小心物品的体积为0的情况。


zoj 3450 

这题极角排序的话好像中间计算结果会溢出,用double可能会有浮点误差,所以直接用数学方法做了

#include<cstdio>
#include<cstring>
#include<vector>
#include<cstdlib>
#include<map>
#include<cmath>
#include<algorithm>
using namespace std;
const int N_MAX = 510;
int dp[10010];
int gcd(int a,int b){
     return !b ? a : gcd(b,a%b);
}
int main()
{
	int T,n;
	int x,y,z,t,val,x0,y0;
	while(scanf("%d%d%d%d",&x0,&y0,&n,&T)!=EOF)
	{
	    map<pair<int,int>,vector<pair<int,pair<int,int> > > > mp;
		for(int i=0;i<n;i++)
		{
			scanf("%d%d%d%d",&x,&y,&t,&val);
			x-=x0;y-=y0;
			int z=gcd(abs(x),abs(y));
			x/=z;y/=z;
			mp[make_pair(x,y)].push_back(make_pair(z,make_pair(t,val)));
		}
		vector<pair<int,int> > edge[N_MAX];
		int tot=0;
		for( map<pair<int,int>,vector<pair<int,pair<int,int> > > >::iterator it=mp.begin();it!=mp.end();it++)
		{
			sort(it->second.begin(),it->second.end());
			int sum1=0,sum2=0;
			for(vector<pair<int,pair<int,int> > >::iterator i = it->second.begin();i!=it->second.end();i++)
			{
				sum1+=i->second.first;
				sum2+=i->second.second;
				edge[tot].push_back(make_pair(sum1,sum2));
			}
			tot++;
		}
        fill(dp,dp+T+1,0);
		for(int i=0;i<tot;i++)
		{
			for(int j=T;j>=0;j--)
			{
				int sz=edge[i].size();
				int tmp=dp[j];//体积可能会等于0,所以要用一个临时变量记录最优值
				for(int k=0;k<sz;k++)
				{
					if(j>=edge[i][k].first)
					{
						tmp=max(tmp,dp[j-edge[i][k].first]+edge[i][k].second);//否则会导致 dp[j]=max(dp[j],dp[j-0]+second),自己更新自己
					}
				}
				dp[j]=tmp;
			}
		}
		printf("%d\n",dp[T]);
	}
	return 0;
}

hdu 4341

#include<cstdio>
#include<cstring>
#include<vector>
#include<algorithm>
using namespace std;
struct Point{
    int x,y,t,val;
    Point(int _x,int _y):x(_x),y(_y){}
    Point(){}
}in[210];
int dp[40010];
Point cent;
inline int det(Point a,Point b,Point c){
    return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
}
int dis(Point a,Point b){
    return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}
int cmp(Point a,Point b){
    if( det(cent,a,b)==0) return dis(a,cent)<dis(b,cent);
    return det(cent,a,b)>0;
}
int main(){
    int n,T,ca=1;
    while(scanf("%d%d",&n,&T)!=EOF){
        vector<pair<int,int> > edge[210];
        for(int i=0;i<n;i++)
            scanf("%d%d%d%d",&in[i].x,&in[i].y,&in[i].t,&in[i].val);
        cent=Point(0,0);
        sort(in,in+n,cmp);
        int tot=0;
        for(int i = 0, j ; i < n ; i = j){
            int sum1=in[i].t,sum2=in[i].val;
            edge[tot].push_back(make_pair(sum1,sum2));
            for(j = i+1; det(cent,in[j],in[i]) == 0 && j < n; j++){
                sum1 += in[j].t,sum2 += in[j].val;
                edge[tot].push_back(make_pair(sum1,sum2));
            }
            tot++;
        }
        memset(dp,0,sizeof(dp));
        for(int i=0;i<tot;i++)
            for(int v=T;v>=0;v--)
                for(int j=0;j<edge[i].size();j++)
                    if(v>=edge[i][j].first)
                       dp[v]=max(dp[v],dp[v-edge[i][j].first]+edge[i][j].second);
        printf("Case %d: %d\n",ca++,dp[T]);
    }
    return 0;
}


hdu 4345

将一个数分成若干份,求这若干份数的不同的最小公倍数的总数

注意 到每个数最终都可以分解为若干质数比如4 * 5 就相当于2^2 * 5,所以就可以分组了

2  2^2  2^3  2^4  2^5 .....
3  3^2  3^3 3^4   3^5.......

.....

每一组最多选一个物品,最小公倍数肯定不会重复(因为都是质数的不同组合)

dp[i]表示各个数的和小于等于i时 总的最小公倍数的个数,因为dp[0]=1,用0的容量装体积为0,所以dp[1~V]=1

最后的答案是dp[V]

而如果状态定义为各个数的和恰好为i时的总最小公倍数的个数dp[0]=1,dp[1~V]=0;

最后的答案是sigma(dp[i])

#include<cstdio>
#include<cstring>
#include<vector>
#include<cmath>
#include<algorithm>
using  namespace std;
__int64 dp[1010];
vector<int> edge[200];
bool isp(int num){
    for(int i=2;i<=sqrt(num*1.0);i++){
        if(num%i==0) return false;
    }
    return true;
}
int main(){
    int n;
    while(scanf("%d",&n)!=EOF){
        int tot=0;
        for(int i=0;i<200;i++) edge[i].clear();
        for(int i=2;i<1000;i++){
            if(isp(i)){
                for(int j=i;j<=1000;j*=i){
                     edge[tot].push_back(j);
                }
                tot++;
            }
        }
        for(int i=0;i<=n;i++) dp[i]=1;
        for(int i=0;i<tot;i++)
            for(int v=n;v>=0;v--)
                for(int j=0;j<edge[i].size()&&edge[i][j]<=v;j++)
                        dp[v]+=dp[v-edge[i][j]];
        printf("%I64d\n",dp[n]);
    }
    return 0;
}


文章评论

要嫁就嫁程序猿—钱多话少死的早
要嫁就嫁程序猿—钱多话少死的早
如何区分一个程序员是“老手“还是“新手“?
如何区分一个程序员是“老手“还是“新手“?
程序员眼里IE浏览器是什么样的
程序员眼里IE浏览器是什么样的
10个调试和排错的小建议
10个调试和排错的小建议
十大编程算法助程序员走上高手之路
十大编程算法助程序员走上高手之路
亲爱的项目经理,我恨你
亲爱的项目经理,我恨你
科技史上最臭名昭著的13大罪犯
科技史上最臭名昭著的13大罪犯
看13位CEO、创始人和高管如何提高工作效率
看13位CEO、创始人和高管如何提高工作效率
团队中“技术大拿”并非越多越好
团队中“技术大拿”并非越多越好
老程序员的下场
老程序员的下场
Java程序员必看电影
Java程序员必看电影
“懒”出效率是程序员的美德
“懒”出效率是程序员的美德
鲜为人知的编程真相
鲜为人知的编程真相
不懂技术不要对懂技术的人说这很容易实现
不懂技术不要对懂技术的人说这很容易实现
“肮脏的”IT工作排行榜
“肮脏的”IT工作排行榜
Java 与 .NET 的平台发展之争
Java 与 .NET 的平台发展之争
如何成为一名黑客
如何成为一名黑客
程序员应该关注的一些事儿
程序员应该关注的一些事儿
程序员必看的十大电影
程序员必看的十大电影
我跳槽是因为他们的显示器更大
我跳槽是因为他们的显示器更大
代码女神横空出世
代码女神横空出世
编程语言是女人
编程语言是女人
中美印日四国程序员比较
中美印日四国程序员比较
60个开发者不容错过的免费资源库
60个开发者不容错过的免费资源库
程序员的鄙视链
程序员的鄙视链
Web开发者需具备的8个好习惯
Web开发者需具备的8个好习惯
当下全球最炙手可热的八位少年创业者
当下全球最炙手可热的八位少年创业者
旅行,写作,编程
旅行,写作,编程
每天工作4小时的程序员
每天工作4小时的程序员
为什么程序员都是夜猫子
为什么程序员都是夜猫子
5款最佳正则表达式编辑调试器
5款最佳正则表达式编辑调试器
什么才是优秀的用户界面设计
什么才是优秀的用户界面设计
做程序猿的老婆应该注意的一些事情
做程序猿的老婆应该注意的一些事情
程序员都该阅读的书
程序员都该阅读的书
10个帮程序员减压放松的网站
10个帮程序员减压放松的网站
程序员最害怕的5件事 你中招了吗?
程序员最害怕的5件事 你中招了吗?
我的丈夫是个程序员
我的丈夫是个程序员
写给自己也写给你 自己到底该何去何从
写给自己也写给你 自己到底该何去何从
老美怎么看待阿里赴美上市
老美怎么看待阿里赴美上市
程序猿的崛起——Growth Hacker
程序猿的崛起——Growth Hacker
我是如何打败拖延症的
我是如何打败拖延症的
漫画:程序员的工作
漫画:程序员的工作
 程序员的样子
程序员的样子
程序员和编码员之间的区别
程序员和编码员之间的区别
那些争议最大的编程观点
那些争议最大的编程观点
软件开发程序错误异常ExceptionCopyright © 2009-2015 MyException 版权所有