\(hash\)+容斥

但是看到这个令人愉快的数据范围还是直接枚举子集吧

首先我们发现\(6\)这个东西简直是小的可怜,复杂度里肯定有\(2^6\)的

于是我们可以直接先枚举子集,把所有状态的对应相等的位置有多少对搞出来

因为一个答案在一个集合里被算过那么就一定在子集里被算过

之后我们倒着枚举子集,顺便容斥就好了

代码

#include<iostream>
#include<cstring>
#include<cstdio>
#include<tr1/unordered_map>
#define re register
#define maxn 100005
#define uint unsigned long long
#define pt putchar(1)
#define lowbit(x) ((x)&(-x))
#define LL long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
using namespace std::tr1;
inline int read()
{
char c=getchar();
int x=0;
while(c<'0'||c>'9') c=getchar();
while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-48,c=getchar();
return x;
}
int a[maxn][7];
int n,m;
int dp[65];
LL ans;
uint base=666233;
inline int cnt(int x) {int tot=0;while(x) tot++,x-=lowbit(x);return tot;}
unordered_map<uint,int> ma;
inline int count(int x)
{
int tot=0,ct=0;int b[10];
while(x){ct++;if(x&1) b[++tot]=ct;x>>=1;}
ma.clear();
int now=0;
for(re int i=1;i<=n;i++)
{
uint k=0;
for(re int j=1;j<=tot;j++)
k=k*base+a[i][b[j]];
now+=ma[k];
ma[k]++;
}
return now;
}
int main()
{
n=read();m=read();
for(re int i=1;i<=n;i++)
for(re int j=1;j<=6;j++) a[i][j]=read();
int N=(1<<6);N--;
for(re int i=0;i<=N;i++) dp[i]=count(i);
for(re int i=N;i;i--)
for(re int t=i;t;t=(t-1)&i) dp[i^t]-=dp[i];
for(re int i=0;i<=N;i++)
if(cnt(i)==m) ans+=dp[i];
printf("%lld\n",ans);
return 0;
}

【[SDOI2013]泉】的更多相关文章

  1. [SDOI2013]泉

    题目描述 作为光荣的济南泉历史研究小组中的一员,铭铭收集了历史上x个不同年份时不同泉区的水流指数,这个指数是一个小于. 2^30的非负整数.第i个年份时六个泉区的泉水流量指数分别为 A(i,l),A( ...

  2. [SDOI2013]泉(容斥)

    /* 容斥加上哈希 首先我们可以2 ^ 6枚举相同情况, 然后对于这些确定的位置哈希一下统计方案数 这样我们就统计出了这些不同方案的情况, 然后容斥一下就好了 */ #include<cstdi ...

  3. 题解 洛谷 P3298 【[SDOI2013]泉】

    考虑到年份数很小,只有 \(6\),所以可以 \(2^6\) 来枚举子集,确定流量指数对应相同的位置,然后通过哈希和排序来计算相同的方案数. 但是这样计算出的是大于等于子集元素个数的方案数,所以还需要 ...

  4. Noip前的大抱佛脚----赛前任务

    赛前任务 tags:任务清单 前言 现在xzy太弱了,而且他最近越来越弱了,天天被爆踩,天天被爆踩 题单不会在作业部落发布,所以可(yi)能(ding)会不及时更新 省选前的练习莫名其妙地成为了Noi ...

  5. NOIP前的刷题记录

    因为这几天要加油,懒得每篇都来写题解了,就这里记录一下加上一句话题解好了 P4071 [SDOI2016]排列计数   组合数+错排 loj 6217 扑克牌 暴力背包 P2511 [HAOI2008 ...

  6. 【JZOJ3295】【SDOI2013】泉(spring)

    ╰( ̄▽ ̄)╭ 济南市"泉历史研究小组"依据济南特有的泉脉关系将济南的泉水分为六个区域,分别是市中区.历下区.天桥区.槐荫区.历城区.长清区. 作为光荣的济南泉历史研究小组中的一员 ...

  7. BZOJ3130: [Sdoi2013]费用流[最大流 实数二分]

    3130: [Sdoi2013]费用流 Time Limit: 10 Sec  Memory Limit: 128 MBSec  Special JudgeSubmit: 960  Solved: 5 ...

  8. Centos 下安装 文泉驿 字体 Odoo

    刚装完centos下的odoo的字体 文泉驿 ,一万头草泥马呼啸而过.....劝君如非必要,千万别再centos下折腾odoo..... 正题,文泉驿官网 只提供 deb包和源码包的字体安装 ,想在c ...

  9. freebsd|odoo - 为odoo报表 安装文泉译中文字体

    来源: Odoo8.0中使用文泉译中文字体         http://blog.csdn.net/vnsoft/article/details/40056935 FreeBSD wkhtmltop ...

随机推荐

  1. Timer控件

    Timer控件是定期引发事件的控件,时间间隔的长度由interval属性定义,其值以毫秒为单位吗,若启用了该组件,则每个事件间隔引发一个Tick事件,Timer组件的主要方法包括start和stop, ...

  2. MySQL---4、语句规范

    1.命名规范 (1)库名.表名.(按现在的规范类似; PromoHayaoRecord),数据库名使用小写,字段名必须使用小写字母,并采用下划线分割.关键字与函数名称全部大写.(2)库名.表名.字段名 ...

  3. string转xml

    //string转xml XmlDocument doc = new XmlDocument(); doc.LoadXml(temp); XmlNodeList nodes1 = doc.GetEle ...

  4. mac下如何找到hosts文件(转)

    打开Finder 在菜单中选择[前往][前往文件夹] 或使用快捷键Command+Shift+G        进入跳转路径,输入: /private/etc/ 点击[前往]      即可找到hos ...

  5. java web工程web.xml介绍

    转载自:http://blog.csdn.net/believejava/article/details/43229361 Web.xml详解: 1.web.xml加载过程(步骤) 首先简单讲一下,w ...

  6. pollard_rho 算法进行质因数分解

    //************************************************ //pollard_rho 算法进行质因数分解 //*********************** ...

  7. sql:PostgreSQL9.3 Using RETURNS TABLE vs. OUT parameters

    http://www.postgresonline.com/journal/archives/201-Using-RETURNS-TABLE-vs.-OUT-parameters.html http: ...

  8. arguments 参数

    下面要写的是知识梳理的一个案例: 写一个求和的方法sumFn,不管传递的参数有什么,都能将最终的和算出来,并且返回给函数外部使用.(要求:一个参数都不传默认结果为0,对于传递的非正常数字的参数不与累加 ...

  9. [SD2015]序列统计——solution

    http://www.lydsy.com/JudgeOnline/problem.php?id=3992 很容易得出DP方程: f[i][c]=f[i-1][a]*f[1][b]① 其中a*b%M=c ...

  10. 旋转/非旋转treap的简单操作

    treap(树堆) 是在二叉搜索树的基础上,通过维护随机附加域,使其满足堆性质,从而使树相对平衡的二叉树: 为什么可以这样呢? 因为在维护堆的时候可以同时保证搜索树的性质: (比如当一棵树的一个域满足 ...