Codeforces 8C Looking for Order(状态压缩DP)

Girl Lena likes it when everything is in order, and looks for order everywhere. Once she was getting ready for the University and noticed that the room was in a mess — all the objects from her handbag were thrown about the room. Of course, she wanted to put them back into her handbag. The problem is that the girl cannot carry more than two objects at a time, and cannot move the handbag. Also, if he has taken an object, she cannot put it anywhere except her handbag — her inherent sense of order does not let her do so.

You are given the coordinates of the handbag and the coordinates of the objects in some Сartesian coordinate system. It is known that the girl covers the distance between any two objects in the time equal to the squared length of the segment between the points of the objects. It is also known that initially the coordinates of the girl and the handbag are the same. You are asked to find such an order of actions, that the girl can put all the objects back into her handbag in a minimum time period.

Input

The first line of the input file contains the handbag's coordinates xs, ys. The second line contains number n (1 ≤ n ≤ 24) — the amount of objects the girl has. The following n lines contain the objects' coordinates. All the coordinates do not exceed 100 in absolute value. All the given positions are different. All the numbers are integer.

Output

In the first line output the only number — the minimum time the girl needs to put the objects into her handbag.

In the second line output the possible optimum way for Lena. Each object in the input is described by its index number (from 1 to n), the handbag's point is described by number 0. The path should start and end in the handbag's point. If there are several optimal paths, print any of them.

Example
Input
0 0
2
1 1
-1 1
Output
8
0 1 2 0 
Input
1 1
3
4 3
3 4
0 0
Output
320 1 2 0 3 0 

 【题解】 这题就不多说了,说实话我自己还不太会,照着大佬博客敲的,敲过还是不太懂,先放这慢慢理解,给自己提个醒。

 【代码】

 

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#include<vector>
using namespace std;
#define mem(a,b) memset((a),(b),sizeof(a))
const int maxs=1<<24;
const int inf=0x3f3f3f3f;
int n,x[25],y[25];
int vis[maxs];
int dp[maxs];
int mapp[25][25];
int d[maxs];
vector<int>ans;

void Init()
{
    mem(dp,0x3f);
    mem(vis,-1);
    mem(d,0);
}

void out(int s)
{
    if(~vis[s])
        out(vis[s]);
    if(d[s])
    {
        if(d[s]>=100)
            ans.push_back(d[s]/100);
        ans.push_back(d[s]%100);
        ans.push_back(0);
    }
}

int main()
{
    while(~scanf("%d%d",&x[0],&y[0]))
    {
        Init();
        scanf("%d",&n);
        for(int i=1;i<=n;++i)
            scanf("%d%d",&x[i],&y[i]);
        int M=1<<n;
        for(int i=0;i<=n;++i)
        {
            for(int j=i;j<=n;++j)
                mapp[i][j]=mapp[j][i]=(x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]);
        }
        dp[0]=0;
        for(int i=0;i<M;++i)
        {
            if(dp[i]<inf)
            {
                for(int j=0;j<n;++j)
                {
                    if(!(i&(1<<j)))//位运算  第j个物品还没用过
                    {
                        int next_i=i|(1<<j),next_val=dp[i]+2*mapp[0][j+1];
                        if(next_val<dp[next_i])
                        {
                            dp[next_i]=next_val;
                            vis[next_i]=i;
                            d[next_i]=j+1;
                        }
                        for(int k=j+1;k<n;++k)
                            if(!(i&(1<<k)))
                        {
                            int next_i=i|(1<<j)|(1<<k) ,next_val= dp[i]+mapp[j+1][k+1]+mapp[0][j+1]+mapp[0][k+1];
                            if(next_val<dp[next_i])
                            {
                                dp[next_i]=next_val;
                                vis[next_i]=i;
                                d[next_i]=(j+1)*100+(k+1);
                            }
                        }
                        break;
                    }
                }
            }
        }
        printf("%d\n",dp[M-1]);
        out(M-1);
        printf("0 ");
        for(int i=0;i<ans.size();++i)
        {
            printf("%d%c",ans[i], i==ans.size()-1?'\n':' ');
        }
    }
    return 0;
}

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值