HDU4352 XHXJ's LIS(LIS 状压)
程序员文章站
2022-04-09 18:40:23
题意 "题目链接" Sol 刚开始的思路是$f[i][j]$表示到第$i$位,LIS长度为$j$的方案。 然而发现根本不能转移,除非知道了之前的状态然后重新dp一遍。。 题解,,,挺暴力的把,直接把求LIS过程中的单调栈当成一个状态压进去了。。 自己真是不长记性,明明已经被这个单调栈坑过一次了。。 ......
题意
sol
刚开始的思路是$f[i][j]$表示到第$i$位,lis长度为$j$的方案。
然而发现根本不能转移,除非知道了之前的状态然后重新dp一遍。。
题解,,,挺暴力的把,直接把求lis过程中的单调栈当成一个状态压进去了。。
自己真是不长记性,明明已经被这个单调栈坑过一次了。。
考虑到$k$非常小,于是直接对$k$进行状压
设$f[i][sta][j]$表示长度为$i$,单调栈内状态为$sta$, lis长度为$k$的方案数
最后一维如果是单组数据的话是不必要的。但是考虑到时多组数据,就一起加进来吧。
转移的时候枚举一下这一位放了什么,然后暴力的改一下lis的状态。
#include<bits/stdc++.h> #define ll long long #define int long long using namespace std; const int maxn = 1e5 + 10; ll t, l, r, k; int f[64][1 << 10][11];//长度为i,lis为k的方案数 int a[maxn], num; int change(int s, int x) {//向状态s中加入一个数x for(int i = x; i <= 9; i++) if(s & (1 << i)) {s ^= (1 << i); break;} return s | (1 << x); } int dfs(int now, int lim, int lead, int s) { // printf("%d %d %d %d\n", now, lim, lead, s); if(!now) return (__builtin_popcount(s) == k); if(!lim && f[now][s][k] != -1) return f[now][s][k]; int ans = 0; for(int i = 0; i <= (lim ? a[now] : 9); i++) ans += dfs(now - 1, lim && i == a[now], lead && (!i), (lead && (!i)) ? 0 : change(s, i)); if(!lim) f[now][s][k] = ans; return ans; } ll solve(ll x) { num = 0; while(x) a[++num] = x % 10, x /= 10; // cout << num << endl; dfs(num, 1, 1, 0); } main() { memset(f, -1, sizeof(f)); cin >> t; for(int i = 1; i <= t; i++) { cin >> l >> r >> k; printf("case #%d: ", i); cout << solve(r) - solve(l - 1) << endl; } return 0; }
上一篇: cf314E. Sereja and Squares(dp)
下一篇: word